The effects of diagenesis on the petrophysical and geochemical attributes of the Asmari Formation, Marun oil field, southwest Iran

The distribution of good reservoir quality and its causes is the main challenges in carbonate reservoir characterization. This study investigates the effects of diagenetic processes on the reservoir quality of the carbonate successions of the Asmari Formation, in the Marun oil field, southwest Iran. The study applies an integrated approach including core petrography, petrophysical rock typing, stable carbon and oxygen isotopes as well as major and trace elements analyses. Petrographic studies and geochemical analyses express that the Asmari limestones have been affected mainly by compaction, dissolution, recrystallization, calcite and anhydrite cementation and dolomitization. Among those diagenetic overprints, dolomitization and dissolution played an important role to enhance the reservoir quality of the formation. Moreover, four types of dolomites were recognized and the rate of dolomitization increases toward the top of the Asmari carbonate successions. Possible models for dolomitization include mixing zone, brine reflux, seepage reflux and tidal pumping of seawater. Employing Flow Zone Index and Discrete Rock Type concepts led to classification of the Asmari reservoir into seven reservoir rock types. Integrating reservoir rock typing with petrographic studies and geochemical analyses also confirms that reservoir quality of the Asmari Fm. would have been mainly controlled by diagenetic processes. Moreover, stable isotopes, trace elements and facies analyses support the idea that carbonate intervals of the Asmari Formation were deposited in a warm, shallow-water environment under a saline condition.

Carbonate intervals of this formation are important hydrocarbon reservoirs, which have been producing oil and gas more than a century and reaching to the latter stages of their producing lifecycle. Therefore, an increasing interest in enhanced oil recovery methods to optimize the production has been raised recently. One of the initial stages in this process is obtaining an improved understanding of the carbonate reservoir characteristics (Van Buchem et al. 2010).
The aim of carbonate reservoir characterization is to define porosity types, reservoir heterogeneities and flow units for the purposes of reservoir simulation. Rock typing is a major step to evaluate carbonate reservoir characterization. However, this is a complex approach in carbonate reservoirs due to the effects of a variety of factors such as diagenesis, depositional texture and fracturing on the heterogeneity and fluid flow of these types of reservoirs (Askari and Behrouz 2011;Xu et al. 2012;Bize-Forest et al. 2014).
Reservoir rock type quality depends on the either type of porosity or pore texture and the modality size of the pore throat radius for carbonate rocks (Bize-Forest et al. 2014). In the Marun oil field, texture and porosity of the carbonate successions have been significantly changed due to digenetic overprints (compaction, cementation, dissolution, dolomitization) and tectonic fracturing (Aqrawi et al. 2006;Ranjbaran et al. 2007). Although digenesis essentially affects the petrophysical attributes of the carbonate reservoirs, the importance of diagenetic processes with respect to relative sea level changes, and their effects on the flow units in the Marun oil field, has not been studied in detail.
Using routine core analyses, petrographic study of cores and thin sections as well as geochemical properties, the present paper investigates the effects of diagenetic processes, (especially, compaction, dissolution, cementation, dolomitization), and depositional texture on the petrophysical attributes of the Asmari Fm. carbonate intervals (porosity types and permeability), in order to define reservoir quality distribution and rock types.

Structural setting
The Marun anticline is located in the central part of the of the Zagros thrust belt (NE sector of the Arabian plate) (Fig. 1a). Motion of the Arabian plate toward the northeast started in the Lower Miocene, about 25 Ma (Homke et al. 2010) or 22 Ma (ArRajehi et al. 2010) and resulted in gentle folding of the Mesozoic to Eocene sedimentary strata in the Zagros basin (Fig. 1a, b). The first episode of folding took place during deposition of the Gachsaran Fm. evaporates, probably as a result of the initiation of continent-continent collision in the Early Miocene (Sherkati et al. 2005). Following this first phase, a period of tectonic quiescence during the Middle to Late Miocene accompanied by the deposition of the Mishan Fm. marls and lower Aghajari Fm. sandstones. The main episode of folding resulted in growth strata in upper Aghajari sandstones during the Miocene-Pliocene. The last tectonic event took place during the Pliocene-Pleistocene due to activation of reverse basement faults and resulted in the construction of the Zagros Mountains (Sherkati and Letouzey 2004; Leturmy et al. 2010). The Zagros fold-and-thrust belt (ZFTB) comprises three geographic provinces: Lurestan, Fars and the Dezful Embayment (Sepehr and Cosgrove 2004) (Fig. 1a). The ZFTB can be divided into two distinct domains from the SW to NE, with major changes in structural style and topography: (a) the Simply Folded Belt (SFB), which shows fairly regular wavelength folds that extend for hundreds of kilometers in the vicinity of the Persian Gulf (Falcon 1974;Sepehr and Cosgrove 2004;Mouthereau et al. 2006); the High Zagros, with higher elevation folding, a sharp increase in elevation and kilometer-scale throws on major thrusts. The Marun oil field is a supergiant (6 km × 70 km) fault-related, asymmetric anticline, located in the central part of the Dezful Embayment province, in the SFB of the Zagros Mountains (Fig. 1a,  c). The general trend of the anticline is NW-SE, parallel with the Zagros folded belt. Structural dip reaches 45° and 65° in the north and south flanks, respectively. It is bounded by a high angle reverse fault in the south flank.

Stratigraphic framework
Deposition of the Asmari Fm. occurred on the margins of a preexisting Eocene platform which surrounded a deep-water basin in which marls, shales and deep-water limestones of the Pabdeh Fm. were deposited (James and Wynd 1965) (Fig. 1a, b). In other words, Oligocene moderate to shallow-water carbonate-dominated sediments that were deposited around the edge of the preexisting Eocene platform are marked as the Asmari Formation, whereas deep-water basinal facies are assigned to the Pabdeh Formation in the center of basin (Ehrenberg et al. 2007). By progressive infilling of this basin, the Asmari platform prograded over the Pabdeh Formation so that the Asmari/Pabdeh boundary is diachronous, (Ehrenberg et al. 2007) and becoming younger basin-ward. Sedimentation of fine clastics continued in the troughs and shallow marine carbonates and evaporites to the northeast margin of the NW-SE trending Zagros basin, while coarse clastic sediments (Ahwaz Sandstone Member) were shed from the southwest (Arabian plate outcrops) into the basin during the Oligocene (Motiei 1993). The Asmari Fm. was dated Oligocene-Lower Miocene (34-18 Ma) using strontium isotopes (Ehrenberg et al. 2007). The formation comprises 314 m of cream-to brown-colored limestones, with thin-bedded sandstones and shales in the lower section, which overlies the Pabdeh Fm. basinal mudstone at the type locality (Motiei 1993). In the Marun oil field, the Asmari reservoir comprises 450 m of carbonate rocks (shallow marine, lagoonal, intertidal to sabkha mudstone to grainstones) in the middle and upper sections, and shale with sandstone and carbonate layers in the lower section (Seyrafian et al. 2011;Avarjani et al. 2015). The Asmari Fm. Carbonate sedimentation continued mainly during the Lower Miocene, followed by deposition of the lower Fars Group (Gachsaran Fm.) evaporites. Mishan and Aghajari Fm. (Middle and Upper Fars Groups) clastic sediments were deposited through the Middle Miocene to Lower Pliocene. The Pasadanian orogenic phase caused folding and thrusting 1 3 toward the southwest of the entire sedimentary cover. Coarse clastic sediments, the Bakhtyari conglomerate, filled the syncline and other low lying areas during the Upper Pliocene to Pleistocene (James and Wynd 1965). The Bakhtyari conglomerate and Quaternary alluviums were transported from highs to the inter-mountain areas and foreland of the Zagros Orogenic Belt due to the uplift of the Zagros basin (Sherkati and Letouzey 2004) (Fig. 1a, Sepehr and Cosgrove (2004)]. b Cenozoic stratigraphy columns of the Zagros basin. The most representative stratigraphy column for the study area has been highlighted in light red [from Beydoun et al. (1992)]. c location of the studied wells in Marun oil field S7 S10 S11 S12 S13 S14 S15 S16 S17 S18 S19

Materials and methods
Petrographic studies were completed on 1800 thin sections from 7 wells (well numbers 8, 342, 124, 330,131, 338 and 30) (Fig. 1c). Thin sections have been stained by using alizarin red. Core samples from three wells (450 m in Mn#8, Mn#342 and 124) were described in terms of lithology, texture and sedimentary structures, visible diagenetic features (stylolites, calcite veins, solution seams, vugs and anhydrite nodules) and ichnofossils (Fig. 2). Carbonate lithofacies were described based on Dunham (1962), Embry and Klovan (1971) and Archie (1952) classifications. Twelve samples of dolostone and limestone were examined using a LEO 1455VP scanning electronic microscope in order to discriminate texture and porosity types and sizes. Different generations of cements and dolomites were identified using a CITL-8200 MK4 cathodoluminiscence microscope, with a gun current of 400 µA and 14 kV potential difference. After detailed petrographic descriptions (Table 1), forty samples (26 samples of limestone and dolomite from well 8, six samples from well 342 and eight samples of dolomite from wells 124, 330 and 30) were selected for geochemical analyses from cores of the Asmari Fm. These samples were analyzed for major and trace element compositions (Mg, Ca, Sr, Na, Mn and Fe). Stable isotope (δ 13 C and δ 18 O) analyses were performed for 35 samples by G.G Hatch Stable Isotope Laboratory, University of Ottawa, Canada.
Samples were weighed into exetainers, 0.1 ml of H 3 PO 4 was added to the side, and exetainers were capped and helium-flushed while horizontal. Reaction at 25 °C for 24 h (calcite) or 50 °C for 24 h (dolomite) was followed by extraction in continuous flow mode. The measurements were performed on a Thermo Finnigan Delta XP and Gas Bench II. Analytical precision (2σ) is ± 0.1 per mil. Data for carbon and oxygen were normalized using international standards NBS-18, NBS-19 and LSVEC (C only). Porosity and permeability data of core plugs were obtained from NIOC data base in order to investigate reservoir quality and rock type. Winland R35 and FZI methods were used for reservoir rock type classifications.

Petrographic evidences of diagenetic features
Petrographic descriptions of thin sections show that the Asmari Fm. is composed mainly of limestones (mudstones to grainstones with calcite cement), and very fine to coarse crystalline dolomites, with intervals of sandy dolomites and sandstones with calcite, dolomite and anhydrite cements (Figs. 2,3,4). Stratigraphy and brief petrographic descriptions of cores and selected samples are provided in Table 1 and Figs. 1b and 2.
The main diagenetic processes that have affected the Asmari limestones are compaction, stylolitization, dissolution, calcite cementation and dolomitization. Additional processes including micritization of bioclasts, anhydrite cementation, bioturbation, pyrite formation and neomorphism are common. Vertical distributions of diagenetic processes in the Asmari reservoir are presented in Fig. 4.
The main diagenetic processes that affected the Asmari reservoir quality are described below.

Compaction
Stylolitization as firm evidence for compaction can be seen throughout the Asmari Fm. Macro-and micro-scale stylolites occurred mainly parallel to the bedding plane in the upper portions of the Asmari Fm., with some oblique stylolites present in lower sections (Fig. 2). Clay seams, silt size quartz grains and organic materials are visible along stylolite surfaces (Fig. 3a1, a2, d1). The stylolites can be closed or semi-open, creating permeability barriers for hydrocarbon migration within the dense, argillaceous limestone intervals of the Asmari reservoir. Evidence for compaction is also provided by interlocked grain contacts seen in grainstones and sandstones (Fig. 3a3). The amount of porosity reduces in limestones and dolomites with increasing burial depth and overburden pressure generally (Fig. 5).

Cementation
Calcite cementation comprises both an earlier phase of finely crystalline, isopachous cement that lines pore walls and grain surfaces (Fig. 3b1), and a later phase of coarser, blocky spar that fills former pore volumes not completely filled by the fine cement linings that are common in limestone intervals (Fig. 3b1, b2, b6). Blocky calcite cement also occurs, filling former moldic pores in some dolostone samples (Fig. 3b2, b6). Concave meniscus, isopachous equant spary calcite and syntaxial overgrowth cements are very common in sandstones (Fig. 3b4, b6, b8). Fine to medium, crystalline dolomite cement has filled the intergranular porosity of sandstones as well as moldic and fracture porosities (Fig. 3b3, b6). In sandstones, two phases of siliceous cement are visible: (a) pore filling; and (b) quartz overgrowth, both of which have resulted in porosity reduction ( Fig. 3a3, b8). Anhydrite cement occurs mainly as filling fractures and porosities in carbonates, and intergranular porosity of sandstones especially in the upper sections of the Asmari Fm (Fig. 3b5). Anhydrite can also occur as nodules and patches (Fig. 3d2).

Dissolution and fracturing
Skeletal and non-skeletal grains within the carbonate matrix have been widely affected by dissolution processes, which have resulted in increasing porosity. The evidence for this can be seen as solution vuggy porosity in mudstones and dolostones ( Fig. 3c1, c2, c5). Dissolution has affected essentially all bioclasts of presumed original aragonite composition: coral fragments, mollusk shells and green algae. The moldic porosity formed by dissolution varies greatly in abundance from bed to bed, due to variations in bioclasts proportions and calcite cementation. Extensive dissolution of all types of bioclasts (including non-aragonite types) is especially common in dolostones, producing a distinctive micro-vuggy fabric (Fig. 3c5 (Fig. 3c4).

Dolomitization
Dolomitization has widely affected the limestones in the middle and upper sections of the Asmari Fm. The rate of dolomitization increases from the lower to upper units of this formation; however, the crystal size decreases in an upwards direction (Figs. 2, 6). Dolomite crystal sizes range from very fine to very coarse, have a unimodal to polymodal distribution and are euhedral planar to planar-s and xenomorphic (e.g., Sibley and Gregg 1987

Trace elements and stable isotope data
The results of trace elements and stable isotopes (δ 13 C and δ 18 O) are shown in Tables 1 and 2. The δ 18 O and δ 13 C values of the Asmari Fm. dolomites range from − 2.4 to 2.1‰ VPDB and − 0.2 to 4.3‰ VPDB, respectively. The values of δ 13 C and δ 18 O in limestones (mudstone to grainstone) range from − 6.8 to − 0.2‰ VPDB and − 3.3 to 2.5‰ VPDB, respectively. The Mg concentration ranges between 5.5 and 11.2% by weight in dolomites and increases upward (Fig. 8), but is less than 3.5% in dolomitic limestones.
A general positive trend can be seen between Mn and Fe molar concentrations ( Fig. 9), with Mn and Fe values generally increasing with depth; however, there is no relationship between dolomite type and molar concentrations of Mn and Fe ( Fig. 8 and
In the present paper, hydraulic units are calculated based on the FZI method (Kozeny 1927;Amaefule et al. 1993), and the complementary DRT method using the following equations (Abedini 2011). where RQI is the rock quality index (μm), k is permeability (md), ϕ is effective porosity (v/v), ϕz is normalized porosity, and FZI is the Flow Zone Index (μm). The results of calculations are presented in Table 3.
The log-log plot of RQI versus normalized porosity shows six regression lines and one point, which correspond to seven reservoir rock types (Fig. 11a, b). The FZI, RQI, DRT values and extracted reservoir rock types are shown in Table 3.
Rock types 1, 2 and 3 are mainly comprised of dolomite types of I, II and III. Rock types 4, 5 and 6 encompass a variety of limestone textures (mudstone, wackestone, packstone, grainstone, boundstone and some dolomite samples (Table 3).

Implications for paleoenvironments
The middle section of the Asmari Fm. is composed mainly of packstones, grainstones and boundstones, with mudstone  Table 3).
The microfacies which have been recognized in the Asmari Fm. suggest that deposition occurred in the spectrum of different depositional environments from coastal and peritidal settings to shallow marine environment. (e.g., Flügel 2004;Vaziri-Moghaddam et al. 2006). The depositional texture of the Asmari Fm. was affected by several digenetic processes during marine to burial and subaerial conditions (Figs. 2, 3, 4).
Stable isotope values vary between − 6.8 and 2.0‰ PDB for δ 18 O and − 3.3 to 4.3‰ PDB for δ 13 C (Table 2). Samples were divided into two groups (A and B) based on stable isotope values in order to recognize depositional and diagenetic paleoenvironments (Fig. 12). The characteristics of these two paleoenvironments are described below.

Depositional environment
The first group (group A) of samples comprises the least altered limestones and dolomites, with near zero or positive values of δ 18 O and δ 13 C. These samples suggest warm shallow marine and evaporitic conditions for the Asmari carbonates deposition (Fig. 12) (Adabi and Rao 1991;Swart 2015).
Modern tropical aragonitic sediments have low Mn (< 20 ppm), moderate Na (2500 ppm) and high Sr (10,000 ppm) concentrations. In contrast, modern temperate calcite sediments contain low Sr (~ 3000 ppm), high Na (~ 5000 ppm) and Mn (~ 150 ppm) concentrations (Rao 1986). In the present study, the average concentrations for trace elements (Sr, Na and Mn) are ~ 471 ppm, ~ 1175 ppm and ~ 124 ppm, respectively (Table 2, Fig. 9). The Sr/Na ratio is about 0.4, while it is ~ 0.7 for modern temperate calcite sediments and high (~ 4.0) in aragonitic carbonates. Thus, the geochemical properties of the Asmari carbonate sediments in Marun oil field differ from recent tropical carbonates, and there is no exact correlation with modern temperate carbonate sediments (Fig. 9). The Sr-Mn data support the results of stable isotope analysis which suggested that the Asmari Fm. was deposited in warm, shallow marine conditions (Fig. 9).
The general increase in Fe and Mn concentrations at drilling depths more than 3450 m might be due to redox conditions and pore fluid sources (Al-Aasm 2000) (Fig. 8). Widespread high Fe values (> 1000 ppm) and/or high Mn (> 50 ppm) could indicate that the pore waters were reducing in a media with significant source of Fe and/or Mn, such as soil horizons, underlying volcanic or detrital siliciclastics was present (Budd 1997;Azmy et al. 2001). Although some authors believe that high concentrations of Fe and Mn may occur due to recrystallization in dolomites (cf. Banner et al. 1988), however, results from the present study show that Fe and Mn values of analyzed samples have not been affect by the type and crystal size of dolomites.

Diagenetic processes sequences
The Zagros folding episode started in the early Miocene, but the main orogenic phase is presumed to take place in the Plio-Pliestocene (Sherkati and Letouzey 2004; Leturmy   (Ala 1982;Beydoun et al. 1992); hence, it can be concluded that much of the porosity evolved because of diagenetic processes before oil emplacement (Fig. 13). A discussion about the sequence of digenetic processes now follows.
Carbonates that form along continental margins and on carbonate platforms can eventually be buried to deeper depth. Assuming a normal geothermal gradient of 30 °C/ km, at burial depths of 2-3.5 km, the formation temperature increases to 60-115 °C. Hence, dissolution and precipitation reactions are enhanced (Swart 2015). Well data show that the Asmari Fm. is buried to depths more than 2 km (Fig. 2), so compaction as an early diagenetic process was an important factor that affected the Asmari Fm. from the early stages of deposition until the main orogeny phase during the Pliocene (Fig. 13). This process resulted in the porosity reduction with increase in drilling depth (Fig. 5). The abundance of stylolitization (pressure solution) in limestones and interlocked contacts of quartz grains in sandstones are the main evidences of compaction (Figs. 3a, 4). Stylolitization (pressure solution) has occurred, mainly parallel to the bedding with some oblique to vertical stylolitization in the upper and lower sections of the formation (Fig. 2). Dissolution of the host carbonates and loss of porosity occur during compaction (Rittenhouse 1971), and insoluble materials (including clay minerals, organic materials, pyrite, quartz grains) are left behind, leaving a distinct planar or ragged surface (Swart 2015). Stylolites usually act as permeability barriers; thus, they are of particular interest in reservoir characterization studies (Heap et al. 2014). However, stylolites are mainly oil stained, suggesting that stylolitization may have contributed to permeability enhance in the Asmari reservoir (Figs. 3a1, a2).
The most depleted values of δ 18 O and δ 13 C in limestones (Table 2 at Table 3 Table 2) (Budd 1997). The δ 18 O values that are lighter than the postulated marine calcite value (− 2.6‰ PDB) are due to equilibration with non-marine waters (Adabi and Rao 1991), Hence, the depleted δ 18 O values of the above samples could be due to subaerial exposure and the effect of meteoric water. Based on the core study, some intervals are located below the Type I or Type II sequence boundaries (samples 6,7,11,14 and 21;Figs. 2,3d). The depletion of δ 13 C is characteristic of vadose diagenetic media, where the soil CO 2 enrichment in 12 C occurs because of decomposition of organic matter (Adabi and Rao 1991). The depleted δ 18 O and δ 13 C values could also be attributed to temperature increase and preservation of organic materials (Azmy et al. 2001) (Fig. 8). The δ 13 C values in phreatic, meteoric calcite approach marine levels due to increasing rock-water interaction (Meyers and Lohmann 1985;Lohmann 1988). Trace element (e.g., Mn and Fe) concentrations show a co-variant relationship with δ 13 C values, but a weak correlation with δ 18 O values (Fig. 8). This correlation could be due to effect of vadose diagenetic zone. Some samples of more positive δ 13 C (up to + 4.28‰ PDB) and δ 18 O values close to marine equilibrium (Fig. 12) may reflect marine diagenetic reactions (Braithwaite and Camoin 2011).
Cementation as a very important process has occurred in the Asmari Fm. in different environments from subaerial, freshwater vadose, sabkha saline, marine and shallow to deep burial. Calcite and anhydrite cementation (Figs. 3b1, b2, b5) are more frequent in limestones, but dolomite and silica cements (Figs. 3b3, b4) can also be seen in sandstone intervals. Syntaxial overgrowth, concave meniscus and isopachous equant spary calcite cements (Fig. 3b4, b6, b8, c2) suggest vadose and freshwater phreatic zones (Halley and Harris 1979). Solution of quartz grains and precipitation of overgrowth silica cements could be evidence of shallow to burial diagenesis of sandstones (Fig. 3b8) (Walderhaug et al. 2009;Swart 2015). The blocky cement may be comprised of partly early, near-surface cement (related to exposure) and partly late (burial) cement, which may have formed as a result of dissolution along stylolites and fractures (Fig. 3b6,  c4). This process could have started at the beginning of the  (Adabi and Rao 1991) 1 3 Asmari Fm. deposition and lasted until oil migrated into the reservoir (Fig. 13). The Asmari Fm. has been affected by four phases of dolomitization (Fig. 6d, CL). The rate of dolomitization increases from lower to the upper parts of the Asmari carbonate succession. This is supported by increasing Mg/Ca ratios with a reduction in depth (Fig. 9). The dolomites of the Asmari Fm. are categorized as evaporitic because of their high concentrations of Na (2500 ppm > Na > 700 ppm) and sulfate (10,000 ppm > S > 900 ppm) (Staudt et al. 1993). The first phase (Phase I) of dolomitization is an early stage, which is documented by very fine, crystalline, unimodal dolomite (Type I with d < 20 μm) (Fig. 6a) and expresses pervasive replacement dolomitization that occurred in supratidal to sabkha and restricted supersaline environments (Budd 1997;Warren 2000).This phase of dolomitization occurred in muddy facies, such as mudstones, wackestones and mud-dominated packstone dolomites in grainy facies (cf. Al-Aasm 2000; Aqrawi et al. 2006). Some dolostone samples (e.g., No. 1, 4 and 11, Table 2) from the upper section of the Asmari Fm. show characteristics of this type of dolomite, with red luminescence, positive values of δ 13 C (+ 2.4‰, + 2.5‰ and + 4.1‰), and near zero values of δ 18 O PBD (− 0.2‰, − 0.3‰ and + 0.1‰) (Tables 1, 2, Fig. 6a). This phase mainly took place during lower Miocene (Burdigalian or later).
Pervasive dolomitization phases (II and III) were formed in a burial environment and high temperature conditions  Table 2). These two phases of dolomitization may have occurred during burial of the Asmari Fm. at depths of more than 2 km after deposition of the Aghajari Fm. in late Miocene. Phase IV is the precipitation of medium to coarse dolomite crystals in fractures, and solution vuggy porosities. This phase comprises planar polymodal dolomite with positive values (+ 1.7 to + 4.3) of δ 13 C and δ 18 O (0.0 to + 2.3), and orange to brown luminescence (Fig. 6d, Table 2). The enriched values of δ 18 O may show the evaporative origin of Mg-bearing fluids (Al-Aasm et al. 2009). Phase IV of dolomitization presumably occurred after the main Zagros folding and fracturing events, during late Miocene to Pliocene.
It is suggested that the reflux model can explain the fine crystalline, pervasive dolomites found mainly in the upper half successions of the Asmari Fm. This model also is supported by the presence of anhydrite nodules with dolomite (Figs. 2, 3a1, d2), and the texture of the muddy facies includes tidal flat (supratidal and intertidal) and subtidal/ lagoonal facies of the precursor limestone. Another possible model for dolomite formation is probably related to flush over of platform grainstones from the previous depositional cycle (particularly within the Middle Asmari Fm.) by condensed and evaporative solutions. Dolomitization of open-marine facies in the lower and middle sections of the Asmari Fm. may have occurred during the later burial phase, when the compaction of the underlying Pabdeh marls and shales generated Mg-rich fluids. The presence of concentric overgrowth of dolomite crystals (crystal zoning) defined by cathodoluminescence can mark the replasive burial origin (Fig. 6d) (cf. Budd 1997).
In addition to invoking the reflux model to explain the fine crystalline dolomite (Type I), the mixing zone model is also suggested for the origin of dolomitization in the Asmari Fm. because of the covariance of δ 18 O‰ and δ 13 C‰ (Figs. 2, 8 and Table 2) (Budd 1997;Swart 2015). This fluid mixing can be the main mechanism for dolomitization near the subaerial exposure surfaces and sequence boundaries (Al-Aasm 2000). Petrographic studies show evidence of freshwater dissolution and development of dolomitization below the subaerial exposures (Figs. 2, 3c1, c2, d) (proposed sequence boundaries by Avarjani et al. 2015). During lowstand systems tract, percolation of freshwater into the underlaying carbonate beds and formation of a mixing zone with the seawater resulted in the pervasive dolomitisation (types II and III) of the Asmari Fm. Dolomite type IV may have precipitated because of percolation of the Mg-bearing fluids through fractures and vuggy porosities (Fig. 6d).
The variation of Sr concentrations (90-2195 ppm) in dolomites can be related to either depositional and digenetic conditions of mixed meteoric realms, or burial digenetic media with Sr-rich fluids (Table 2) (White and Al-Aasm 1997;Al-Aasm 2000) that could indicate dolomite formation in a water-buffered diagenetic system (Budd 1997  Permeability, porosity and bulk density are used from routine analysis of core plugs ( Fig. 9), it can be inferred that the original chemistry of the precursor carbonate is well preserved during diagenesis (Azmy et al. 2001).
Compaction and cementation are the main causes for reservoir quality destruction (Fig. 3a, a2, b5). Reduction in the Asmari reservoir porosity versus depth can be seen generally (Fig. 5). A large proportion of the primary intergranular porosities in grainstones and packstones are occluded by calcite and anhydrite cements (Fig. 3c2, b6), while anhydrite and dolomitic cements are the main cause of porosity reduction in sandstones (Fig. 3b3, b4, b5).
The petrophysical properties of selected limestone and dolomite core plug samples are presented in Table 3. This dataset shows that the average porosity and permeability, respectively, are 4.4% and 2.16 md for limestones, and 11.80% and 7.29 md for dolomites. Although a large increase in porosity (approximately 7.4%) is related to dolomitization, there is a fairly good correlation between porosity and permeability for dolomites (Fig. 10b). This may be due to the closing of some pore throats by over dolomitization or precipitation of anhydrite cement (Lucia 1983).
On contrary to dolomites, limestone samples represent a good correlation between porosity and permeability ( Fig. 10a) that reflects diagenesis is more positive for permeability preservation in limestones rather than dolomites, or it can be evidence of permeability preservation in limestones during diagenesis.
A cross-plot of porosity versus permeability, overlaid with standard curves of the relationship between petrophysical properties and particle size for uniformly cemented, nonvuggy carbonates (Lucia 1983), shows that Type II dolomites (20 < d < 100 µm) match with the standard area for this particle size. However, Types I and III do not match with the standard area of the corresponding particle size (Fig. 10c). This deviation to the left of Type 1 dolomites may be due to porosity and permeability enhancement by dissolution and development of vuggy porosity (Fig. 10c). Interparticle porosities in Type I dolomite samples are mainly connected to each other by touching vugs that led to high effective porosities (Figs. 3c1, c2). The porosity and permeability of Type III dolomites were reduced due to cementation and over dolomitization, as if the plots of these samples deviated to the right (Figs. 6c, 10c).
The cross-plot of porosity and permeability in limestones (Fig. 10d) reveals that the petrophysical properties of the Asmari reservoir differ from the original rock texture.
Reservoir rock typing, based on FZI and DRT methods, illustrates that the reservoir hydraulic flow units do not  (James and Choquette 1983), (Adabi and Rao 1991 Fig. 13 Sequence of diagenesis processes and porosity evolution relative to each other, folding of Zagros belt and oil emplacement necessarily follow the geological classification and texture of carbonate reservoir rocks (Table 3). Therefore, it seems that diagenesis is the main controlling factor of the carbonate reservoir quality. Reservoir rock type 3 (RRT3) consists of Type I and II dolomites as well as sandstones with dolomitic cement (Table 3), indicating the best reservoir quality (Fig. 11).

Conclusions
Seven hydraulic flow units were distinguished using FZI (Flow Zone Index) and DRT (Discrete Rock Type) methods. Rock types 1, 2 and 3 mainly consist of dolomite, while rock types 4, 5, 6 and 7 encompass a variety of limestone textures (mudstone, wackestone, packstone, grainstone, boundstone) and some dolomite samples. Cross-plots of petrophysical properties and analysis of FZI show that the main controlling factors of the reservoir heterogeneity and quality in the Asmari Fm. are diagenetic processes, especially, cementation, dolomitization and dissolution. Core plug data reveal that the porosity and permeability of the Asmari reservoir increased due to dolomitization (average increase of 7.4% and 5 md, respectively) as the most effective diagenetic factor influencing reservoir quality enhancement.
Four types of dolomites have been recognized in the Asmari Fm.: (1) Type I, very fine to fine crystalline dolomite (crystal size < 20 μm); (2) Type II, fine to medium crystalline (< 100 μm), polymodal, planar to planar-s secondary pervasive dolomite; (3) Type III, medium to coarse crystalline, polymodal, subhedral to xenotopic; and (4) Type IV, medium to very coarse crystalline, dolomitic cement or saddle dolomite. Type I dolomite would have been formed in restricted sabkha and supersaline environments, while the other types would have been formed in medium to deep burial and high temperature conditions. Based on dissolution porosity, the development of dolomitization under subaerial exposures and Type 1 sequence boundaries the mixing zone model is proposed for the main dolomitization phase in the Asmari Fm.
However, brine reflux, seepage reflux and tidal pumping of seawater models are also proposed to explain dolomitization in some other intervals of the Asmari Fm.
Stable isotope and trace elements analyses indicated that the carbonate successions of the Asmari Fm. were deposited in warm, shallow marine environments under a saline evaporative condition.