The Eocene El Garia Formation in the offshore Hasdrubal Field was originally a nummulitic limestone in which subsequent burial dolomitization has significantly enhanced permeability. Identification of the reservoir's petrophysical property distributions requires knowledge of the spatial extent of its dolomitization, in turn requiring understanding of the processes that caused the dolomitization. Some of this understanding can be derived from measurements but others need to be simulated. In this study, the former are used as guides and we focus on the latter, evaluating the character of the dolomitizing fluid's movement and temperature patterns by using basin modelling to develop heat-flux simulations to represent the time of dolomitization. Basin modelling reconstructs the region's geology at the time of dolomitization, while heat-flux simulations recreate the appropriate conductive and convective heat and mass transport through these systems. Potential key drivers are rock mass and fault-zone permeability, and the position and shape of any salt domes. The results suggest that salt dome shape and position is the dominant control, the salt dome localizing convective systems which also use convenient faults so that hotter upwelling fluids pass through the Hasdrubal reservoir and are instrumental in the development of burial dolomitization.@en

We have investigated the role of high-frequency cycles (HFCs) in carbonates during water flooding. These are normally sub-seismic and therefore below grid resolution and must be upscaled, both in terms of single- and multiphase flow behaviour, to predict time to water breakthrough, recovery factors, and the habitat of remaining oil properly. Particularly the latter could be important when designing appropriate enhanced oil recovery schemes that target the remaining oil. The importance of selecting a high vertical resolution grid has been demonstrated for HFCs with continuous petrophysical log gradients and discontinuities, and the consequence these high permeability ranges have on waterflood velocity has been presented: capturing all of a continuous gradient requires increasing the vertical resolution so as to predict a fast enough water front, and a realistic distribution of the residual oil within the HFCs. In addition, the importance of selective diagenesis has been investigated by considering the effect of early-stage diagenesis at the top of HFCs in end-member Greenhouse and Icehouse climates. HFCs are very common in carbonates, even if they are nested within larger-scale heterogeneous geobodies, and the value of this investigation is in its application to real reservoir models.@en

Porosity and permeability of carbonate sediments evolve markedly with burial depth, reflecting the combined effects of mechanical compaction, chemical compaction, dissolution and cementation. While trends in porosity change with depth can be qualified, the evolution of permeability remains problematic. Here, we create a theoretical series of 2D images of major pore-occluding and pore enhancing diagenetic processes linked to their depth of occurrence. These images were then used to create 3D pore architecture models using Markov Chain Monte Carlo simulation, from which pore network were extracted to obtain multiphase fluid flow properties. The modelled porosity and permeability evolution from three different diagenetic pathways display several tipping points where the decrease in permeability is significantly larger than the associated drop in porosity. Such diagenetic pathway models can provide constraints on the predicted behaviour of carbonates during burial and/or uplift scenarios.@en

Basin modelling tools are widely used to predict reservoir charging but, unfortunately, their application to carbonate reservoirs is far from straightforward. For carbonates, the relationships between porosity/ effective stress/depth, must be addressed before a reliable basin model can be produced. We propose using a synthetic diagenesis approach, where Pore Architecture Models (PAMs) are used to generate representative porosity values for the different diagenetic (backstripped) stages that are recognised from cores or thin sections. The synthetic diagenesis approach is then used to inform more conventional basin modelling. The viability of this method is the topic of this contribution.@en