Petroleum Geoscience - Volume 17, Issue 3, 2011

The carbonate-dominated Albian to Turonian succession is one of the key petroleum systems of the Arabian Plate. It is dominated by shallow water platform carbonates that were deposited in a shallow epeiric sea on the margins of Neotethys. In general, the reservoirs in this succession have high porosities but exhibit heterogeneous permeabilities. This study reviews published data for the region and attempts to unravel the key diagenetic controls on the porosity and permeability of the reservoirs. The results demonstrate that a spectrum of diagenetic processes created highly heterogeneous multimodal pore networks. Intense boring and micritization of skeletal allochems, differential cementation of a pervasive burrow network and preferential dissolution of aragonitic skeletal allochems are ubiquitous. Locally, particularly on the northern and eastern Arabian Plate, deep-penetrating epikarst can be tied to a differential response to global sea level fluctuation and local tectonism. The development of a peripheral bulge in late Cenomanian–Turonian times, halokinesis, localized influx of channelized clastic material and sub-regional climatic variability contributed to a heterogeneous pattern of meteoric diagenesis across the Arabian Plate. The succession was then buried to up to 10 km during the Late Cretaceous–Tertiary. Where deep-penetrating fault systems were reactivated by Alpine tectonism, flushing by hydrothermal brines resulted in highly localized patterns of hydrothermal dolomitization and leaching, associated with hydrocarbon emplacement.

Investigation of core and outcrop samples of the Cenomanian Natih-B Member (North Oman) indicates that the different lithofacies present experienced rather different early diagenesis shortly after deposition. Transmitted-light, cathodoluminescence and backscattered scanning-electron microscopy, as well as stable-isotopic, X-ray diffraction and total organic carbon (TOC) analyses were employed to delineate the major controls on the cyclic pattern of early diagenesis and hydrocarbon source potential.

The Natih-B intrashelf basinal carbonates are composed of pelagic sediments that exhibit high-frequency cyclicity marked by decimetre-thick lithofacies alternations, mainly between: Lithofacies A compacted, partially bioturbated, skeletal, organic carbon-rich mudstone to wackestone; and Lithofacies B uncompacted, extensively bioturbated, skeletal, sparry-calcite rich wackestone to packstone. Individual units are composed variously of authigenic and biogenic calcite (58.1–97.6%, average 78.5%) and organic carbon (0.3–13.7% TOC, average 3.6%), together with minor quartz, clay, pyrite, dolomite and phosphatic material (fish debris). Lithofacies A contains relatively more organic carbon, clay, pyrite and dolomite than Lithofacies B and constitutes an excellent source rock. Diagenetic textures of Lithofacies A are dominated by compactional deformation of burrow fabrics, faecal pellets and solution seams, in addition to zoned/bright luminescent, non-ferroan sparry and isopachous calcite cement in and around uncompacted foraminifer tests, in an uncemented matrix. In contrast, Lithofacies B does not show any signs of compaction other than microstylolites and is dominated by zoned/dull luminescent, non-ferroan calcite microspar replacement, in addition to pore-filling, predominantly dull-luminescent, non-ferroan, sparry calcite cement. Moreover, Lithofacies B shows evidence of isopachous and meniscus-style cementation, together with geopetal structures and mictritic peloids. Stable-isotopic compositions of both lithofacies were determined from whole-rock samples (δ13C = −0.9 to +0.9‰, average +0.3‰; δ18O = −5.6 to −3.7‰, average −4.8‰) and sparry calcite (both cement and matrix) subsamples (δ13C = −0.6 to +1.2‰, average +0.6‰; δ18O = −5.7 to −3.7‰, average −4.3‰); all results being relative to Vienna Pee Dee Belemnite.

These petrographic and isotopic characteristics suggest that the Natih-B abundant calcite cements and replacements were precipitated early, prior to compaction, mainly from ‘normal’ (open, oxic) seawater at slightly elevated depositional temperatures. Some of the slightly negative δ13C values, however, indicate an addition of isotopically light carbon, probably derived from organic-matter oxidation by organisms living in marine pore waters. Based on evidence of extensive seafloor bioturbation and cementation, and their position within the depositional succession, the tops of Lithofacies B (wackestones to packstones) are interpreted as ‘discontinuity surfaces’ that cap shallowing-upward, fifth-order cycles, formed as a function of sediment starvation and increased bottom-current activity during relative sea-level stillstand/turnaround. In contrast, Lithofacies A (mudstones to wackestones) is believed to reflect high organic production coupled with high sedimentation rate and rapid burial. These conditions limited total infaunal colonization and extensive calcite precipitation, and preserved organic matter together with some escape burrows and in-place fauna, suggesting episodic sediment influx when more accommodation was available and seafloor diagenesis was minimized during relative sea-level rises. The relatively higher amounts of pyrite and dolomite in Lithofacies A likely indicate organic-matter degradation by bacterial sulphate reduction in anoxic pore waters during shallow burial.

The Simply Folded Belt of the Zagros Mountains, Iran, is a spectacularly well-exposed example of a foreland fold and thrust belt. A regional analysis of the Cenomanian–Coniacian Sarvak and Ilam Formations, exposed in the southern Lurestan Province, is presented as a case study for sub-seismic fracture development in this type of compressive setting. The area is characterized by gentle to tight anticlines and synclines parallel to the NW–SE trend of the belt. In the Lurestan Province, the Cenomanian–Coniacian interval is exposed in the core of most of the outcropping anticlines. Fold style is intimately related to both vertical and lateral facies distribution.

Geometry, kinematics and timing of sub-seismic fractures were characterized through extensive fieldwork, interpretation of orthorectified QuickBird imagery and interpretation of 3D photorealistic models derived from LiDAR. Data were collected from 12 anticlines covering an area of approximately 150 × 200 km. Key outcrops for understanding fracture geometry, kinematics and timing are presented.

Field observations and interpretation of QuickBird and 3D photorealistic models reveal the complexity of fracture geometry and timing. Fractures record pre-, syn- and post-folding stages of deformation. Pre-folding structures include synsedimentary normal faults, and subsequent small-scale thrusts, systematic veins and stylolites. During folding, pre-existing fracture planes were re-activated and through-going fractures and reverse faults developed. Strike-slip faults typically postdate pre- and syn-folding structures and are probably related to the late stages of fold tightening. All structures are geometrically and kinematically consistent with the trend of the Arabian passive margin and its subsequent tectonic inversion.

Outcrop data derived from fieldwork, remote sensing satellite and LiDAR-derived 3D models were used to build an integrated dual porosity-permeability static reservoir model which captured stratigraphic, diagenetic and structural heterogeneities. The study focuses upon the Mishrif-Mauddud/Sarvak interval, one of the most prolific reservoir units in the Middle East. The study area is exceptionally well exposed in deep gorges which cut transversally across anticlines of the Simply Folded Belt of the Zagros Mountains. The outcrops reveal volumetrically significant dolomitization of the latest Albian to Turonian carbonates of the Lower and Upper Sarvak formations. Three different dolomite bodies, which are spatially connected and genetically linked to the same fluid flow event, were recognized and mapped: (1) a thick dolomite body replacing the Lower Sarvak and forming a massive dolomite core; (2) horizontally extensive stratabound dolomite bodies (sheets), emanating laterally from the massive dolomite; and (3) vertically elongated dolomite pipes, rooted in the massive dolomite and typically replacing slope facies of the Upper Sarvak Formation.

The widespread development of tight, non-planar dolomite textures (a typical feature of high temperature dolomitization) drastically reduces the reservoir potential of the dolomitized geobody in the study area. In particular, this is present in the massive dolomitized body. Vuggy porosity seems to increase porosity only locally and to a limited extent, developing a non-connected pore network. The dominant porous dolomite textures are more abundant in the peripheral part of the geobody (dolomite sheets), where they are strongly controlled by precursor facies and diagenesis. Three main dolomite pore types were identified (intercrystalline, interparticle and mouldic), linked to the depositional environment of the precursor limestone. These pore types were used for petrophysical modelling. The approach adopted in this study allowed the distribution of rock properties in the dolomitized geobody to mimic the main depositional facies architecture.

The study area is characterized by a simple fracture network. Two main fracture sets and two major sets of conjugate normal faults were recognized in the field and mapped on 3D virtual outcrop data. Non-stratabound fracture density varied according to stratigraphic unit and/or dolomite body type (pipes/massive/sheets), showing a general increase from precursor limestone to dolomite. Fracture density also varied according to distance from faults (fault damage zone). This was particularly true in the limestone. The data also showed a prominent increase in fracture height from limestone to dolomite bodies, indicating that the dolomitized geobodies are likely sites for high production and early water breakthrough.

In Jabal Madmar in the Sultanate of Oman, Cretaceous epeiric carbonate platform architectures were characterized by employing a digital outcrop modelling workflow. A framework model for Natih Sequence I (Natih E member) was established, which embeds a meticulously studied platform-top incision and shoal complex. Outcrop-scale clinoforms are recognized in these shoals by hectometre-scale (100 m long) medium to high-angle (1–5°) inclined stratal surfaces comprising texture-based facies transitions. These clinoforms are usually beneath the resolution of seismic data and as such are not easily recognized and correlated between wells. Geologically realistic clinoform models were built using a well-defined stratigraphic model that incorporated inclined surfaces in the model grid and if available, data on lateral facies transitions. Waterflood simulations demonstrated improved sweep efficiency in these models. In contrast, simple models without clinoform heterogeneities resulted in less efficient piston-like patterns of sweep. The study presented in this paper demonstrates an outcome contrary to previous studies, as in this study, barriers to flow are absent. Complex clinoform models must be considered in reservoir modelling workflows to correctly derive static and dynamic rock properties. This is because outcrop-scale clinoforms have a potential impact on reservoir behaviour under secondary and tertiary recovery mechanisms.