GeoArabia, 2013, v. 18, no. 3, p. 135-178 Gulf PetroLink, Bahrain

Mid-Permian Khuff Sequence KS6: Paleorelief-influenced facies and sequence patterns in the Lower Khuff time-equivalent strata, Oman Mountains, Sultanate of Oman

Daniel Bendias, Bastian Koehrer, Michael Obermaier and Thomas Aigner

ABSTRACT

Khuff Sequence KS6 was studied in the Al Jabal al-Akhdar, Oman Mountains, in an area of 30 x 50 square kilometers by means of detailed sedimentological logging of five time-equivalent outcrop sections of the Saiq Formation. KS6 represents one transgressive-regressive, third-order sequence, and is composed of four facies associations each representing particular environments of (backshoal, , foreshoal and offshoal) with distinct sedimentological characteristics. Facies stack to form cycles and cycle sets that were used for correlation at a subregional scale and to reveal the KS6 stratigraphic architecture. During the initial phase of basin-fill, clastic (“Basal Saiq Clastics”) were deposited in paleolows above the “Sub-Saiq Unconformity”. In contrast to younger Upper Khuff sequences KS4 to KS1, the underlying paleorelief strongly affects the thickness and facies composition of KS6. The correlation strategy to follow paleolandscape surfaces using all available sedimentological, biostratigraphic and lithostratigraphic data resulted in a stratigraphic architecture with subtle shingle geometries.

Sequence KS6 shows a strong facies partitioning resulting in the necessity of two separate facies models for the transgressive (crinoidal ramp) versus regressive hemisequence (oolitic/peloidal carbonate ramp). This study revealed potential reservoir units in KS6, commonly regarded as non-reservoir in the subsurface of Oman and other parts of the Gulf region. The abundance, nature and lateral extent of reservoir facies strongly varies with stratigraphic position. In the transgressive part of KS6, crinoidal grainstones are concentrated around the margin of a gentle paleohigh. They might have the best reservoir potential, although early diagenetic cementation is common in most settings. Oolitic/peloidal grainstones in the upper regressive part have a much higher diagenetic reservoir potential and are laterally much more widespread. Thus, Khuff Sequence KS6 differs profoundly in its stratigraphic architecture from the more “layer-cake”-like KS4 to KS1 sequences. Facies and thickness patterns are controlled by a marked paleohigh to paleolow configuration, resulting from the antecedent uneven topography during the Neo-Tethyan syn-rift setting, in contrast to the post-rift setting with low tectonic activity during KS4 to KS1.

INTRODUCTION

The Khuff Formation represents one of the most important carbonate reservoir units across the Middle East (Figure 1). Six sequences, KS1 to KS6 in descending order, have been identified in Khuff- equivalent strata in the Oman Mountains (Figures 2 and 3; Koehrer et al., 2010). The Upper Khuff sequences KS1 and KS4 contain several reservoir zones and are well explored across the Gulf region. The lower sequences KS5 and KS6 are so far only very sparsely investigated, and fewer penetrations leave these units relatively underexplored.

This study on Sequence KS6 is part of a larger research project on Khuff time-equivalent strata in the Al Jabal al-Akhdar area (Oman Mountains) in the Sultanate of Oman (Figures 1 and 2). Its aim is to unravel the geometries and distribution of Khuff grainstones as potential reservoir bodies. Initially a one-dimensional (1-D) facies and sequence-stratigraphic framework was proposed by Koehrer et al. (2010) for the Khuff time-equivalent strata at a type locality on the Saiq Plateau. Subsequently,

135

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al.

LATE PERMIAN (CHANGHSINGIAN) PALEOGEOGRAPHY

45°E 50° Inner shelf 55° carbonates

IRAQ 30°

30°N N KUWAIT Middle and outer 0 250 shelf carbonates Kuh-i-Mand km Dalan

North Pars Kangan SAUDI ARABIA Karan IRAN Khursaniyah Abu Sa'fah Zagros South Berri Pars High Abqaiq BAHRAIN North Field Ghawar AbuAl Bukhoosh Nasr 25° Khurais Awali Zakum 25° QATAR Gulf of Oman Riyadh Hail Study Area Harmaliyah UAE Muscat Arzanah Umm Arabian Shaif Shield Yibal

Marginal-marine/ OMAN deltaic deposits Erosional limit 20°

20°

45° 50° 55°

Figure 1: Paleogeographic map of the Arabian Gulf area during the Late Permian showing location of the study area and major hydrocarbon fields within the Khuff reservoir (red) (modified from Maurer et al., 2009, after Ziegler, 2001). The study area (Al Jabal al-Akhdar) is located about 130 km west of Muscat.

Khuff grainstone geometries were documented on the Saiq Plateau from near-well-scale (< 2 km) by Zeller et al. (2011) to field-scale (< 10 km) by Koehrer et al. (2011), and to subregional scale (< 60 km) by Koehrer et al. (2012). These studies revealed an overall “layer-cake”-type geometry of shoal grainstone bodies in Upper Khuff sequences KS4 to KS1.

Only subtle lateral heterogeneity of facies associations is observed in Khuff sequences KS1 to KS4 on a subregional-scale (60 x 40 km) in the Al Jabal al-Akhdar by Koehrer et al. (2012). Cycle sets and sequences were found to be highly correlatable, pointing towards a rather uniform gross depositional environment and absence of significant tectonic activity in this area during the post-rift phase of the Neo-Tethys Ocean. Post-depositional deformations and structurally-related breccias and discontinuity surfaces are however described in sequences KS1 to KS4 from the Saih Hatat Window (Figure 2), paleogeographically located more proximally to the Arabian platform margin (e.g. Chauvet et al., 2009; Weidlich and Bernecker, 2011; Weidlich and Bernecker, 2012).

However, heterogeneities markedly increase in KS5 and KS6, possibly due to deposition during a syn-rift phase. Walz and Aigner (2012) and Walz et al. (2013) report significant thickness variations and apparent downlaps on cycle and cycle set scale, which are interpreted to result either from differential or initial topography. This paper focuses on the architecture of the grainstones in Khuff Sequence KS6 on a scale of 30 x 50 sq km. It aims at assessing the possible influence of paleotopography on potential reservoir facies distribution and continuity. It follows a systematic 1-D/2-D/3-D approach: after a detailed documentation of 1-D outcrop sections, 2-D stratigraphic cross-sections are constructed, leading to 3-D geocellular models.

136 136

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman

OUTCROP MAP, OMAN MOUNTAINS

57°E 57°30' 58°558°30' 9° Gulf of Oman Q Jabal N Bawshar 0 25 Muscat km Semail Wadi Ophiolite S Tt Aday 23°30'N J Wadi Hulw 23°30' Wadi Bani Awf PTr Mayh PTr Wadi Sahtan Jabal Jabal Abu Da’vd Tayin Pc Wadi Mistal J Wadi CO CO Mijlas Tt P Al Jabal al-Akhdar JK Semail Ja ba Jabal Pc Ophiolite l J Ab Aswad Wadi Bani Hajir ya J d

23° JK 23° Saiq Plateau Semail Ophiolite Tt

Tt

Jabal S Hamrat ad Duru H Range am m S ah Hawasina Q Jabal Nappes Safra 57° 57°30' 58°558°30' Q 9°

BAHRAIN 54° 58° (Q) Quaternary Kawr Group (Triassic – Cretaceous) 26° QATAR Gulf of Oman (Tt) Tertiary Al Aridh Group Abu (Triassic – Late Cretaceous) Dhabi Location UAE (J) Sahtan Group (Jurassic) Umar Group (Triassic – Cretaceous) Muscat (JK) Kahmah Group Hamrat Duru Group 22° OMAN 22° (end Jurassic – mid-Cretaceous) (Late Permian – Late Cretaceous) SAUDI Baid Formation Wasia Group (mid-Cretaceous) N ARABIA (Late Permian – Jurassic) 0 200 (PTr) Akhdar Group Muti Formation (mid-Late Cretaceous) km (Late Permian – Triassic) Arabian Sea (CO) Haima Group Aruma Group (end Cretaceous) (Cambrian – Ordovician) 18° 18° (S) Semail Ophiolite (P) upper Huqf Group YEMEN (Proterozoic – Cambrian) (mid-Late Cretaceous) 54° 58° Metamorphic sole (Pc) lower Huqf Group (Proterozoic) (mid-Late Cretaceous) Correlation lines (Figs. 28 to 30) Figure 2: Geological map of the Oman Mountains showing location of the studied sections. Note that the Saiq and Mahil formations of the Akhdar Group are shown together (PTr) (after Béchennec et al., 1993).

GEOLOGICAL SETTING

The study area is located in the Oman Mountains, about 130 km west of Muscat (Figures 1 and 2). The northern Oman region experienced several phases of long-term subaerial exposure between the Proterozoic and Early Permian, which resulted in times of non-deposition and (e.g. Forbes et al., 2010). In the Oman Mountains, a distinctive angular unconformity between Neo-Proterozoic and Middle Permian strata is preserved, which is here termed as the “Sub-Saiq Unconformity” (Figures 3 and 4).

The initial deposition of Khuff sediments started as a result of the drift of Cimmerian terranes away from Gondwana and the accompanying opening of the Neo-Tethys Ocean (e.g. Pillevuit, 1993; Stampfli and Borel, 2002). Subsequent flooding of the Arabian Plate resulted in the development of an epeiric carbonate ramp and the unconformable deposition of Khuff and equivalent sediments on folded Proterozoic strata (e.g. Glennie et al., 1974; Rabu et al., 1993; Searle, 2007).

137 137

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al.

According to the paleogeographic reconstructions of the Permian–Triassic by Konert et al. (2001a, b) the Arabian Plate drifted northwards, from 30° to 15° south. Sea-level oscillations of moderate amplitude and wavelength, and a transitional icehouse to greenhouse climate (Al-Jallal, 1995), created conditions very similar to those observed along the present-day Arabian Gulf coast (Strohmenger et al., 2006). According to paleogeographic reconstructions from Ziegler (2001), the deposition of Khuff strata took place on a shallow-marine to open-marine epeiric carbonate ramp (Figure 1).

Stratigraphic Framework

Permian–Triassic strata in the Oman Mountains (Al Jabal al-Akhdar) are subdivided into two formations by Glennie et al. (1973) (Figure 3):

Outcrop ethys) Stage Er a

MF S Marker Beds -Order Group Epoc h Period (Oman ) rd Member Lithology Reservoir 3 Formation Equivalen t Rift-stages Subsurface (Neo-T

Olenekian Sudair Middle 250.0 iassi c Earl y Mahil

Tr "Top Breccia" Mesozoic K1 Induan Tr20 Upper Lower KS1 Tr10 K2 252.2 to "Saiq/Mahil KS3

Boundary" Post-rif t Changhsingian P40 K3 "Coral Marker" (Dorashamian) "Microbial Marker 3" 254.2 Middle Lopingian Wuchiapingian KS4 K4 (Dzhulfian) P30 Akhdar 259.8 "Microbial Marker 2" f Khuf Upper Capitanian KS5 (Midian) Saiq "Chert Marker" Permian Paleozoic

265.1 "Microbial Marker 1" K5 P20 Lower "Muddy Marker" Syn-rift Wordian Guadalupian (Murgabian)

KS6 This study 268.8 This study Roadian

(Kubergandian) Lower Basal Saiq Clastics Pre-Permian Metamorphic Basement Clastics Figure 3: Chronostratigraphic and sequence-stratigraphic framework of Limestone the studied interval (Koehrer et al., 2010). Age dates from Baud and Dolomite Bernecker (2010). Geological time scale from Gradstein et al. (2012). Tentative correlation to MFS postulated by Sharland et al. (2004). Shale

138 138

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman

• Saiq Formation (Permian), time-equivalent to the Lower and Middle Khuff Formation (K7–K3 reservoir units), • Mahil Formation (Triassic), time-equivalent to the Upper Khuff (K2 and K1 reservoir units), Sudair and Jilh formations.

Baud and Richoz (2013) explained that some authors apply a different lithostratigraphic scheme, also attributed to Glennie et al. (1974), in which the Saiq Formation corresponds to the entire Permian– Lower Triassic Khuff Formation, and the Mahil Formation to the overlying Triassic formations. In the range of this research project (this paper as well as Koehrer et al., 2010, 2011, 2012; Pöppelreiter et al., 2011; Zeller et al., 2011; Obermaier et al., 2012; Walz et al., 2013), we refer to the lithostratigrahic Lower Saiq/Mahil Boundary of Glennie et al. (1974), defined in the type locality of the Saiq Formation on the Saiq Plateau.

In the Al Jabal al-Akhdar outcrops, the Khuff-equivalent strata are divided into six sequences: Khuff Sequence KS6 to Khuff Sequence KS1 from bottom to top (Koehrer et al., 2010). These sequences are further subdivided into fourth-order cycle sets and fifth-order cycles.

Based on the Lower Saiq/Mahil Boundary of Glennie et al., 1974, the following lithostratigraphic units can be identified in Khuff-equivalent strata in the Oman Mountains:

• Sub-Saiq Unconformity: The contact between Saiq Formation (KS6) and pre-Permian metamorphic basement strata, represented by the Mistal, Hajir and Mu’aydin formations, forms a rather spectacular angular unconformity in the Oman Mountains (Figure 4). • Lower Saiq Member: Conglomerates, sandstones, siltstones with paleosoils and ostracods (Rabu et al., 1986) transitioning into marine carbonates of the Upper Saiq Member. This unit might be a time equivalent to the lowermost Khuff in the subsurface or the Gharif Formation.

KHUFF SEQUENCE KS6, WADI BANI AWF

KS6 limestones

Sub-Saiq Unconformity

Basal Saiq Clastics

Proterozoic deposits

Figure 4: The angular Sub-Saiq Unconformity separates the pre-Permian basement from Permian Basal Saiq Clastics and the KS6 limestones. Outcrop coordinates: 0547837E; 2575791N.

139 139

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al.

• Upper Saiq Member: Shallow to open-marine carbonate deposits of an epeiric ramp. Above the first 120 m of limestones, the Upper Saiq Member is entirely dolomitized. It is time-equivalent of the Lower and Middle members of the Khuff Formation in subsurface Oman (Osterloff et al., 2004), and KS6 to basal KS2 (Koehrer et al., 2010). • Lower Mahil Member Koehrer et al. (2010, 2012): Shallow- to open-marine carbonate deposits of an epeiric ramp (entirely dolomite). Time-equivalent of the Upper Khuff Member in subsurface Oman (Osterloff et al., 2004) or middle part of Sequence KS2 and Sequence KS1 (Koehrer et al., 2010). Baud and Richoz (2013) propose that this unit be named the “Saiq unit C”.

Biostratigraphy

First biostratigraphic studies of KS6 deposits which have been carried out on the Saiq Plateau by Montenat et al. (1976) suggested the lower part of Sequence KS6 to be of a Middle Murgabian (Wordian) age due to the presence of Neoschwagerina schuberti. Data from basal strata in the deep- water sections of the Hawasina Basin indicate a Roadian (Henderson and Mei, 2003) or Wordian (Kozur and Wardlaw, 2010) age according to conodonts Hindeodus excavatus and Hindeodus wordensis. Recent biostratigraphic studies from Forke et al. (2012) report the presence of Verbeekina grabaui and primitive (Afghanella? cf. tereshkovae) species in the basal part of KS6 in the Al Jabal al-Akhdar area, which indicates a stratigraphic range from the base of the Murgabian (late Roadian?/Wordian) to the lower part of middle Murgabian (Leven, 1997, 2009).

METHODS AND DESCRIPTION OF OUTCROP SECTIONS

This study focuses on Sequence KS6 in the wadis on the northern rim of Al Jabal al-Akhdar region (Wadi Sahtan, Wadi Hajir, Wadi Bani Awf and Wadi Mistal) and one section from the southern flank (Saiq Plateau) (Figure 2). In all logged sections the lower part of KS6 consists of limestones with some meters of basal clastics at the base, while the upper interval of about 50 m is entirely dolomitized. Each outcrop, where sections have been logged, is located by UTM coordinates at its base and representative figures (Table 1).

Table 1 Location, coordinates and figures for all outcrop sections Location Longitude Latitude Figures Saiq Plateau 0572491E 2553176N Figures 6 and 7 Wadi Sahtan 0532062E 2580402N Figures 8 and 9 Wadi Bani Awf 0547837E 2575791N Figures 10 and 11 0553398E; 2570306N (coordinates are taken in Wadi Bani Wadi Bani Hajir/ Kharuz Hajir, where the upper 90 m of the section have been logged). Coordinates for lower part: 0555626E; 2570098N Figures 12 and 13 Wadi Mistal 0570359E 2574936N Figures 14 and 15

A standardized logging sheet was used to describe five sections of Sequence KS6 in detail. GPS coordinates were captured for all sections (Figure 2). The following properties were recorded and digitized in WellCAD 4.2: texture (Dunham, 1962), (biogenic and physical), lithology, components and spectral gamma-ray, integrated into the facies types and facies associations according to Koehrer et al. (2010). The color code used in WellCAD, correlation panels and 3-D models is presented in Figure 5.

More than 250 samples were taken for microfacies analyses in thin sections, which were carried out with a transmission light microscope. The integration of thin section analysis and field observations led to a highly reliable facies characterization in KS6 with calibrated facies logs and a detailed microfacies atlas (Figures 17 to 24).

A spectral gamma-ray (GR) survey was run in outcrops using a portable spectral GR spectrometer (model GS-512, manufactured by Geofyzika, Czech Republic). The spectrometer is equipped with a 3 x 3’’ NaI(TI) scintillation detector, which measures natural gamma-radiation of rocks. To detect overall GR trends usable for stratigraphic correlations and sequence interpretations, a sampling time

140 140

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman

Dunham Texture Lithofacies Types Lithology 2c (Burrowed to vertically rooted Mudstone Dolomite mud- to wackestone) Wackestone 3a (Microbial laminites) Limestone

Packstone 5a1 (Graded wacke- to mudstone) Clastics Pack- to Grainstone 5a2 (Graded pack- wackestone) Grainstone Cycle 5c1 (Peloidal-rich pack- grainstone) Rudstone Regressive Boundstone 5c2 (Bioclast-rich pack- grainstone) hemicycle

5d1 (Well sorted oolitic grainstone) Transgressive Lithofacies Association hemicycle LFA 4 (Backshoal) 5d2 (Well sorted peloidal grainstone) LFA 5 (Shoal) 5e (Skeletal floatstone) LFA 7 (Foreshoal) LFA 8 (Offshoal)

Figure 5: Color code for KS6 lithologies, texture, lithofacies types, lithofacies associations and cycles.

interval of 15 seconds every 0.5 m was found to be sufficient after test measurements of 180, 90, 30 and 15 seconds (Koehrer et al., 2010). The authors are aware that quantitative analysis cannot be carried out using 15 second measuring time. However it is sufficient to capture and reproduce clear GR trends. The concentration of each of the elements are automatically calculated by the instrument and displayed in ppm (U, Th), % (K) or counts per second (cps, total GR). Although spectral GR was recorded only total GR was used for correlation and cycle interpretation since trend curves of U, Th and K seem not reproducible with a measuring time interval of 15 seconds.

Five KS6 sections in the Oman Mountains were correlated on a scale of 30 x 50 km (Figure 2). All correlation scenarios use the top of KS6 as the datum since this surface is well traceable throughout the study area. Various correlation scenarios were tested, focussing on the gamma-ray logs, carbonate facies and cycles. The preferred correlation strategy was that correlation lines follow paleolandscape surfaces, i.e. elevation changes associated with depositional profile, integrating gamma-ray, facies, and sequence-stratigraphic interpretations.

Facies associations were used to create a 3-D facies model with standard industry software (Petrel). Several modeling techniques were applied based on different correlation strategies in order to provide a wide spectrum of possible reservoir geometries (cf. modeling section).

SEDIMENTOLOGICAL DEVELOPMENT

In general, Sequence KS6 is subdivided into four units from the base to the top: (1) Basal Saiq Clastics (clastics, Figures 16 to 18); (2) transgressive hemi-sequence (limestone, Figures 19 and 20), (3) maximum flooding interval (limestone, Figures 21 and 22), and (4) and regressive hemi-sequence (dolomite, Figures 23 and 24).”

Basal Saiq Clastics Description: In all outcrops where the base Khuff is exposed, the Sub-Saiq Unconformity is overlain by clastic deposits, 2–10 m thick, which are mixed upwards with an increasing amount of carbonate clasts (Figures 16 to 18). They are termed Basal Saiq Clastics and vary in thickness, composition and grainsize from outcrop to outcrop. On the Saiq Plateau they are composed of rooted siltstones with some ostracods (Rabu et al., 1986), whereas cross-bedded sandstones and conglomerates with imbrications dominate in the Wadi sections on the northern flank of the Oman Mountains. A variety of sedimentary structures, grain sizes and lithologies can be observed at the different outcrops of the Basal Saiq Clastics throughout the Oman Mountains.

141 141

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al.

KHUFF SEQUENCE KS6, SAIQ PLATEAU d Dunham Texture microbial laminites (meter) (top KS6) Gamma- Ray (API) ckeston e -order -order th rd 3 4 Rudstone Floatston e Boundstone Grainstone Packston e Wa Mudstone Thickness Lithology 10 40 180

170

160 2 cm

150 c massive mudstones (muddy marker, maximum flooding)

140

130

120

110

100 b crinoidal pack- to grainstones

90

80

70

60

50 a coral heads

40

30

20

10

0 Figure 6: Log of Khuff Sequence KS6 on the Saiq Plateau with typical facies: (a and b) Fossil-rich carbonate beds with coral heads and crinoidal pack- to grainstones. (c) At around 100 m above the base the interpreted maximum flooding interval is marked by fossil-poor, mud- and wackestones (Muddy Marker). (d) The upper, dolomitic portion gradually develops into oolitic and peloidal rich facies and is topped by microbial laminites forming the upper sequence boundary (Microbial Marker 1).

142 142

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman

LOWER KS6, SAIQ PLATEAU SECTION NEAR HAIL AL JEMEN

Mud- to wackestones

First KS6 carbonates, packstones 2 m Basal Khuff Clastics

Figure 7: Basal Saiq Clastics (mainly reddish siltstones) and the basal 50 meters of KS6 carbonates on the Saiq Plateau. Dashed line marks location of logged section.

Interpretation: The basal clastics are interbedded with fossil-rich Permian carbonates, which are interpreted to represent the initial Khuff transgression, and a transition from a terrestrial (lacustrine/ alluvial) to a shallow-marine environment. The fossil-rich carbonate clasts within the Basal Saiq Clastics contain typical Khuff fossils (Forke et al., 2012). They are interpreted as basal Khuff deposits and not pre-Khuff sediments (e.g. Gharif-equivalent strata). The Early Permian continental clastics of the Gharif Formation present in southwest interior of Oman pinch-out towards the Oman Mountains in the northeast (Blendinger et al., 1990).

Bioclastic Facies - Transgressive Hemisequence Description: The transgressive hemisequence of KS6 is characterized by bioclastic-rich beds (bioclastic wackestones, packstones, grainstones and floatstones (Figures 19 and 20). They contain a very diverse marine fauna (e.g. crinoids, fusulinids, coral heads, bivalve, brachiopod and gastropod shells) and beds are often graded or show bioturbation. Well-sorted and sometimes trough cross-bedded crinoidal pack- to grainstones are common especially in the lowermost part. They are interbedded with coral- rich units such as coral floatstones and coral framestones.

Interpretation: The very diverse marine fauna indicates fully open-marine conditions in a foreshoal environment. The cross-bedded, well-sorted crinoid columnals can be interpreted as high-energy bioclastic deposits which were probably deposited within range of wave energy in a shoal-like environment. The graded beds are likely to represent storm sheets (tempestites) forming meter-scale coarsening-upward cycles of a foreshoal and shoal margin environment. This interpretation attributes the transgressive hemisequence of KS6 to a shoal to foreshoal environment.

Deeper-water Facies - Maximum Flooding Interval Description: In this interval dark colored mud-rich carbonates become more frequent (Figures 21 and 22). Graded storm sheets with abundant skeletal debris interfinger with thinly bedded to very massive dark-blue mudstones, which contain diverse ichnofabrics such as Zoophycus burrows. In general this interval is characterized by a very low biodiversity. The main fossils are bivalve shells, gastropods and gymnocodiacean algae as well as some calcispheres; corals and crinoids are absent.

143 143

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al.

KHUFF SEQUENCE KS6, WADI SAHTAN e Dunham Texture Peloidal/oolitic grainstones and Top KS6 (meter) Gamma- Ray (API) ckeston e -order -order th rd 3 4 Thickness Rudstone Floatston e Boundstone Grainstone Packston e Wa Mudstone Lithology 10 40 170 Onset of 160 dolomitization

150 d Massive and thinly bedded mudstones 140 (maximum flooding)

130

120

110

2 cm 100

b c Bioclasts 90 Coral head

80

70

60

20 cm 2 cm 50 a KS6 limestones 40

30

20

10

0 Basal Saiq Clastics Figure 8: Log of Sequence KS6 in Wadi Sahtan with typical facies. (a) The contact of Proterozoic and Permian strata (Sub-Saiq Unconformity) is characterized by the erosive base of some 1−2 meter thick Sandstones sandstones. (b and c) The sandstones are Proterozoic overlain by increasingly fossil-rich carbonates with abundant coral heads and 2 m crinoids. (d) Around the maximum flooding interval massive and thinly bedded mudstones dominate whereas the upper dolomitic part (e) is completely dominated by peloidal and oolitic grainstones. The top of the KS6 is marked by microbial laminites. 20 cm

144 144

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman

KHUFF SEQUENCE KS6, WADI SAHTAN

Dunham Texture

(meter) Top KS6 Gamma- Ray (API) ckestone -order -order th rd Thickness Lithology 3 4 Rudstone Floatstone Boundstone Grainstone Packstone Wa Mudstone 040 170

160

150

140

130

Maximum Flooding Interval

120

110

100

90

80

70

60 10 m 50

40 Figure 9: Cliff showing most parts of Khuff KS6 Sequence in Wadi Sahtan. Note that in Wadi Sahtan the onset of 30 dolomitization (marked by yellow dashed line) does not appear to be parallel to bedding. 20

10 Interpretation: The increasing lime mudstone volume

0 indicates that with progressing transgression in the lower part of KS6, the overall deepening-upward trend from shallow-marine to open-marine carbonates continued. The lack of fauna, as well as the dark color, suggests less oxygenated waters and an azoic environment. Texture (mudstones, wackestones) and possible calcispheres, together with abundant gymnocodiacean algae (Permocalculus), indicate deep water far below the storm wave base (SWB) characterized by suspension settling, sediment starvation and muddy background sedimentation. This mudstone interval forms the so called “Muddy Marker” (Koehrer et al., 2010) and is interpreted as the maximum flooding interval of the lowest Khuff Sequence KS6 (Figure 3).

Oolitic/Peloidal Grainstones - Regressive Hemisequence Description: The upper part of KS6 is about 50 m thick and entirely dolomitized. Low biodiversity, a high percentage of ooids and peloids, and a lack of crinoids are observed (Figure 23 and 24). Fossil content is restricted to bivalve and brachiopod shells and few foraminifers. Massive, cross-bedded

145 145

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al.

KHUFF SEQUENCE KS6, WADI BANI AWF f g Dunham Texture Rooted mudstone at the

(meter) top of Khuff KS6 Gamma- Ray (API) ckestone -order -order th rd 3 4 Thickness Lithology Rudstone Floatstone Boundstone Grainstone Packstone Wa Mudstone 040

180

170 2 cm Side view 2 cm Top view

160 e Zoophycus burrows

150

140

130 5 cm

120 c Coral head in crinoidal d Crinoidal grainstones grainstone matrix

110

100

90 1 cm 1 cm

80 a b Mix of clastics and Basal Khuff carbonates Clastics

70

60 Proterozoic

50 10 cm

40 Figure 10: (a and b) Up to 10 m-thick unit of conglomerates form the Basal Saiq Clastics. They show an erosive base 30 and become upwards mixed with fossil-rich carbonates. (c and d) Especially the section in Wadi Bani Awf is in the lower part dominated by up to meter-thick crinoidal 20 grainstones (d). They are interbedded with coral-rich intervals which sometimes appear as coral framestones (c).

10 (e) Zoophycus burrows in mudstones are common around the maximum flooding interval and the uppermost 45 meters are again dominated by peloidal and oolitic 0 grainstones. (f and g) The top of the KS6 is marked by a rooted mudstone.

146 146

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman

and graded oolitic and peloidal grainstones are abundant. The KS6/KS5 sequence boundary (SB KS5) is placed on top of an up-to-50-cm-thick microbial laminite, which can laterally develop into rooted mudstone on a 10 km scale. This interval was termed “Microbial Marker 1” by Koehrer et al. (2010). Baud et al. (2012) described from the upper KS6 on the Saiq Plateau paleokarst features and ferrigenous crusts; however, in this study such features were not observed in the KS6 section on the Saiq Plateau.

KHUFF SEQUENCE KS6, WADI BANI AWF b

Top KS6

MFI

Top corals and crinoids 20 m

LOCATION OF LOGGED SEQUENCE KS6, WADI BANI AWF a Mudstones

Wacke- and packstones

Corals and crinoids

Figure 11: (a) Photo shows the logging location (marked by a dashed line) and the main facies units in Wadi Bani Awf. (b) Outcrop wall showing overview of KS6 Sequence in Wadi Bani Awf.

147 147

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al. KHUFF SEQUENCE KS6, WADI BANI HAJIR/KHARUZ f Dunham Texture Top KS6 “Microbial Marker” (meter) Gamma- Ray (API) ckestone -order -order th rd Thickness Lithology 3 4 Rudstone Floatstone Boundstone Grainstone Packstone Wa Mudstone 040

100

5 cm 90 e 80 Cross-bedded grainstones

70

60

50

10 cm

40 d Bioturbated mudstone

30

20

10

10 cm 0

a c Carbonates with Proterozoic clasts

Conglomerates 10 cm

b

Conglomerates

Proterozoic

20 cm 10 cm

Figure 12: (a, b and c) Meter-thick conglomerates showing an erosive base mix and intercalated with carbonates. (d) The fossil-rich transgressive lower part of the section is followed by a muddy interval (Muddy Marker) around the maximum flooding interval, (e) turning upwards into peloidal/oolitic grainstones. (f) Microbial laminites (Microbial Marker 1) form the top of the section. (Due to accessibility the lower 16 m of this composite section have been logged in Wadi Bani Kharus, whereas everything above was logged in Wadi Bani Hajir).

148 148

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman Proterozoic

adi Bani Kharus is pre served as a

Maximum Flooding Interval Flooding Maximum Basal Saiq Clastics Saiq Basal op KS6 ? T

ADI BANI KHARUS

m 10 KHUFF SEQUENCE KS6, W 40 Ra y (API) adi Bani Kharus. The contact of Proterozoic and Permian sediments in W

Gamma-

10

3 -order

rd

-order 4

th

Lithology

Mudstone

ckestone Wa Packstone

exture

Grainstone

Boundstone

Floatstone

Dunham T Rudstone

Thickness Thickness (meter) 0 90 80 70 60 50 40 30 20 10 spectacular angular unconformity (blue line marks location of maximum flooding interval). The base of Basal Saiq Clastics is ero sional. spectacular angular unconformity (blue line marks location of maximum flooding interval). The base Basal 100 Figure 13: Overview of the KS6 in W

149 149

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al.

KHUFF SEQUENCE KS6, WADI MISTAL c Dunham Texture Cross-bedded oolitic and peloidal grainstones (meter) Gamma- Ray

Lithology (API) ckestone -order -order th rd 3 4 Rudstone Floatstone Boundstone Grainstone Packstone Wa Mudstone Thickness 040

120

110

100 b Massive mudstone around maximum flooding

90

80

70

60 a Bioclastic packstones 50

40

30

20

10 ?

0 Figure 14: (a) The 130 meter thick KS6 in Wadi Mistal starts with bioclastic packstones and (b) some grainstones, develops into massive mudstones around the maximum flooding interval and (c) consists of cross-bedded peloidal/oolitic grainstones in the regressive part.

150 150

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman

KHUFF SEQUENCE KS6, WADI MISTAL a

Top KS6

l Maximum Flooding Interva 30 m

LOCATION OF LOGGED SEQUENCE KS6, WADI MISTAL b

Figure20 m 15: (a) Cliff with marked location of top KS6 and the zone of maximum flooding. (b) Location of the section in Wadi Mistal (dashed white line). Note that the contact between Proterozoic and Permian strata in Wadi Mistal is not exposed. Due to observations close to the logged section it can be assumed that only few meters at the base are missing.

151 151

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al.

LOWER KS6, SAIQ PLATEAU SECTION NEAR SAIQ VILLAGE

KS6 carbonates

Mix of carbonates and clastics

Proterozoic 2 m

Figure 16: Basal Saiq Clastics (Saiq Plateau). Proterozoic deposits (lower left of the picture) are overlain by a mix of carbonates and clastics (Basal Saiq Clastics).

Interpretation: Proximal storm deposits with some bioclasts but a generally low biodiversity, indicate high-energy conditions in a more restricted shallow-marine environment. The lack of crinoids, typical for regressive systems, can be observed in similar settings like in the Middle Triassic Muschelkalk in Germany (Palermo et al., 2010). Abundant cross-bedded peloidal and oolitic grainstones probably represent shoal bodies, which prograde over open-marine and foreshoal facies from the maximum flooding interval below. In general the regressive hemisequence can probably be placed into a backshoal to shoal environment with some influence from a rather proximal foreshoal.

Facies Models

In order to describe the development of the KS6 carbonate ramp two facies models are required. One for the transgressive part, a system with high biodiversity and no real shoal where the crinoidal grainstones form the most proximal facies (Figure 25). The other is for the regressive part with very typical shoal and backshoal facies (Figure 26). Thus Sequence KS6 shows a high degree of facies partitioning resulting in the need for separate facies models for the transgressive and regressive hemi-sequences.

CYCLE HIERARCHY

Facies types show vertical stacking patterns that indicate a hierarchical cyclicity. The terminology to describe this cycle hierarchy was adopted from Kerans and Tinker (1997). Accordingly, cycles are stacked to form cycle sets, which can be grouped into high-frequency and composite sequences, which in turn build super-sequences. In general three different cycle orders can be detected in Sequence KS6, 1–5 m thick cycles, 6–15 m thick cycle sets and the overall Sequence KS6. Recent biostratigraphic studies indicate that the timeframe of KS6 ranges from late Roadian to late Wordian (Forke et al., 2012).

152 152

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman

BASAL SAIQ CLASTICS

General features Texture Siltstones and quartz grains interbedded with small carbonate packstone beds

Mineralogy Siltstones, limestone Main Silt, bivalve shells, brachiopods, Components quartz grains

Environment Terrestial - shallow-marine

Interpretation Basal siltstones with rootlets indicate terrestrial/lacustrine environment. Interfingering small packstone beds with quartz grains represent fast transgression.

1 m

Basal siltstones First carbonate beds within siltstones

Rootlets

0.5 m 0.5 m

Figure 17: Basal Saiq Clastics on the Saiq Plateau. Near the village Hail al Jemen (lower left picture) and near Saiq village (lower right picture).

153 153

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al.

Siltstone Packstone

Organic material Ostracod

Quartz grain

Quartz grain

Organic material

Brachiopod shell

0.5 cm 0.5 cm

Packstone Packstone

Organic material

Quartz grains

Gastropod

Brachiopod Quartz grain Brachiopod shell

Organic material

Possible ostracod

1 cm 1 cm

Figure 18: Appearance of the Basal Saiq Clastics in thin sections from Saiq Plateau, Oman Mountains.

154 154

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman

BIOCLASTIC TRANSGRESSIVE FACIES

General Features Top KS6

Texture Bioclastic wackestones, packstones, grainstones and floatstones

Mineralogy Limestone MFI Main Crinoids, corals (rugose), Components foraminifers, gastropods, bivalves, skeletal debris

Environment Shoal to foreshoal

Interpretation High energy crinoidal above fair weather wave base (FWWB) and moderate energy bioclastic storm deposits above storm wave base (SWB)

20 m

Transgressive System

Graded storm beds and Coral blankets and Zoophycus muds and Crinoidal shoal bioclastic sheets crinoidal gardens azoic mudstones FWWB SWB Shoal Proximal foreshoal

Distal foreshoal

Offshoal

Crinoidal/Bioclastic Grainstone Bioclastic Packstone

Crinoid

Bivalve shell

Gastropod Coral

5 cm 2.5 cm

Figure 19: Bioclastic transgressive facies of Sequence KS6 with pictures of typical facies from outcrops and their interpreted position on the carbonate ramp.

155 155

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al.

Bioclastic Grainstone Rugose Coral

Crinoid

Foraminifer Bryozoa

1 cm 1 cm

Bioclastic Packstone Bioclastic Packstone

Foraminifer Foraminifer Crinoid

Gastropod

Likely brachiopod Shell (Bivalve)

Bivalve shell

Rugose coral

1 cm 1 cm

Figure 20: Bioclastic transgressive facies of the KS6 in thin sections from Wadi Bani Awf, Oman Mountains.

156 156

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman

DEEPER-WATER FACIES

General Features Texture Bioturbated mudstones Top KS6 and wackestones, azoic mudstones Mineralogy Limestone Main Foraminifers, gastropods, Maximum Components bivalves, skeletal debris flooding KS6 Environment Offshoal

Interpretation Low to very low energy sediments mostly below storm wave base (SWB). Lack of bioturbation in some intervals points to azoic environment.

20 m

Transgressive System Graded storm beds and Coral blankets and Zoophycus muds and Crinoidal shoal bioclastic sheets crinoidal gardens azoic mudstones FWWB SWB Shoal Proximal foreshoal

Distal foreshoal

Offshoal

Bioturbated Mudstone Azoic Mudstone

Zoophycus burrows

Stylolites

5 cm 2.5 cm Figure 21: Deeper-water facies around the interpreted maximum flooding of Sequence KS6 with pictures of typical facies from outcrops and their interpreted position on the carbonate ramp.

157 157

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al.

Storm Sheet Mudstone

Possible calcisphere

Mudstone overlain by distal storm sheet

1 cm 1 cm

Bioturbated Mudstone Bioclastic Wackstone

Bivalve shell Gastropod

Bioturbation

1 cm 1 cm

Figure 22: Deeper-water facies of Sequence KS6 in thin sections from Wadi Bani Awf, Oman Mountains.

158 158

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman

REGRESSIVE FACIES

General Features

Texture Cross-bedded oolitic/peloidal Top KS6 grainstones, packstones and wackestones

Mineralogy Dolomite Main Ooids, peloids, foraminifers, Maximum Components gastropods, bivalves, skeletal flooding debris, microbial laminites KS6

Environment Backshoal - foreshoal

Interpretation High energy oolitic/peloidal shoal sediments above FWWB and fossil poor fore- and backshoal deposits (above SWB) are indicating restricted conditions.

20 m

Regressive System Bioturbated mw/ X-bedded oolite/ Graded peloidal/bioclastic Mudstone microbial mats peloid shoal storm sheets (partly bioturbated) FWWB

SWB Backshoal Shoal Foreshoal Offshoal

Foraminifer Crinoid Gastropod Peloid Bivalve Shell Rugose Coral

Cross-bedded Oolitic/Peloidal Grainstone

Zoophycus burrows

5 cm Figure 23: The regressive facies below the interpreted top of Sequence KS6 with pictures of typical facies from outcrops and their interpreted position on the carbonate ramp.

159 159

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al.

Oolitic Grainstone Oolitic Grainstone

Foraminifer

Peloid Stylolites

Ooid Bioclast

Ooid

1 cm 1 cm

Mudstone and Bioclastic Packstone Grainstone and Packstone

Crinoid

Ooid filled burrow

Gastropod

Foraminifer

Stylolite

Stylolite

1 cm 1 cm

Figure 24: Regressive facies of Sequence KS6 with thin sections from Wadi Bani Awf, Oman Mountains.

160 160

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman

Since first possible Roadian fossil samples originate from around 15 m above the base of the KS6 the lowermost part could be even older. Thus it can be assumed that the time interval assigned to the ca. 190 m-thick KS6 approximately lies in between 4 and 7 million years (Myr) (Gradstein et al., 2012).

Cycles

With a thickness of ca. 1–5 m these cycles represent the smallest scale of cyclicity within the studied section. The most significant characteristics of these small cycles are regular vertical changes in texture, grain-size and fossil content. Four basic cycle types defined by Koehrer et al. (2010) can be identified in the KS6 (Figure 27).

Cycles consist of a transgressive and regressive hemicycle. The facies stacking pattern varies along the depositional gradient. Due to the predominant foreshoal environment in KS6 they tend to have high grain content at cycle boundaries and highest mud content around the zone of maximum flooding. A maximum of 61 cycles are recorded in Sequence KS6, and given the 4 to 7 Myr duration, they would represent ca. 65–115 Kyr (thousand year) in duration. As in other Khuff sequences (Koehrer et al., 2012; Walz and Aigner, 2012) these 1–5 m cycles possibly represent a fifth-order or 100 Kyr Milankovitch signal.

Foreshoal Cycle Description: This cycle type is variably composed of stacked muddy, normally graded mid- to outer- ramp facies types. The thinner lower part usually consists of bioturbated mudstones to wackestones with various ichnofabrics. The thicker upper part consists of graded, commonly bioturbated packstones to mudstones and bioclastic packstones showing erosive bases and scour surfaces. These may turn upwards into massive, low-angle laminated peloidal packstones to grainstones.

Interpretation: Dark burrowed mudstones at the base of the cycle type indicate fully open-marine/ deeper intra-shelf offshoal conditions and maximum relative water depth in a low-energy depositional environment around SWB. The rise to fall turnaround (interval of maximum accommodation) is defined at intervals with a maximum percentage of mudstone textures. The facies stacking pattern of the regressive hemicycle suggests a transition from the outer ramp to distal to proximal foreshoal. This typical shallowing-upward trend is associated with an increase of energy indicators, grain-size and sorting.

Shoal Margin Cycle Description: The lower part usually consists of graded storm beds or bioturbated mudstones to wackestones. The upper part mainly starts with amalgamated, graded bioclastic packstones to grainstones with erosive bases and frequent scouring. These may be overlain by thicker, low- angle laminated peloidal packstones to grainstones or coarse-grained intra-clastic grainstones and rudstones. They show an increase in sorting compared to the lower, bioturbated bioclastic packstone.

Interpretation: This cycle type represents the transition from a storm-dominated foreshoal environment to a higher-energy, shallower shoal flank setting. Common bed amalgamation and scouring reflects lower accommodation, accumulations of peloids and low-angle lamination suggest high-energy conditions and sediment input from the adjacent shoal complex. Storm influence is interpreted from normally graded beds. The regressive maximum occurs at the top of the packstone to grainstone that represents the time of maximum depositional energy and minimum accommodation.

Shoal Cycle Description: The thin lower part of the cycle type is represented by sheets of muddy mid- to outer- ramp facies such as graded/bioturbated storm beds. The thicker upper part starts with thick beds of bioclastic packstones to grainstones, thin layers of scoured graded beds or graded low-angle laminated peloidal packstones to grainstones. Upward these sediments may pass into massive amalgamated intra-clastic grainstones/rudstones. In most cases, facies grade into well sorted and cross-bedded peloidal or oolitic grainstones. In some cases the grainstones are overlain by microbial laminites or thin muddy lagoonal deposits.

161 161

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al.

as rather SWB FWWB FWWB possibly reworking. (FWWB) by indicates between base In 0.5 cm mudstones sheet s wave KS6. azoic mudstones Offshoal bioclasts Zoophycus muds and massi ve lo wer weather bioclastic the other in in fair the facies various fossils produced of 0.5 cm above and proximal lack events occurs corals The most storm the with facies burrows. these, energy represent together crinoidal gardens Coral blankets and Form 0.5 cm Zoophycus highest deposits thrived. The with Distal foreshoal grainstones Rugose Coral these have KS6. partly TRANSGRESSIVE SYSTEM crinoidal developed blankets of Sequence 0.5 cm Bivalve Shell mudstones not of bioclastic sheets coral in part Graded storm beds and Proximal foreshoal and abundance resulted apparently Peloid is The gardens transgressive setting conditions Model. lower 0.5 cm Gastropod crinoidal the Facies in backshoal (SWB) a low-energy Crinoid base Since SWB conditions ransgressive T wave below shoals. Shoal 25: Foraminifer 0.5 cm storm Crinoidal shoal fshoal, Figure open-marine crinoidal and Of indicates anoxic conditions.

162 162

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman and SWB FWWB foreshoal corals of Offshoal Most Mudstone 0.5 cm absence (partly bioturbated) the , dev eloped. are biodiver sity backshoal low by 0.5 cm muddy a indicated and as shoal prevail to 0.5 cm seem peloidal/oolitic storm sheets Rugose Coral Foreshoal Graded peloidal/bioclastic conditions REGRESSIVE SYSTEM cross-bedded A restricted 0.5 cm Bivalve Shell facies. system Peloid peloidal regressive and the In oolitic 0.5 cm Gastropod peloid shoal of X-bedded oolite/ model. Shoal facies Crinoid abundance the Regressive and Bioturbated 26: microbial mats Foraminifer 0.5 cm mudstone-wackestone/ Backshoal deposits show grading (storm deposits) and abundant bioturbation. crinoids, Figure

163 163

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al.

Interpretation: The shallowing-upward trend, associated with an increase of energy, sorting and the change from skeletal to peloidal or oolitic grains, is interpreted as a prograding shoal body. The regressive maximum mainly occurs at the top of the peloidal-oolitic grainstone that represents the time of maximum depositional energy. Cycle caps are partly composed of muddy beds that represent further shoaling into a lower accommodation lagoonal/intertidal setting.

Shoal-to-Backshoal Cycle Description: Cycles of this type are dominated by grainy textures, i.e. stacked physically stratified shoal flank to shoal complex beds. Layers of densely packed, often imbricated rip-up clast packstones pass rapidly into graded packstones that turn into high-angle trough cross-bedded peloidal/oolitic grainstones during an overall coarsening up trend. These rather massive grainstone bodies make up the gross of this cycle type. They are covered, often above a sharp boundary, by microbial laminites. These “microbial-caps” may be replaced by bioturbated “muddy-caps” in some cases.

Interpretation: A coarsening up trend from shallow-water muddy beds to thick shoal beds in the lower part of the cycle records a clear increase in accommodation. It is associated with a rapid landward stepping of shoal bodies. The arrival of laminites or bioturbated wackestones heralds the seaward stepping of low-energy, shallow-water backshoal facies deposits during falling relative sea level.

Cycle Sets

In many outcrop sections, 6–15 m-thick cycle sets are the most obvious and easiest scale of cycles to recognize. They are often well-reflected in the gamma-ray signal, which makes them the smallest possible units to correlate (cf. 2-D correlation section). Cycle sets display the lateral movement of facies associations or belts triggered by medium-term changes of relative sea level according to Walther’s Law. While the lower, transgressive part of most sections shows offshoal (LFA 8) to foreshoal (LFA 7) transitions within medium cycles, the upper regressive part of Khuff Sequence KS6 shows the prevalence of repeated and retrogradation of shoal complexes (LFA 5). Sequence KS6 contains a maximum of 17 cycle sets. Dividing the time interval of 4 to 7 Myr by this number would result in 235–412 Kyr per cycle set. Thus cycle sets possibly represent the fourth-order or 400 Kyr Milankovitch signal.

Sequence

The entire KS6 shows an overall, 190 m-thick transgressive-regressive sequence. The lower part starts with terrestrial facies (Basal Saiq Clastics, Figures 16 to 18) and transitions into foreshoal carbonates (Bioclastic Facies, Figures 19 and 20). Around the maximum flooding, offshoal facies predominates (deeper-water facies, Figures 21 and 22) whereas the regressive part is composed of shoal to backshoal facies (oolitic/peloidal grainstones, Figure 23 and 24). High-frequency sequences as detected in Sequence KS5 (Walz and Aigner, 2012) and KS4 (Koehrer et al., 2012) could not be identified unequivocally.

Bounding Surfaces

Base Sequence Boundary: The base of Sequence KS6 is the Sub-Saiq Unconformity, where Permian sediments (Basal Saiq Clastics and Saiq carbonates) overlie Proterozoic strata. The unconformity below the Basal Saiq Clastics is erosive where high-energy sediments such as sandstones (Wadi Sahtan) or conglomerates (wadis Bani Awf and Bani Hajir) predominate. On the Saiq Plateau the contact is not erosive but the siltstones that are present there show intense rooting. In the context of the increasingly open-marine overlying sediments, the Basal Saiq Clastics represent the most proximal, most shallow facies in KS6 and their base is a clear sequence boundary.

Top Sequence Boundary: Microbial laminites (Saiq, Wadi Sahtan and Hajir) and rooted mudstones form the top of Sequence KS6.

164 164

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman SWB FWWB Seaward Sea Level

Peloidal packstone to grainstone Graded packstone to wackestone (tempestite) Bioturbated mudstone

A LF Offshoal PG extur e T

MW

s Cycle Foreshoal cycle Foreshoal cycle

Intraclastic rudstone Peloidal packstone to grainstone Graded packstone to wackestone (tempestite)

A LF PG Foreshoal erent cycle types, shoal- to backshoal, shoal, extur e Shoal margin cycle T

MW Cycles Shoal margin cycle CYCLE TYPES 2 m 0 OUTCROP ell-sorted

Bioturbated mudstone Graded packstone to wackestone (tempestite) W oolitic grainstone

A LF Shoal cycle PG extur e T Shoal

Shoal cycle MW Cycles shoal margin and foreshoal cycle, have been defined by Koehrer et al., 2010 . ell sorted

Burrowed mud- to wackestone W oolitic grainstone Microbial laminites Bioclastic/peloidal packstone to grainstone

A LF Figure 27: Outcrop cycle types along a carbonate ramp. The 4 di ff Shoal- to backshoal cycle PG xtur e Backshoal Te

MW Cycles Landward Shoal- to backshoal cycle

165 165

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al.

Cycle Set Boundaries: In the context of the shoal to foreshoal environment of Sequence KS6, cycle set boundaries are often erosive surfaces. Higher-energy facies, such as grainstones or proximal storm deposits, erode lower energy facies, like distal storm and mud rich deposits.

2-D CORRELATION

Five sections in the Oman Mountains were correlated on a scale of 30 x 50 sq km. A mainly W-E transect runs from Wadi Sahtan to Wadi Mistal, and a N-S transect from the Saiq Plateau to Wadi Mistal. In general the Khuff epeiric carbonate ramp dips towards the north (Ziegler, 2001) (Figure 1). The deposition of the lower part of KS6, however, was probably strongly controlled by the initial topography created by the Sub-Saiq Unconformity.

Various correlation scenarios were applied based on sequence-stratigraphic interpretations as well as on facies and gamma-ray trends. The effect of the different correlation scenarios on potential reservoir geometries are discussed, modeled and displayed in the 3-D modeling section. Most correlations use the top KS6 (“Microbial Marker 1”) as datum (Figures 28 to 30b), some use the maximum flooding surface (Muddy Marker) (Figures 31 and 32).

The thickness of KS6 varies considerably from west to east (Figures 28 to 32) with the lowermost part of the KS6 seems to be missing in the east. These thickness changes can be explained by a paleohigh in the east. Consequently correlation lines in all correlations onlap onto the interpreted paleohigh in the east (Figures 28 to 32).

Correlation based only on Gamma-ray Data Although GR values cannot be directly translated into Dunham texture, the values tend to increase around muddy facies and decrease around grainy facies. Mudstones form in low-energy conditions and are therefore more likely to contain clay, which creates higher GR values in the sedimentary record. Figure 28 illustrates how the cyclic pattern of the GR signal allows a simple peak-to-peak correlation. The correlation based on GR only results in a layer-cake pattern with the maximum flooding surface being interpreted at a GR high with good correlatability. Although the GR signal seems cyclic throughout most of Sequence KS6 a pure GR correlation probably reveals only lithostratigraphic time lines and does not follow paleolandscape surfaces.

Correlation based on Dunham Texture and Components This correlation scenario focused on the Dunham carbonate textures, key components and marker beds (Figure 29). The main marker bed is formed by an interval that can be found in all logged sections: the so-called “Muddy Marker”. Although its appearance varies with respect to the abundance of mudstones, it is traceable through all sections; its thickness increases towards the west. A second marker is the microbial bed at the top of Sequence KS6 (“Microbial Marker 1”) in Wadi Sahtan and Wadi Bani Awf, which laterally changes into a rooted mudstone in Wadi Hajir and Wadi Mistal.

The occurrence of specific fossils was further used to correlate the KS6 sections. Examples are crinoidal floatstones (rather distal setting), which are correlated with crinoidal grainstones (more proximal setting). In the lower KS6 facies, texture and fossils apparently follow a proximal-to-distal trend from east to west. The upper portion of Sequence KS6 generally shows less lateral heterogeneity than the interval beneath the maximum flooding.

Compared to a pure GR correlation, correlation lines follow proximal-to-distal trends from east to west. However, the resulting correlation shows that some of the timelines dip in different directions. Such a pattern would be relatively unlikely since the time lines do not consistently follow paleo- landscape surfaces.

Integrative Correlation In this scenario, correlation lines were chosen to honor both GR and textural trends (Figures 30a and 30b). In addition, time lines strictly mimic proximal-to-distal paleolandscape surfaces, and

166 166

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman

follow sequence-stratigraphic principles. During the transgressive hemi-sequence, an aggradational- retrogradational stacking pattern with onlaps onto the paleohigh is apparent. In contrast, the regressive hemi-sequence shows a progradational stratigraphic architecture with shingles above the maximum flooding, with correlation lines dipping towards Wadi Sahtan in the west (dipping angle < 0.06°). Additional fifth-order cycles appear below the top of the KS6 towards the more distal sections in Wadi Sahtan to the west and the Saiq Plateau to the south forming toplaps. A fifth-order cycle pinch-out (downlap) around the maximum flooding surface can be observed in Wadi Sahtan. Toplap, downlap features and the pinch-out of cycles are indicated by correlation and can not be observed on an outcrops scale.

Influence of Sub-Saiq Unconformity In the context of a syn-rift phase, cycles pinch-outs and thickness variations are probably influenced by a combination of differential subsidence and initial topography. The key factor for the correlation of KS6 appears to be the paleorelief associated with the Sub-Saiq Unconformity. The first deposition of Khuff carbonates occurred in paleolows forming onlapping geometries on paleohighs. This depositional pattern explains the strong thickness variations in KS6 with thicknesses of up to 190 m in paleolows and about 90 m around potential paleohighs. Initial topography also contributed to facies variations between different sections (e.g. rooted mudstones in Wadi Mistal and Wadi Hajir turn laterally into microbial laminites in Wadi Bani Awf and Wadi Sahtan).

The undulating initial paleotopography probably affected different depositional trends and the shapes of facies belts. More grainstones appear around the paleohigh in the east and more muddy facies types in Wadi Sahtan to the west. The maximum dipping angle of correlation lines of 0.06° appears steep compared to KS4 to KS1 (Koehrer et al., 2012) where the maximum dip angle is around 0.001°. The steeper-dipping correlation time lines in KS6 are probably a result of the underlying topography.

Correlation of Potential Reservoir Bodies Previous studies in Khuff sequences KS1 to KS3 on the Saiq Plateau in Oman (Zeller et al., 2011) revealed that single-grainstone bodies may show pinching and swelling geometries on a 1 x 2 sq km scale, but reflect an overall layer-cake geometry. Correlations on an 8 x 8 sq km scale for sequences KS1 to KS3 on the Saiq Plateau (Koehrer et al., 2011) indicates single shoal bodies (LFA 5) to be traceable if they reach a certain thickness. The KS6 correlation in the Oman Mountains is on a subregional scale (30 x 50 sq km). Sparse data coverage (five sections with an average spacing of ca. 15 km) creates high uncertainty for the correlatability of single grainstone bodies. The upper, regressive part of KS6 contains more than 50% of grainstones throughout all sections. Based on this observation it can be assumed that the percentage of grainstones in between the sections is about as high. However this assumption does not reveal the lateral geometry of single grainstone beds.

In order to visualize the potential differences between laterally extensive and laterally limited grainstone bodies, two correlation scenarios for grainstone units are displayed in Figures 31 and 32. These scenarios visualize two possibilities on how laterally extensive or confined grainstone bodies could be. They mirror how the interpretation on this scale by the correlating geologist can influence reserve estimations and production forecasts. Observations in the field could only prove that individual grainstones beds are traceable over hundreds of meters. Both illustrations do not incorporate the basal part of KS6, and therefore the number of fourth-order cycle sets is different to the number in other figures. In Figure 31 laterally continuous grainstones (ca. 15 km correlation lengths) form large, connected reservoir bodies. In the lower part, grainstones mainly composed of bioclasts, do not correlate over wide areas; they are confined by the paleohigh to the east and pinch-out to the west. In the upper part of KS6 where peloidal/oolitic grainstones predominate in all sections, larger volumes and a high connectivity can be assumed.

In the correlation shown in Figure 32, single potential reservoir bodies were assumed to form laterally more confined patches with horizontal correlation lengths of about 5 km. Even in upper KS6, where there is an overall very high percentage of grainstones, the single units seem often laterally disconnected in 2-D. Although single-grainstone units might not be very extensive, the overall chance

167 167

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al.

PURE GAMMA-RAY CORRELATION

West 16 km Wadi Bani 8 km Wadi Bani 14 km East Wadi Sahtan Awf Hajir Wadi Mistal 190 Top KS6 180

170 160 150 140

130 MFS 120 110 100 90

Depth (meter) 80 70 60 50 Proterozoic Basement 40 Gamma-ray 0.06° API 30 10 40 m 20 Dip angle 10 km 0 Figure 28: Cycle set correlation using gamma-ray data only (Top KS6 was used as datum).

West PURE DUNHAM CORRELATION East Wadi Wadi Bani Wadi Bani Wadi Sahtan 16 km Awf 8 km Hajir 14 km Mistal 190 Top KS6 180

170 160 150 140 MFS 130 120 110 100 90 Depth (meter ) 80 70 60

50 Proterozoic 40 Basement 0.06° 30

20 Dip angle m 10 km 0 Figure 29: Facies-based cycle set correlation (Top KS6 was used as datum).

168 168

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman

(a) West INTEGRATIVE CYCLE SET CORRELATION East

Wadi Sahtan Wadi Bani Wadi Bani Wadi Mistal 16 km Awf 8 km Hajir 14 km 190 Top KS6 180

170 160 150 140 130 MFS 120 110 100 90 Depth (meter) 80 70 60 Proterozoic Basement 50 40 0.06°

30 m Dip angle 20 10 km 0

(b) North South Texture Wadi Mistal 22 km Saiq Gamma- 190 Ray (API) ckestone Top KS6 -order

180 th Cycles 4 Rudstone Floatstone Boundstone Grainstone Packstone Wa Mudstone 10 40 170 160 150 MFS 140 130 120 110 100 90

Depth (meter) 80 70 60 Figure 30: (a) Integrative cycle set correlation comprising GR, texture and 50 fossil content in east-west direction 40 0.06° (Basal Saiq Clastics are displayed in 30 orange, top KS6 was used as datum). m 20 Dip angle (b) Integrative north-south oriented cycle 10 set correlation (Basal Saiq Clastics in km Proterozoic Basement 0 orange, top KS6 was used as datum).

169 169

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al. East 40 0.04 ° 0.03 ° 0.02 ° 0.01 ° 0.00 °

di Wa Mistal Ra y (API )

Gamma-

10 Cycles

-order -order 4

th Mudstone

8 km ckestone Wa

Dip angle

Packstone

Grainstone

Reservoir Facies Boundstone extur e

14 km T A Correlation: Large Extent of Floatstone

LF Rudstone 5 m l di Bani Hajir Shoa l Fore - Offshoal Backshoa Wa Proterozoic Basement 8 km f Aw di Ban i Wa 16 km MFS op KS6 T est-east cross-section. Reservoir facies association correlates over 15 km. This figure highlights the potential ge n eral reservoir architecture est-east cross-section. Reservoir facies association correlates over 15 km. This figure highlights the di Sahtan Proterozoic Basement Wa st We Figure 31: W 0 0

90 80 70 60 50 40 30 20 10 and does not include the two additional cycles in the lowermost clastic part of Sequence KS6 (the maximum flooding surface, "MFS " was used as datum). and does not include the two additional cycles in lowermost clastic part of Sequence KS6 (the maximum flooding

11 100 150 140 130 120 Depth (meter) Depth

170 170

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman r was East 40 0.04 ° 0.03 ° 0.02 ° 0.01 ° 0.00 ° reservo i di

"MFS" Wa Mistal Ra y (API)

Gamma-

10 Cycles

general

-order -order 4

th

surface, Mudstone

8 km ckestone Wa

Dip angle Packstone

potenti al Grainstone

Reservoir Facies floodin g Boundstone

the extur e

14 km T A Correlation: Short Extent of Floatstone

LF Rudstone 5 m maximum highlights (the l figure KS6, di Ban i Hajir Shoa l Fore - Offshoal Backshoa This Wa km. Sequence 5 of Proterozoic Basement ca. part of 8 km clastic lengths have to lowermost f Aw di Ban i the Wa in assumed are cycles bodies additional two 16 km Reservoir MFS the op KS6 T section. include not cross does Proterozoic Basement est-east and W di Sahtan 32: Wa st We used as datum). Figure architecture 0 0

90 80 70 60 50 40 30 20 10

11 15 0 14 0 13 0 12 0 10 0 Depth (meter) Depth

171 171

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al.

of amalgamation is expected to be very high when the overall percentage of grainstones is sufficiently high (upper KS6 > 50%). This implies for the upper KS6 that even if grainstones are restricted in their lateral extent, they may still form connected potential reservoirs.

3-D FACIES MODELING

To illustrate possible reservoir facies distribution and depositional trends, facies associations (backshoal, shoal, foreshoal and offshoal) from outcrops logs were imported into Petrel to create a 3-D model.

Framework: Grid, Zones and Layers Since outcrop observations have shown that the reservoir facies seem laterally continuous over kilometers, the horizontal cell size was set to 1 km x 1 km. After the correlation of reservoir bodies in 2-D (Figures 30a and 30b), tops of cycle sets were correlated in Petrel. The resulting surfaces where transformed into horizons and zones were assigned in between. Therefore most zones represent one complete fifth-order cycle. Only around the maximum flooding interval (“Muddy Marker”) down- lapping hemicycles are interpreted, which led to the described shingle-like architecture (Figure 30a, in between Wadi Sahtan and Wadi Bani Awf).

Subsequently zones were subdivided into layers. Depending on thickness, a high number of layers were allocated to each zone to ensure good vertical resolution. The resulting models comprise 16 zones, 1,380 layers and a total number of about 2 million cells. Each cell has an average size of 1.0 by 1.0 sq km and a thickness of 0.13 m. The smallest recorded bed in the field measures about 0.1 m in thickness.

Facies Modeling Each cell in the model was populated with a value for each facies association by upscaling the well logs using the “most-of averaging” method, so the predominant facies fills the whole cell. Vertical facies proportions derived from the outcrop data were used for all models to automatically define the percentage of each facies in each zone. Thus the proportion of each facies association within one zone (equivalent to one cycle set) is derived from the proportions in outcrops.

The method “Truncated Gaussian with Trends” (TGS) modeling algorithm was used to areally distribute facies associations within the model. The TGS method was selected because it combines statistical methods with geological concepts such as facies trend maps. Specific trend maps were generated for each zone. Since a variogram cannot be calculated with the presented data density, an isotropic variogram was chosen (5 and 15 km range). As the correlation revealed (e.g. Figures 30a and 30b) facies belts in KS6 follow a proximal-to-distal trend from the paleohigh in the northeast (Wadi Hajir and Wadi Mistal) to south and west (Saiq Plateau, Wadi Bani Awf and Wadi Sahtan).

Modeled Horizontal Facies Ranges and their Impact on Connectivity Previous studies suggest that thicker grainstone bodies are laterally continuous for more than 8 km (Koehrer et al., 2011). Other studies, e.g. the Triassic Muschelkalk in SW-Germany, (Palermo et al., 2010), showed that grainstone bodies could be traceable over distances of up to 10–20 km. To take the different possibilities into account, facies modeling was performed with different lateral facies variograms. Two simulations were run and compared with respect to the resulting connectivity of potential reservoir facies (Figures 33 and 34).

1. Lateral extent of grainstone bodies 15 km: isotropic variogram with 15 km horizontal facies range 2. Lateral extent of grainstone bodies 5 km: isotropic variogram with 5 km horizontal facies range

Both models portray the general vertical deepening-shallowing trend in the KS6 (Figures 33 and 34): Some shoal grainstone facies is present in the lower KS6, barely any in middle KS6 around the MFS while the upper part of KS6 is dominated by shoal facies association. Most offshoal facies can be observed near the maximum flooding interval. Some backshoal facies can be found at the top of

172 172

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman

the KS6. Laterally the most proximal facies concentrates around Wadi Mistal with a deepening trend towards southern and western sections (Saiq Plateau, Wadi Sahtan). Shingles can be observed around the maximum flooding and at the top of the KS6 (Figure 33 and 34, marked in yellow).

Volumes and Connectivity A potential reservoir volume was delineated, calculating the percentage of potential reservoir facies association in each zone. Approximately 40% of KS6 consists of potential reservoir facies (i.e. shoal LFA, grainstones). The regressive part of KS6 consists of more than 50% potential reservoir facies, the lowermost part of 33%, and the middle transgressive around the maximum flooding of only 6% (Table 2). Although the overall percentage of potential reservoir facies stays the same in all models because the vertical proportion is honored, the lateral extent of grainstones can make a difference when it comes to the connectivity between single grainstone shoal bodies.

Wadi Bani Wadi Bani Wadi Mistal Hajir Awf Figure 33: Model with a Wadi Sahtan Saiq Plateau facies range of 15 km. Yellow -0 -0 lines mark shingles. KS6 -40 facies model shows realistic -40 facies successions and a -80 clear geological trend from -80 -120 proximal facies around the -120 paleohigh (Wadi Mistal) to -160 distal facies in paleolows -160 -200 (Wadi Sahtan and Saiq -200 530000 Plateau). 540000

550000 2560000 Facies 560000 Backshoal Foreshoal 2570000 N Shoal Offshoal 0 10 570000 Clastics 2580000 km

Wadi Bani Wadi Bani Wadi Mistal Hajir Awf Saiq Plateau Wadi Sahtan -0 -0

-40 Figure 34: Model with a -40 facies range of 5 km. Facies -80 -80 belts show a more patchy -120 appearance. Yellow lines -120 mark shingles. -160 -160 -200

200 530000 540000

550000 2560000 Facies 560000 Backshoal Foreshoal 2570000 N Shoal Offshoal 570000 0 10 Clastics 2580000 km

173 173

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al.

Table 2: Overall potential reservoir facies volumes and percentages.

Reservoir Volumes

Overall model volume of KS6 180.8 km3

Overall reservoir volume (grainstones) 34.8 km3 19.2%

Regressive part (zone 52) 26.4 km3 52%

Around maximum flooding surface (zone 6-14) 5.0 km3 6%

Lower part (zone 14-16) 4.8 km3 33%

Figure 35 illustrates how the lateral size affects the reservoir connectivity in single zones with different grainstone percentage. The models with 5 and 15 km facies range are filtered so that they show exclusively the shoal facies association (red). Two zones (A and B, see Figure 35) are presented in top view (A) from the top of the KS6, where grainstones have an abundance of over 50%; and (B) from the interval of maximum flooding where grainstones account for less than 10% of the facies. It can be observed that in zone A, where grainstones are abundant (> 50%), the lateral extent of single shoal bodies apparently does not have a significant effect on connectivity since grainstones amalgamate and form one, laterally connected unit. In zone B, where grainstones are rare they form single isolated volumes which are not inter-connected. Thus the lateral extent of single grainstone units does not effect the connectivity if the overall percentage of grainstones is high enough (approx > 50%).

15 KM FACIES RANGE 5 KM FACIES RANGE

AA A Saiq

Hajir Hajir Mistal Awf Mistal Awf N N

BBSaiq

Hajir Hajir

Mistal Awf Mistal Awf N N Figure 35: Zones from the top of the KS6 (A) and from the maximum flooding interval (B) are presented in top view for 5 km (right side) and 15 km (left side) lateral facies range. Potential reservoir facies is displayed in red. A: Zone from top KS6 with overall high percentage of reservoir facies (> 50%). Even smaller sized potential reservoir bodies (5 kms) seem to be connected (amalgamation). B: Overall low percentage of reservoir facies (< 10%). Lateral range of 15 km results in few larger bodies, lateral range of 5 km results in multiple, smaller and disconnected bodies.

174 174

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman

CONCLUSIONS

1) Five outcrop sections of Lower Khuff (KS6) time-equivalent strata (Saiq Formation) were investigated in the Oman Mountains (Al Jabal al-Akhdar, Sultanate of Oman). The KS6 can be subdivided into four different facies associations (backshoal, shoal, foreshoal, offshoal). The KS6 forms a third-order transgressive-regressive sequence, built by smaller-scale cycle sets and cycles.

2) The initial paleorelief apparently controlled the thickness and distribution of the Basal Saiq Clastics, and that of the overlying KS6 carbonates as well as their composition. The correlation strategy to follow paleolandscape surfaces using all available data shows that the stratigraphic architecture is aggradational with onlaps against the paleohigh during the transgressive hemi-sequence and displays subtle shingles during the regressive hemi-sequence.

3) This study revealed potential reservoir units in the KS6, commonly regarded as non-reservoir in the subsurface of Oman. In the transgressive part, the predominant reservoir facies are bioclastic crinoidal grainstones (with only poor diagenetic potential), concentrated around the margin of a paleohigh. In contrast, oolitic/peloidal grainstones in the upper regressive part (with a much higher diagenetic potential) have a much more widespread distribution.

4) A range of 3-D models were generated using different correlation scenarios and varying lateral extents of potential reservoir grainstones. The lateral continuity of shoal grainstones bodies seems to have no to little effect on their inter-connectivity if the overall percentage of reservoir facies is high enough (> 50%) as single reservoir units amalgamate.

ACKNOWLEDGEMENTS

This study is part of an extra-mural research project of the University of Tuebingen with Qatar Shell and Petroleum Development Oman. We would like to thank Shell and PDO for their financial support and M. Poeppelreiter, J. Amthor, A. Brandenburg, J.-M. Dawans, G. Forbes and J. Schreurs for their assistance. PDO and the Omani Ministry of Oil and Gas are thanked for permission to publish the paper. We are grateful to our Sedgeo members of the University of Tuebingen: L. Walz (now Shell), M. Haase (now ExxonMobil) and M. Bartenbach (now Statoil), in addition to P. Jeisecke who prepared the thin sections. H. Forke is thanked for biostratigraphic analysis of the thin-sections. Shuram Oil and Gas (Muscat) is acknowledged for fieldwork logistics. We are also very grateful to ALT and Schlumberger for providing access to WellCAD and Petrel software packages. GeoArabia’s Assistant Editor Kathy Breining is thanked for proofreading the manuscript, and GeoArabia’s Production Co- manager, Nestor “Nino” Buhay IV, for designing the paper for press.

REFERENCES

Al-Jallal, I.A. 1995. The Khuff Formation: Its regional reservoir potential in Saudi Arabia and other Gulf countries; depositional and stratigraphic approach. In M.I. Al-Husseini (Ed.), Middle East Petroleum Geosciences, GEO ‘94. Gulf PetroLink, Bahrain, v. 1, p. 103-119. Baud, A. and M. Bernecker 2010. The Permian-Triassic transition in the Oman Mountains. Gutech Geoscience Workshop Publication 1, IGCP 572 Field Guide Book 2, 109 p. Baud, A. and S. Richoz 2013. Stratigraphic Note: The Permian–Triassic transition and the Saiq/Mahil Boundary in the Oman Mountains: Proposed correction for lithostratigraphic nomenclature. GeoArabia, v. 18, no. 3, p. 87-98. Baud, A., B. Beauchamp and S. Richoz 2012. Paleokarst in the Saiq Formation (Saiq Plateau, Oman). Abstract book, International Conference on the of the Arabian Plate and the Oman Mountains 7-9th January, 2012. Sultan Qaboos University, Muscat, p. 51-52. Béchennec, F., J. Le Métour, J.-P. Platel and J. Roger 1993. Explanatory notes to the Geological map of the Sultanate of Oman, Scale 1:1,000,000. Directorate General of Minerals, Oman, Ministry of Petroleum and Minerals, 93 p.

175 175

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al.

Blendinger, W., A. Van Vliet and M.W. Hughes-Clarke 1990. Updoming, rifting and continental margin development during the late Palaeozoic in northern Oman. In A.H.F. Robertson, M.P. Searle and A.C. Ries (Eds.), The Geology and Tectonics of the Oman Region. Geological Society Special Publication no. 49, p. 27-37. Chauvet, F., T. Dumont and C. Basile 2009. Structures and timing of Permian rifting in the central Oman Mountains (Saih Hatat). Tectonophysics, v. 475, p. 563-574. Dunham, R.J. 1962. Classification of Carbonate Rocks According to Depositional Texture. In W.E. Hamm (Ed.), Classification of Carbonate Rocks, A Symposium. American Association of Petroleum Geologists, p. 108-121. Forbes G.A., H.S.M. Jansen and J. Schreurs 2010. Lexicon of Oman subsurface : Reference guide to the stratigraphy of Oman’s hydrocarbon basins. GeoArabia Special Publication 5, Gulf PetroLink, Bahrain, 371 p. Forke, H., M. Pöppelreiter, T. Aigner, B. Koehrer, L. Walz, D. Bendias and M. Haase 2012. Integrated biostratigraphy of the Saiq Formation (Al Jabal al-Akhdar, Oman Mountains) and its implication for the regional correlation of Khuff time-equivalent deposits. In The Permo–Triassic Sequence of the Arabian Plate, Abstracts of the EAGE’s Third Arabian Plate Geology Workshop, Kuwait. Abstract, GeoArabia, v. 17, no. 1, p. 230-234. Glennie, K.W., M.G. Boeuf, M.H.W. Hughes-Clarke, M. Moody-Stuart, W.F. Pilaar and B.M. Reinhardt 1973. Late Cretaceous nappes in Oman Mountains and their geologic evolution. American Association of Petroleum Geologists Bulletin, v. 57, p. 5-27. Glennie, K.W., M.G. Boeuf, M.H.W. Hughes-Clarke, M. Moody-Stuart, W.F. Pilaar and B.M. Reinhardt 1974. Geology of the Oman Mountains. Verhandelingen Koninklijke Nederland geologisch mijnbouwkundig Genootschap, v. 31, 423 p. Gradstein, F.M., J.G. Ogg, M. Schmitz and G. Ogg 2012. (Eds.) The Geological Time Scale 2012. Elsevier, 1,144 p. Henderson, C.M. and S.L. Mei 2003. Stratigraphic versus environmental significance of Permian serrated conodonts around the Cisuralian-Guadalupian boundary: New evidence from Oman: Palaeogeography Palaeoclimatology Palaeoecology, v. 191, p. 301-328. Kerans, C., and S.W. Tinker. 1997. and characterization of carbonate reservoirs. Society of Economic Paleontologists and Mineralogists Short Course Notes, no. 40, 130 p. Koehrer, B., M. Zeller, T. Aigner, M. Poeppelreiter, P. Milroy, H. Forke and S. Al-Kindi 2010. Facies and stratigraphic framework of a Khuff outcrop equivalent: Saiq and Mahil formations, Al Jabal al-Akhdar, Sultanate of Oman. GeoArabia, v. 15, no. 2, p. 91-156. Koehrer, B., T. Aigner and M. Pöppelreiter 2011. Field-scale geometries of Upper Khuff reservoir geobodies in an outcrop analogue (Oman Mountains, Sultanate of Oman). Petroleum Geoscience, v. 17, no. 1, p. 3-16. Koehrer, B. T. Aigner, H. Forke and M. Pöppelreiter 2012. Middle to Upper Khuff (Sequences KS1 to KS4) outcrop-equivalents in the Oman Mountains: Grainstone architecture on a subregional scale. GeoArabia v. 17, no. 4, p. 59-104. Konert, G., A.M. Al-Afifi, S.A. Al-Hajri, K. de Groot, A.A. Al Naim and H.J. Droste 2001a. Paleozoic stratigraphy and hydrocarbon habitat of the Arabian Plate. In M.W. Downey, J.C. Threet and W.A. Morgan (Eds.), Petroleum Provinces of the Twenty-first Century, The Pratt Conference II, Association of Petroleum Geologists, Memoir 74, p. 483-515. Konert, G., A.M. Al-Afifi, S.A. Al-Hajri and H.J. Droste 2001b. Paleozoic stratigraphy and hydrocarbon habitat of the Arabian Plate. GeoArabia, v. 6, no. 3, p. 407-442. Kozur, H.W. and B.R. Wardlaw 2010. The Guadalupian conodont fauna of Rustaq and Wadi Wasit, Oman and a West Texas connection. Micropaleontology, v. 56, no. 1, p. 213-233. Leven, E. Ya. 1997. Permian stratigraphy and Fusulinida of Afghanistan with their paleogeographic and paleotectonic implications. Geological Society of America, Special Paper 316, p. 1-133. Leven, E. Ya. 2009. The Upper Carboniferous (Pennsylvanian) and Permian of the Western Tethys: fusulinids, stratigraphy, biogeography. Transactions of the Geological Institute, v. 590, p. 1-238. Maurer, F., R. Martini, R. Rettori, H. Hillgärtner and S. Cirilli 2009. The geology of Khuff outcrop analogues in the Musandam Peninsula, United Arab Emirates and Oman. GeoArabia, v. 14, no. 3, p. 125-158. Montenat, C., A.F. de Lapparent, M. Lys, H. Termier, G. Termier, D. Vachard 1976. La transgression permienne et son substratum dans le Jabal Akhdar (Montagnes d’Oman,Péninsule Arabique). Annales de la Société Géologique du Nord, Lille, France, v. 96, no. 3, p. 239-258.

176 176

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Mid-Permian Khuff Sequence KS6, Oman

Obermaier, M., T. Aigner and H.C. Forke 2012. Facies, sequence stratigraphy and reservoir/seal potential of a Jilh Formation outcrop equivalent (Wadi Sahtan, Triassic, Upper Mahil Member, Sultanate of Oman). GeoArabia, v. 17, no. 3, p. 85-128. Osterloff, P., A. Al-Harthy, R. Penney, P. Spaak, G. Williams, F. Al-Zadjali, N. Jones, R. Knox, M.H. Stephenson, G. Oliver and M.I. Al-Husseini 2004. Depositional sequence of the Gharif and Khuff formations, subsurface Interior Oman. In M.I. Al-Husseini (Editor), Carboniferous, Permian and Early Triassic Arabian Stratigraphy. GeoArabia Special Publication 3, Gulf PetroLink, Bahrain, p. 83-147. Palermo, D., T. Aigner, S. Nardon and W. Blendinger 2010. Three-dimensional facies modeling of carbonate sand bodies: Outcrop analog study in an epicontinental basin (Triassic, southwest Germany). American Association of Petroleum Geologists Bulletin, v. 94, no. 4, p. 475-512. Pillevuit, A. 1993. Les blocs exotiques du Sultanat d’Oman. Evolution paléogéographique d’une marge passive flexurale. Mémoires de Géologie, Lausanne, v. 17, p. 1-249. Pöppelreiter, M.C., C.J. Schneider, M. Obermaier, H.C. Forke, B. Koehrer and T. Aigner 2011. Seal turns into reservoir: Sudair equivalents in outcrops, Al Jabal al-Akhdar, Sultanate of Oman. GeoArabia, v. 16, no. 1, p. 69-108. Rabu, D., F. Béchennec, M. Beurrier and G. Hutin 1986. Explanatory notes to the geological map of the Nakhl Quadrangle, Sultanate of Oman. Geoscience map, scale 1:100,000, sheet NF 40-3E. Ministry of Petroleum and Minerals, Directorate General of Minerals, Sultanate of Oman. 42 p. Rabu, D., P. Nehlig, J. Roger, F. Béchennec, M. Beurrier, J. Le Métour, C. Bourdillon de Grissac, M. Tegyey, J.-J. Chauvel, C. Cavelier, H. Al-Hazri, T. Juteau, D. Janjou, B. Lemiere, M. Villey and R. Wyns 1993. Stratigraphy and structure of the Oman Mountains. Document du Bureau de Recherches Géologiques et Minières, no. 221, 262 p. Searle, M.P. 2007. Structural geometry, style and timing of deformation in the Hawasina Window, Al Jabal al Akhdar and Saih Hatat culminations, Oman Mountains. GeoArabia, v. 12, no. 2, p. 99-130. Sharland, P.R., D.M. Casey, R.B. Davies, M.D. Simmons and O.E. Sutcliffe 2004. Arabian Plate Sequence Stratigraphy. GeoArabia, v. 9, no. 1, p. 199-214. Stampfli, G.M. and G.D. Borel 2002. A plate tectonic model for Paleozoic and Mesozoic constrained by dynamic plate boundaries and restored synthetic ocean isochrones. Earth and Planetary Science Letters, v. 196, p. 17-33. Strohmenger, C.J., A. Ghani, O. Al-Jeelani, A. Al-Mansoori, T. Al-Dayyani, L.J. Weber, K. Al-Mehsin, L. Vaughan, S.A. Khan and J.C. Mitchell 2006. High-resolution sequence stratigraphy and reservoir characterization of Upper Thamama (Lower Cretaceous) Reservoirs of a giant Abu Dhabi oil field, United Arab Emirates. In P.M. Harris and L.J. Weber (Eds.), Giant hydrocarbon reservoirs of the world: From rocks to reservoir characterization and modeling. American Association of Petroleum Geologists, Memoir 88/SEPM Special Publication, p. 139-171. Walz, L. and T. Aigner 2012. Khuff Sequence 5 (KS5), Oman Mountains: Lateral facies and sequence variability – a record of differential subsidence? In The Permo–Triassic Sequence of the Arabian Plate, Abstracts of the EAGE’s Third Arabian Plate Geology Workshop, Kuwait. Abstract, GeoArabia, v. 17, no. 1, p. 250-252. Walz, L., T. Aigner and B. Koehrer 2013. Khuff sequence KS5 outcrop equivalents in the Oman Mountains, Sultanate of Oman: Variations to the simple “layer-cake“ stratigraphy. GeoArabia, v. 18, no. 4, p. 179-218. Weidlich, O. and M. Bernecker 2011. Biotic carbonate precipitation inhibited during the Early Triassic at the rim of the Arabian Platform (Oman). Palaeogeography, Palaeoclimatology, Palaeoecology, v. 308, p. 129-150. Weidlich, O. and M. Bernecker 2012. Earthquake-triggered post-depositional deformation at the rim of the Arabian Platform (Permian–Triassic, Oman Mountains). In The Permo–Triassic Sequence of the Arabian Plate, Abstracts of the EAGE’s Third Arabian Plate Geology Workshop, Kuwait. Abstract, GeoArabia, v. 17, no. 1, p. 253-256. Zeller, M., B. Koehrer, E.W. Adams, M. Pöppelreiter and T. Aigner 2011. Near well-scale heterogeneities in a Khuff outcrop equivalent (Saiq Plateau, Jebel Al Akhdar, Sultanate of Oman). Journal of Petroleum Geology, v. 34, no. 3, p. 241-260. Ziegler, M.A. 2001. Late Permian to Holocene paleofacies evolution of the Arabian Plate and its hydrocarbon occurrences. GeoArabia, v. 6, no. 3, p. 445-504.

177 177

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021 Bendias et al.

ABOUT THE AUTHORS

Daniel Bendias studied Geosciences at the University of Tuebingen (Germany). His diploma thesis (2010) at the University of Tuebingen was on the Paleorelief-influenced facies and sequence patterns in the Lower Khuff time equivalent strata (Sultanate of Oman). He is currently working as Research Associate and PhD student at the Center for Applied Geosciences (University of Tuebingen). His PhD thesis, funded by Petroleum Development Oman (PDO), focuses on sequence stratigraphy, reservoir and seal geobodies of the Jurassic Mafraq Formation, a mixed carbonate siliciclastic system in outcrops and subsurface of Oman. [email protected]

Bastian Koehrer is a Development Geologist in Wintershall’s German business unit, working on mature oil field and tight gas sands development in the German North Sea and Lower Saxony. He has more than five years of E&P project experience in Germany, Oman, Qatar and the UAE with a professional track record in both carbonate and clastic reservoirs. Bastian obtained a PhD degree (2011) in Carbonate from the University of Tübingen (Germany) in research collaboration with Shell (Qatar) and Petroleum Development Oman. For his PhD dissertation on the Khuff Formation, Bastian spent 18 months of outcrop mapping in the Sultanate of Oman. Bastian is a member of the EAGE, AAPG and DGMK and has published several papers on carbonate sequence stratigraphy and reservoir outcrop analogs. [email protected]

Michael Obermaier studied Geosciences at the Universities of Tuebingen (Germany), and Miami (Florida, USA). His PhD thesis (2013), a research cooperation between the University of Tuebingen and Petroleum Development Oman was on Triassic reservoir and seal characterization in outcrops and subsurface of Oman. Since 2013 Michael has been working as a Carbonate Geologist for Shell Global Solutions in Rijswijk, The Netherlands. [email protected]

Thomas Aigner studied Geology and Paleontology at the Universities of Stuttgart, Tuebingen/Germany and Reading/UK. For his PhD dissertation on storm depositional systems (1985) he worked at the Senckenberg- Institute of Marine Geology in Wilhelmshaven (Germany) and spent one year at the University of Miami in Florida (USA). He then became an Exploration Geologist at Shell Research in Rijswijk/Holland and Houston/ Texas focussing on basin analysis and modelling (1985–1990). Since 1991 Tom has been a Professor and Head of the Sedimentary Geology Group at the University of Tuebingen. In 1996 he was a “European Distinguished Lecturer” for the AAPG. His current projects focus is on sequence stratigraphy and reservoir characterisation/modelling in outcrop and subsurface. [email protected]

Manuscript received October 17, 2012 Revised March 20, 2013 Accepted March 27, 2013

178 178

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/18/3/135/4567666/bendias.pdf by guest on 25 September 2021