Oil and gas reside in reservoirs within peritidal and shallow subtidal lagoonal carbonate sediments across the globe. This is a zone of facies heterogeneity, controlled by changes in depositional energy, water depth, clastic influx, and evapotranspiration. Close proximity to evaporitic brine pools means that it is also an environment with the potential for dolomitization during shallow burial. As a result, the original pore system of carbonate sediment can become drastically altered prior to burial, such that reservoir properties may not be predictable from facies models alone. The Miocene Santanyí Limestone Formation, Mallorca, Spain, is well exposed and has undergone minimal burial and therefore presents an excellent opportunity to integrate sedimentology, facies architecture, and diagenesis to determine how porosity evolves within individual facies in the shallow subsurface. From here, the impact on pore type, pore volume, pore connectivity, and petrophysical anisotropy can be assessed. The Santanyí Limestone consists of pale mudstones and wackestones, rooted wacke-packstones, stratiform laminites, and skeletal and oolitic, cross-bedded grainstone. Thin-section analysis reveals a paragenetic pathway of grain micritization, followed by dissolution of aragonite, possibly by meteoric fluids associated with karstification. Subsequently, the unit underwent fracturing, compaction, recrystallization, cementation, dolomitization, and matrix dissolution to form vugs. Petrophysical analyses of 2.54-cm-diameter plugs indicate that these complex diagenetic pathways created petrophysical anisotropy [mean horizontal permeability (Kh)/vertical permeability (Kv) of whole formation = 3.4] and that measured parameters cannot be related directly to either geological facies or pore type. Instead, petrophysical data can be grouped according to the diagenetic pathways that were followed after deposition. The best reservoir quality (i.e., typical porosity 15 to >40% and permeability >100 mD) is associated with pale mudstones, stratiform laminites, and skeletal and oolitic grainstone that have undergone pervasive recrystallization or dolomitization. These rocks have the some of the lowest formation resistivity factor (FRF) values (<200) and thus the simplest pore system. The poorest reservoir properties ( k <10 mD) occur in mudstones and wackestones that have not been recrystallized and, hence, are dominated by a simple network of micropores (FRF <101). Skeletal and oolitic grainstones and rooted and brecciated wacke-packstones that have undergone some cementation and partial recrystallization have moderate reservoir properties and a high FRF (>>1000), reflecting a complex pore system of biomolds, vugs, and microporosity. Consequently, reservoir properties can be predicted based on their primary rock properties and the diagenetic pathway that they followed after deposition.

Abstract Fault-related dolomite subsurface reservoirs are formed from fluid circulation that results in significant transformation of the reservoir properties. The geometry and internal organization of such dolomitic reservoirs remain difficult to image with seismics alone. A multi-scale approach is essential to understand and predict the diagenetic processes that control the exact 3D morphology of the dolomite with spatial precision and true dimensions, and consequently the reservoir properties. In this context, we propose an analytical workflow including field work, LIDAR scanning and numerical geology applied to dolomite outcrops in Mesozoic carbonates (SE France). The exposed dolomite-limestone contact exhibits sinuous, irregular and convolute shapes, which are either fault-parallel, bedding-parallel or chaotic. To characterize this complex distribution, we performed LIDAR scanning on 500 m x 150 m cliffs and road cuts with 4.5 cm to 1–1.5 cm average point spacing. The cloud is composed of 22 millions points comprising X, Y, Z, intensity, red, green, and blue attributes. Digitization of the limestone-dolomite boundary was performed in RiscanPro and GOCAD environments, for extracting the true 3D geometry of the dolomite body for further geostatistical and 3D facies modelling. This approach captures the large-scale geometry of the dolomite bodies. However, single RGB or intensity properties do not unequivocally reproduce small-scale (below ∼ 1 m) heterogeneities of the late diagenetic dolomite. Color changes induced by weathering or climatic conditions are of the same size range as the small-scale heterogeneities, thus they are not unique to allow automated tracking on the point set. As a result, the workflow remains time-consuming, and further work is needed to allow calibration of the LIDAR data points with mineralogy.