Integrated rock typing and pore facies analyses in a heterogeneous carbonate for saturation height modelling, a case study from Fahliyan Formation, the Persian Gulf

Knowledge of initial fluids saturation has great importance in hydrocarbon reservoir analysis and modelling. Distribution of initial water saturation (Swi) in 3D models dictates the original oil in place (STOIIP), which consequently influences reserve estimation and dynamic modelling. Calculation of initial water saturation in heterogeneous carbonate reservoirs always is a challenging task, because these reservoirs have complex depositional and diagenetic history with a complex pore network. This paper aims to model the initial water saturation in a pore facies framework, in a heterogeneous carbonate reservoir. Petrographic studies were accomplished to define depositional facies, diagenetic features and pore types. Accordingly, isolated pores are dominant in the upper parts, while the lower intervals contain more interconnected interparticle pore types. Generally, in the upper and middle parts of the reservoir, diagenetic alterations such as cementation and compaction decreased the primary reservoir potential. However, in the lower interval, which mainly includes high-energy shoal facies, high reservoir quality was formed by primary interparticle pores and secondary dissolution moulds and vugs. Using huge number of primary drainage mercury injection capillary pressure tests, we evaluate the ability of FZI, r35Winland, r35Pittman, FZI* and Lucia’s petrophysical classes in definition of rock types. Results show that recently introduced rock typing method is an efficient way to classify samples into petrophysical rock types with same pore characteristics. Moreover, as in this study MICP data were available from every one meter of reservoir interval, results show that using FZI* method much more representative sample can be selected for SCAL laboratory tests, in case of limitation in number of SCAL tests samples. Integration of petrographic analyses with routine (RCAL) and special (SCAL) core data resulted in recognition of four pore facies in the studied reservoir. Finally, in order to model initial water saturation, capillary pressure data were averaged in each pore facies which was defined by FZI* method and using a nonlinear curve fitting approach, fitting parameters (M and C) were extracted. Finally, relationship between fitting parameters and porosity in core samples was used to model initial water saturation in wells and between wells. As permeability prediction and reservoir rock typing are challenging tasks, findings of this study help to model initial water saturation using log-derived porosity.


Introduction
Calculation of initial water (S wi ) and hydrocarbon (1-S wi ) saturations has important implications in petrophysical analyses of petroleum reservoirs (Cannon 2015;Kennedy 2015;Tiab et al. 2015). Reservoir modellers use initial hydrocarbon saturation in order to model past and future behaviour of the reservoir and consequently to make an appropriate development plan (Ghedan et al. 2004;Esther et al. 2015). Generally, a combination of resistivity log data, porosity, water resistivity and volume of shale is required to calculate the water saturation in each well (Kennedy 2015;Tiab et al. 2015). However, limited log-derived water saturations are not suitable for modelling, especially in reservoir with long production history and thick transition zone.
The common and accepted method for initial saturation modelling is based on special core analysis (SCAL) and capillary pressure (PC) tests (Jamiolahmady et al. 2007;ZHAO et al. 2008;Kumar et al. 2011). The capillary pressure against gravity force controls the fluids distribution above the FWL in a reservoir. When capillary pressure data are available, water saturation can be calculated at any point above the FWL (Tiab et al. 2015). This process is called saturation height modelling (Al-Khaldi et al. 2012;Yong et al. 2012). However, in a real reservoir, water saturation does not show a single idealized vertical profile showing increase in hydrocarbon saturation in predictable pattern upward from free water level (FWL), which can be seen in text books. This is because of the heterogeneous nature of porous media (Lalicata et al. 2012;Moore et al. 2013;Ma 2019). Changes in lithofacies, storage and flow capacities and pore characteristics (i.e. pore types, size, shape and geometry) will affect the water saturation vertical profile and its three-dimensional distribution in subsurface area (Aliakbardoust et al. 2013). Figure 1 represents a schematic diagram showing different rock types with different saturations versus height profile in a heterogeneous reservoir and highlights using appropriate rock typing approach. As it can be seen from this figure, Point "A" located in the transition zone must be classified as rock type 1, which is prone to be classified in rock type 4 with implementing wrong rock typing method.
Carbonate reservoir rocks are heterogeneous by nature in various scales because of their formation processes and complex diagenesis history (Lucia 2007;Moore et al. 2013;Mehrabi et al. 2015;Tavoosi-Iraj et al. 2021). Consequently, the water saturation height functions are more complex in carbonate reservoir rocks than the clastic reservoirs (Lian et al. 2016;Ma 2019). Accordingly, the more reservoir heterogeneity results in higher possibility of inaccurate calculation of initial hydrocarbon saturation (Cuddy et al. 1993;Al Waili 2009;Chudi et al. 2010). Therefore, using an adequate rock typing method is vital for saturation modelling in a reservoir (Askari et al. 2011;Mirzaei-Paiaman et al. 2015). Net-pay cut-off determination, core sample selections for SCAL tests and permeability prediction in uncored intervals are among the other important applications of rock typing (Amaefule et al. 1993;Abbaszadeh et al. 1996;Mirzaei-Paiaman et al. 2015;Mirzaei-Paiaman et al. 2016). FZI, r35Winland, r35Pittman, FZI* and Lucia's petrophysical classes (Kolodzie Jr et al. 1980;Pittman 1992;Amaefule et al. 1993;Lucia et al. 2001;Mirzaei-Paiaman et al. 2018)  are indices that utilize porosity and permeability to separate reservoir into distinct rock type. However, their efficiency needs to be checked in every reservoir. As characterization of pore heterogeneity is critical to the prediction of flow behaviour under reservoir conditions, this study investigates pore geometry in Fahliyan Formation to model initial water saturation.
Over the years, various mathematical formulas are introduced and used to describe variations in Sw with height and reservoir properties. Leverett (1941) introduced J-Leverett function to normalize all the capillary pressure versus saturation curves into a single curve. Aufricht et al. (1957) were first to relate Pc and Sw data from core measurements to porosity and permeability information. Heseldin et al. (1974) modified the Aufricht et al. 1957 by relating the porosity to hydrocarbon bulk volume. Alger et al. (1989) modified the Heseldin's method by proposing a multilinear regression to relate Pc to porosity or permeability and used the term Cap-Log method. Gradually, new saturation functions were introduced which were written based on height and permeability Sondena 1992), height and porosity (Cuddy et al. 1993) or merely based on height (Skelt et al. 1995).
In recent years, researchers have focused on using different saturation functions based on various frameworks such as lithofacies and rock typing of reservoir interval based on different indices (Kamalyar et al. 2013;Omeke et al. 2014;Lian et al. 2016;Kundu et al. 2017). However, it should be noted that permeability prediction and rock typing are challenging issues in heterogeneous formations.
Recently, Lian et al. (2016) have used the idea of capillary height and Pittman r35 to map of water saturation. A similar approach has been adapted by Dakhelpour-Ghoveifel et al. (2019). They have used capillary pressure curves and specific borders between rock types to calculate water saturation and assign rock type to grids in transition zone. Here, we used pore facies and parameter fitting approach to model water saturation.
This paper aims to determine initial water saturation in a heterogeneous carbonate reservoir in northwest of Persian Gulf. First, a detailed petrographic analysis on 388 thin sections, taken from 290 m of Fahliyan Formation, resulted in the definition of six microfacies, dominant diagenetic events and types of pores in this reservoir. Then, in order to evaluate reservoir quality and pore geometry, reservoirs rocks were classified into groups based on the different indices and appropriate methods were selected. As pore system properties controls fluid saturation and the reservoir rock quality, we used pore types as a framework to extract fitting parameters and their relationship with porosity in core plugs. Then, saturation height modelling was conducted using the extracted equation, and its constants were modified based on porosity in core samples. In this way, substitution of core porosity with log-derived initial water saturation can be calculated anywhere in the reservoir.
The advantage of this method is that it uses rock types only in plug samples scale to extract fitting parameters. Therefore, calculated initial water saturation from this method can be used for reservoir simulation and reservoir rock typing in reservoir scale using methods introduced by other researchers (Dakhelpour-Ghoveifel et al. 2019). On the contrary, this method utilizes vast SCAL and RCAL core analysis results and it is not applicable in fields without core samples.
First, the reservoir stratigraphy is presented. Then, used data and methods are described. Thereafter, in "Results" section, microfacies, dominant diagenetic events and types of pores were recognized. Next, pore facies were determined through implementing different methods on our data set and finally discussion and conclusions are presented.

Reservoir stratigraphy
The studied field is located in northwest of the Persian Gulf (Fig. 2a). The main reservoir rock of this field is Berriasian-Valanginian strata known as the Fahliyan Formation, which is partly equivalent to the Yamama, Minagish and Manifa Formations in the Middle East (Fig. 2b). This formation unconformably overlies the Hith evaporites and is conformably overlain by the tight and micritic limestones of the Gadvan Formation (Fig. 2b). In the studied field, based on the reservoir quality, Fahliyan Formation has been divided into three members including lower Fahliyan (LF) with the best reservoir quality that is equivalent to the Manifa Formation in surrounding areas, the middle Fahliyan (MF) representing the major reservoir interval and the upper Fahliyan (UF) with poor reservoir quality (Fig. 4). According to previous studies, this formation and its equivalents were deposited on homoclinal carbonate ramps on the NE margin of the Arabian Plate (Adabi et al. 2010;Jamalian et al. 2011;Jamalian et al. 2014;Noori et al. 2019). It is mainly composed of grey to brown, massive, oolitic to pellety limestones with minor brecciation at the basal parts (James et al. 1965). The extensive development, moderate to good storage and flow capacities and proximity to the organic-rich Sargelu and Garau Formations all make this formation as one of the most important reservoir rocks in the Zagros area and the Persian Gulf region (Setudehnia 1978;Ghazban 2007

N N A A R R I I K KU UW WA AI IT T
Basinal marl/shale Siliciclastics A B

Materials and methods
This research is based on a dataset from the reservoir interval of the Fahliyan Formation in an oilfield selected from the northwest of Persian Gulf (Fig. 2a). In particular, the dataset of this study comprises 290 m of core samples and 388 thin sections for petrographic analyses. All samples were dyed by blue epoxy resin in order to differentiate pore types. In particular, all thin section samples were analysed by petrographic microscope for identification of allochems, textural classification (Embry et al. 1971), depositional facies and environments (Flügel 2013), as well as lithology and diagenetic processes. More than 1460 core plug samples, from three recently drilled wells, were used for routine core analysis (RCAL) parameters including porosity, permeability and grain density. These data were used in order to conduct reservoir rock typing and reservoir quality assessment. Primary drainage mercury injection capillary pressure data of 211 core plug samples, measured using CMS-300™ apparatus, were used in this study. In all samples, pertinent helium porosity and permeability (corrected for the Klinkenberg effect) are measured at overburden pressure condition. Conventional full set of wireline logs including spectral gamma ray, density, neutron, sonic and resistivity logs was available wells. Petrophysical logs were evaluated using optimizing petrophysical approach in this study. Results of petrophysical evaluation were used to calculate porosity in uncored wells and water saturation modelling.
In this study, reservoir rocks were classified into specific pore types using combination of petrographic studies and different rock typing methods. The identified rock types were used as a good indicator of reservoir quality and pore geometry index. In this study, we compare the results from and Winland r35 (Kolodzie Jr et al. 1980), Pittman r35 (Pittman 1992), FZI* (Mirzaei-Paiaman et al. 2018), FZI (Amaefule et al. 1993) techniques in defining rock types in our samples using the following equations: where r35, FZI and FZI* are in micron and permeability (k) and porosity ( ) are in mD and fraction, respectively. Moreover, we used Lucia's petrophysical classes according to the following equation (Lucia et al. 2001): where A = 9.7982, B = 12.0838, C = 8.6711, D = 8.2965 and RFN is Rock Fabric Number. Lucia defined rock classes according to the RFN values. The coarser grained rocks have a lower rock fabric number (RFN) and finer grained rocks a higher RFN. It should be noted that in Lucia's rock typing method similar RFNs correspond to rocks of similar grain size not pore throat radii.
Finally, capillary pressure curves of defined rock types were used to model water saturation by implementing parameters fitting method.

Petrography
Petrographic analysis of the Fahliyan Formation has resulted in the identification of six microfacies based on their allochems, textures and sedimentary environments. Facies characteristics of the studied samples are summarized in Table 1, and the representative photomicrographs are provided in Fig. 3. In terms of depositional setting, the studied interval was deposited in a homoclinal carbonate ramp represented by lagoon, shoal, algal mound and open marine facies ( Table 1). The recognized depositional model is consistent with previous studies on this formation in the Zagros area (e.g. Adabi et al. 2010;Jamalian et al. 2011;Jamalian et al. 2014;Dehkar et al. 2018;Noori et al. 2019). According to the sedimentological log ( Fig. 4) and petrophysical evaluation, this formation is substantially composed of limestone with some thin dolomitized intervals. As shown in Fig. 5, micritization, cementation, dolomitization, dissolution and compaction (chemical and mechanical) have been identified as the main diagenetic features of the studied intervals.
( . 2 a Location map of the Middle East, the orange rectangular is representing study area. (b) Stratigraphic column of the Early Cretaceous strata in south of Iran and surrounding areas (adapted from petroleum geo-service, 2001; http:// www. pgs. com) ◂ Considering facies characteristics, the LF and MF members are mainly consisting of grain-dominated facies with low mud content representing deposition under high-energy condition in shoal setting or algal facies with high values of porosity and permeability (Fig. 4). The sedimentary texture in these intervals varies from packstone to rudstone and floatstone.
In the MF, dominant pore types are interparticle, intercrystalline, intraparticle and mouldic. In the LF, in addition to all abovementioned pore types, dissolution moulds and vugs are frequently recorded that resulted in porosity and permeability enhancement (Fig. 4). Dissolution mainly occurred in shoal facies and led to enlargement of primary pores. However, calcite cementation occluded pore spaces and pore throats and negatively affected the reservoir quality in the studied succession.
The UF member has mainly composed of mud-dominated facies deposited in lagoon setting. In the lower parts, this member includes moderately sorted peloidal packstone with relatively high porosity, which are strongly influenced by calcite cementation (Fig. 4). Toward the upper parts, mudstone becomes the dominant texture. Some mouldic and intercrystalline pores are defined in this part of the formation.

Rock typing
Rock typing in hydrocarbon-bearing intervals, mainly carbonate rocks, has always been a challenging task for geoscientists and engineers. Lack of a universal approach for reservoir rock typing has impelled many researchers to introduce different methods (Kolodzie Jr et al. 1980;Pittman 1992;Lucia et al. 2001;Rezaee et al. 2006;Mirzaei-Paiaman et al. 2015, 2018. Moreover, data limitation and more importantly heterogeneous nature of carbonate reservoir in each hydrocarbon field impose much more difficulties in implementing these standard methods (Kadkhodaie-Ilkhchi et al. 2018;Al-Jawad et al. 2020). Carbonate reservoir rock types are generally complex to identify due to their formation process and sensitivity to diagenesis. Diagenetic processes such as cementation, dissolution, fracturing, dolomitization influence pore structure in carbonate rocks and increase pore network complexity and heterogeneity. As mentioned in the previous section, couple imprint of sedimentary environment and diagenetic processes in the studied formation have caused various pore types and considerable fluctuation of petrophysical properties. In the studied interval, porosity varies from less 1% to more than 26% and permeability ranges from 0.01 mD to more than 2000 mD. Figure 6 shows the cross-plots of porosity versus permeability and the Lorenz plots for horizontal and vertical plug samples. As it can be seen from this figure, the Lorenz coefficients, as an index of reservoir heterogeneity, are equal to 0.89 for horizontal samples and 0.82 for vertical plugs. This index can vary from zero in a homogeneous reservoir to a maximum value of 1 (Tiab et al. 2015;Mirzaei-Paiaman et al. 2019).
In order to find best rock typing method in this reservoir, we compare the results of FZI, r35Pittman, r35Winland, FZI* and Lucia techniques in defining petrophysical rock types in our samples. However, prior to any calculation, in order to exclude erroneous values, a cut-off was applied on the data to define nonreservoir intervals with very low porosity (< 3%) and permeability (< 0.01 mD) (Tiab et al. 2015;Mehrabi et al. 2019). Note that these very low Poro-Perm zones are classified as RT1 in this study. Different methods such as histograms, discrete rock  deh et al. 1996;Corbett et al. 2004;Tiab et al. 2015). In this study, DRT method was used because of its simplicity (Abbaszadeh et al. 1996;Mirzaei-Paiaman et al. 2018).
where A is a constant and should be adjusted for each index (i.e. FZI, Winland and Pittman's r35 and FZI*) in order to have the outcomes starting from 1. Including the rock type representing nonreservoir intervals with very low porosity (< 3%) and permeability (< 0.01 mD), using A equal to 2, 2.3, 2.9 and 3.2 for FZI, r35 Winland, r35 Pittman and FZI*, respectively, our data were classified into four rock types. Furthermore, we used Lucia's petrophysical rock types which classify our data into six rock types (Fig. 7).
In this study, MICP curves and pore throat radius distributions were used in order to evaluate the ability of abovementioned indices. Figure 7 shows the results of applying mentioned indices on our data. As it can be seen from this figure, while capillary pressure curves are not following similar trends, FZI cannot differentiate among them. On the other hand, it is expected that each rock type must have similar pore throat size distributions (Mirzaei-Paiaman et al. 2015, 2018. But this index classifies samples with completely different pore throat radii into same rock types. Also, implementing Lucia's petrophysical classes on our MICP data was not a good idea and this method's efficiency in separating samples into distinct groups is questionable (Fig. 7). Using r35winland and r35 Pittman, studied samples were classified into four rock types. Results show that both indices behave in relatively similar way and better than previously mention FZI and Lucia's methods (Fig. 7). Outcomes show that FZI* performed better than other indices to separate any existing rock types (Fig. 7). This method not only does find more distinct pore throat clusters than other techniques but also groups MICP curves in much more efficient way (Fig. 7). It should be noted that Mirzaei-Paiaman et al. 2015 and 2019 introduced FZI* method and verified it using 34 and 10 MICP samples, respectively. Therefore, the present study which include 211 primary drainage mercury injection capillary pressure data can be considered as a comprehensive study in verification of this valuable newly developed method in defining petrophysical rock types.

Pore facies (MICP)
Pore facies analysis is a practical approach for the reservoir rock typing according to the pore geometry characteristics (Chehrazi et al. 2011;Aliakbardoust et al. 2013). Pore system characterization plays an important role in hydrocarbon reservoir studies because it controls fluid saturation and the reservoir rock quality (Rahimpour-Bonab et al. 2014). Analysis of pore size distribution allows reservoir characterization to be pore system oriented, and therefore, the defined reservoir models will be based on the pore system characteristics (Chehrazi et al. 2011). The mercury injection capillary pressure (MICP) curves are popular tools for pore size distribution characterization in reservoir rocks and have been widely used in numerous case studies (Purcell et al. 1949;Dullien et al. 1974;Schowalter 1979;Hollis et al. 2010;Chehrazi et al. 2011;Skalinski et al. 2013;Xu et al. 2013;Rahimpour-Bonab et al. 2014). In this study, MICP curves were used to classify the carbonate reservoir rock samples based on their pore size distribution characteristics.

Pore facies 1
This rock type is entirely composed of lime mudstone and wackestone with no or very rare macroscopic pores. Almost all pores and microfractures are completely filled with sparry calcite cement (Fig. 8). It shows very low porosity (less than 3%) and permeability (less than 0.01 mD). Figure 8 illustrates an example of thin section image of this rock type and pore throat size distribution (PTSD) for all samples categorized in this class. As it can be seen from this figure, pores are in nanoscale with throat radius smaller than 0.1 micron, this increase entry pressure of this rock type to more than 2000 psi. This rock type mostly belongs to the upper parts of UF member.

Pore facies 2
This rock type has a grain-supported texture that is composed of fine (100-150 microns) peloids and ooids as dominant allochems (Fig. 8). Grains are well-rounded, often spherical and mostly very well sorted. Cementation is major diagenetic process in this rock type, occurring as pore-filling cements. Interparticle porosity was considerably occluded by calcite cements. Remaining porosity is because of the isolated moulds with poor connectivity and, consequently, low permeability. Average porosity and permeability are 12% and 0.4 mD, respectively. In this rock type, dominant pore throat sizes are between 0.1 and 0.5 microns which can classify as micropores. Average entry pressure is 500 psi (Fig. 8). It mainly belongs to the UF and MF members.

Pore facies 3
It represents grain-supported facies, containing fine (150-250 microns) peloids and ooids, as dominant allochems. Grains are well-rounded, often spherical and mostly very well-sorted. Calcite cementation and compaction are two important diagenetic processes in this rock type (Fig. 8). These has reduced pore throats and, consequently, permeability to some extent. Average porosity and permeability are 20% and 55 mD, respectively. Dominant PTS is between 0.1-2 microns. The entry pressure is below 250 psi (Fig. 8). This rock type shows limited distribution and can be detected mostly in the MF member.

Pore facies 4
This rock type is composed of coarse intraclasts grainstone, with large mollusk debris (Fig. 8). Equant calcite cement is main diagenetic feature in this rock type. Porosity and permeability are generally high; porosity values range from 11 to 26% and permeability values vary from 129 to 2100 mD. These high values of porosity and permeability result from the extensive dissolution. Pores dominantly macro-/ megapores and their throat radius range between 1 and 100 microns. However, different PTS are recorded indicating a poor sorting. Large pores with wide pore throat radii cause very low entry pressure (lower than 10 psi) and flat PC curves. This rock type mainly belongs to the LF member.

Saturation height modelling
Capillary pressure (PC) data are crucial in reservoir characterization and distribution during static and dynamic modelling (Tiab et al. 2015;Al-Bulushi et al. 2019). Capillary pressure curve tends to be more complex in carbonates in comparison with the siliciclastic reservoirs (Cuddy et al. 1993;Al Waili 2009;Chudi et al. 2010). Therefore, accurate determination of this property is essential in carbonate reservoirs characterization (Ahr 2011). There are a series of corrections including closure, stress, clay content and fluid corrections which should be applied on laboratory capillary pressure curves before fitting the capillary curves by the saturation height functions (Honarpour et al. 2004;McPhee et al. 2015;Lian et al. 2016;Kundu et al. 2017;Al-Bulushi et al. 2019). In this study, all PC curves are corrected for closure effect. Clay content and stress corrections were not applied because of clean (clay-free) limestone lithology of the Fahliyan Formation and laboratory test under reservoir P-T condition (Lian et al. 2016;Hulea et al. 2018). The last correction is wettability and interfacial tension correction. Because of the difference in wettability and surface tension between laboratory and reservoir fluid system, it is necessary to convert laboratory-derived PC data into reservoir PC (Lian et al. 2016;Kundu et al. 2017).
In this study, we used determined pore types as a basis to estimate water saturation. In order to decrease the uncertainty coming from permeability prediction and subsequent rock typing in uncored wells, parameters fitting approach were applied on capillary pressure curves. Defined pore facies were used to average MICP samples to have a representative capillary pressure curve for each pore facies. Then, we used Eq. (7) to extract fitting parameters in each pore facies as follows (Skelt et al. 1995): where PC is the capillary pressure in PSI, SW is water saturation in fraction and M and C are constant and curve fitting parameters. Using this equation, best fit values of the M and C parameters were determined for capillary pressure curves. Finally, relationship between core porosity and these parameters were investigated. Figure 9 represents relationship between fitting parameters and core porosity.
Substituting extracted correlations in Eq. (7) allows creating a model as a function of two variables: the capillary pressure (PC) or height above free water level (HAFWL) and the core porosity as follows: As porosity is the most straightforward parameter to be calculated from petrophysical curves with low uncertainty, it can be utilized in Eq. (8) to calculate water saturation in well and reservoir scale. This equation, in fact, belongs to a surface that describes the relationship between irreducible water saturation, core porosity and height above free water level in a 3D space. Figure 10 shows the model surface in relation to the capillary pressure curve for each pore type. Using equation of this surface, water saturation can be calculated at any point above free water level by porosity values derived from petrophysical curves.
In order to investigate the accuracy of the model, calculated water saturation was compared with log-derived water saturation in five of the oldest drilled wells and their petrophysical logs had not experienced production in the studied field (Fig. 11). As it can be seen in Fig. 11, there is a high consistency between the results. Figure 12 represents histogram of difference between log-derived water saturation and saturation height function-calculated initial water saturation in some of the oldest drilled wells in the studied field. It shows that error of calculation is dominantly very low and varies between ∓20 percent , with the highest frequency close to zero. Therefore, 3D model of water saturation can be, easily, built, based on the 3D model of porosity.

Discussion
In this study, sedimentary environment and diagenetic processes of the Fahliyan Formation were investigated in an oil field in northwest of Persian Gulf. Pore types were analysed based on detailed petrographic studies. Results show a wide various types of pore spaces ranging from micropores plugged with calcite cement to interconnected megapores which make a precise rock typing necessary for any kind of study in this reservoir. Rock typing analysis was performed using different methods to find the parameters controlling the reservoir quality and introduce the best way for classification of reservoir rocks according to their pore facies characteristics. Four pore facies (PF) were identified according to the FZI* method. The pore size distribution and capillary pressure curve for each pore facies is unique, and each pore facies contains samples with the greatest similarities in MICP and PSD curves. Since water saturation is dependent on the pore type and pore facies, the extracted pore facies were used in order to define the water saturation in the studied field. Pore facies were used to model initial fluid saturation using definition of a saturation height function that is written based on HAFWL and porosity.
Generally, reservoir quality increases from PF1 to PF4. This is evident from thin section studies, pore throat size distribution and entry pressure values for each facies. According to the petrographic analysis and MICP curves, a decrease in irreducible water saturation and entry pressure and increase in pore throat size in capillary pressure-defined classes are compatible with increase in interparticle to mouldic and vuggy pores ratio and a general decrease in the intergranular and porefilling cements from PF1 to PF4.
The results show that the initial rock texture is important controlling factor on reservoir quality. In fact, deteriorating diagenetic overprints act mostly in low-quality textures like mud-dominated intervals. On the other hand, grain-supported textures mostly experienced reservoir quality enhancer events such as dissolution. For example, PF1, which is easily detectable with high entry pressure, belongs to the lagoon environment. In this facies, pore types are dominantly nanopores and larger pores are fully filled by calcite cements. On the contrary, PF4 shows the best reservoir quality in the Fahliyan Formation. This is because of the coarse grain-supported textures and dominant interparticle pores, and extensive dissolution created a connected pore network with various pore throat sizes. The frequency of pore-filling cements decreases from PF1 to PF4. However, isopachous calcite cement is common only in PF4. Pore throat radius changes from nanoscale to megascale from PF1 to PF4. PF1 is limited to the UF member. Formation in the low-energy environment and subsequent cementation in large pore spaces caused low porosity and high percentage of irreducible water saturation intervals in upper parts of Fahliyan Formation. An increase in pore throat radius size to microand mesopores in PF2 and PF3, respectively, resulted in decline in irreducible water saturation and entry pressure rather than PF1. In spite of the preserved porosity, permeability reduced in PF3 due to the diagenetic processes, mainly physical compaction and cementation. PF4 with macro-to megapores is the most frequent facies in LF member. Interconnected pore spaces resulted from dissolution, enhancing the porosity and permeability, created the best reservoir intervals in LF member with flat MICP curves and low irreducible water saturation.
As permeability prediction and subsequent rock typing in uncored wells are challenging tasks in reservoir studies, most specifically in heterogeneous carbonate reservoirs, relationship between pore types and irreducible water saturation was used to model water initial saturation in the pore facies framework. Generally, from PF-1 to PF4 irreducible water saturation decreases, while pore types change from isolated pores to completely interconnected pore spaces with the highest reservoir quality.

Fig. 11
Comparison between the calculated water saturation using petrophysical logs-and saturation height function-derived water saturation in five of the oldest wells drilled in the studied field Results show that a there is a logarithmic correlation between porosity and fitting parameters that can be utilized to generate a formula to model water saturation. From PF-1 to PF-4, M values increase, while C decreases from PF-1 to PF-4. Results show that extracted relationship between porosity and fitting parameters in core scale can be used to model water saturation in well or reservoir scale using log-derived porosity.

Summary and conclusions
In this study, different types of information including thin sections, routine and special core analysis results were integrated to categorize reservoir intervals of the Fahliyan Formation into specific rock types and model the initial water saturation in the framework of pore facies. The following conclusions were obtained through application and consideration of facies studies, diagenesis, reservoir properties and MICP assessments.
• According to the core and thin-section petrographic studies, six facies were recognized in the Fahliyan Formation. Mudstone to wackestone facies deposited in lagoonal settings have the lowest reservoir quality. The main pore types in these intervals are microporosity and interparticle. On the other hand, facies with grainstone to rudstone textures deposited in shoal environment have best reservoir quality formed by dominant interparticle and vuggy pores.
• Petrophysical rock types were defined by implementing different rock typing indices on primary drainage capillary pressure data. In this way, four petrophysical rock types were recognized based on pore size distribution and capillary pressure curves trends. FZI* method is the best way to separate data set into discrete groups in the way that each class represents unique pore geometry and capillary pressure curve (pore facies). • The identified pore facies were used as basis for water saturation modelling. Porosity as the straightforward reservoir properties to model was used to model water saturation.
Acknowledgments The authors wish to thank the NIOC for the data preparation and permission to publish the data. The authors are also grateful to Turgay Ertekin, Editor-in-Chief of Journal of Petroleum Exploration and Production Technology, for editorial handling and three anonymous reviewers of this article for their time and constructive comments.

Conflict of interest
On behalf of all the co-authors, the corresponding author states that there is no conflict of interest.

Funding No funding.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, Frequency ΔSW (v/v) provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.