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.