Received 26 April 2016; accepted 25 July 2016; published 28 July 2016
More than 65 percent of proven petroleum resources in the world and 34 percent of gas fields of the world are located in giant hydrocarbonic fields in the Middle East   . Zagros is one of the most important fold belts and active tectonic in the Middle East that is suited in South-West of Iran, between the Arabian Plate and Center of Iran    . This area containing lucrative reservoir bodies and including enormous hydrocarbonic fields is famous as one of the richest oil belts in the world   . Due to mud-supported nature of more facies which inherited low reservoir potential from their depositional settings, diagenetic process plays an important role in increasing of reservoir quality, but not always. However diagenetic processes have an important control on reservoir characteristic, in which even the quality of reservoirs can be investigated by studying of the diagenesis processes  -  . In fact, various diagenesis processes are able to influence on petrophysics features with different intensity on total and effective porosity, permeability, size of pore type and their distribution for example    -  . Therefore investigations of diagenesis processes are necessary for understanding reservoirs quality. The goal of this study is survey of effect of diagenetic presses on reservoir properties and also flow unites determination in the Ilam Formation of Dezfol embayment using petrography evidence and petrophysics data (such as   , and many others).
2. Geological Setting and Stratigraphy
The Zagros belt is extended along North-Northwest boundary of Arabian Plate, i.e. from Southeast of Turkey, Syria and Iraq to the Minab’s fault in South of Iran     . Dezfol embayment in Zagros sedimentary-structure zone is surrounded via tree structural phenomenon’s, which are included Flexural Balarood zone in Northwest, flexural Kohestan zone in Northeast and a complicated fault zone named Qatar-Kazerun in East and Southeast (Figure 1). More than 45 fields are discovered inside Dezfol embayment  .
The cretaceous deposits of Dezfol embayment such as Zagros embraced significant amount of hydrocarbonic
Figure 1. Approximate position of Dezfol embayment (updated from  ).
reserves.  -  suggested the sedimentary carbonate ramps environment was extended over long time in this area. They also emphasized that the Arabian Plate, in that period of time, has tropical and subtropical climate. During the late Campanian and late Santonian, seawater had covered Zagros fields, and Ilam Formation was deposited. Then relative sea level rise continued and Gurpi Formation which have shale and marl lithology was deposited  . Wynd and Setudehnia (1975) believe that the boundary between two Formations of Ilam and Gurpi formed gradually and continues. Whereas Furst (1970) considered an erosional unconformity nature for mentioned boundary  . The stratigraphy position of this formation on the Arabian Plate and Dezfol embayment has been showed in (Figure 2).
3. Data and Method
In this research, 300 thin sections of cores and cutting from five wells in the studied field have been studied. Polarizing microscope is used to investigate petrography and diagenesis evidences and to describe facies. A modified Dunham (1962)  texture scheme was used. Flugel (2010)  method is applied for micro-facies determination and assortment, and finally method of Choquette and Pray (1970)  is used for classification and ordering of the facies. In order to recognition of depositional settings of facies associations and reconstruction of palaeoenvironmental condition, standard models of ancient and recent carbonate environments are used as analogues    . Detailed petrography analysis and well-logging data which are provided by Geolog professional software, this Formation is divided into five different zones of reservoir units. In term of reservoir quality. Determined reservoir unite are correlated among all five well and flow units have been determined across the studied field.
4. Facies and Depositional Environment
Detailed petrographic analyses, have led to identification of 10 Micro-facies types, which are grouped into six facies associations (Table 1). Depositional environments of facies range from the outer ramp to inner ramp settings. In order to recognition of depositional settings of facies associations and reconstruction of palaeoenvironmental condition, standard models of ancient and recent carbonate-evaporate environments are used as analogues     (Figure 3).
5. Diagenetic Processes
The main diagenetic alterations and their controls on the reservoir properties of Ilam formation are elaborated in details and presented as follow and have been showed in Figure 4.
Figure 2. Ilam formation and its stratigraphy position at Arabian plate (modified from   ).
Table 1. Determined facies and depositional setting of Ilam formation of the studied field.
Figure 3. Conceptual depositional models of Ilam formation of the studied field.
In the studied interval, early stage of marine cements is very rare, but the late stage diagenetic cements are frequent. Based on the time of formation, the following types of cement are recognized.
5.1.1. Isopachous Bladed Cement
This type of cement created in form of a thin and isopachous layer around the grains (Figure 4(a)). Existence of the layer of cement between the grain’s boundaries indicates that these layers are generated on early stage of d genesis as first generation of cement, which mainly observed in the grainstone and packstone facies  . As a result of physical pressure, this type of cement prevents from the deformation and keeps the primary porosities  which meaning positive effect on reservoir quality. This cement is the index of sea water phreatic zone    .
5.1.2. Equant Cement
Approximate homologous dimension is the remarkable common characteristics of this type of cement, which crystals are subhedral to anhedral form and boundary between each crystal is distinctly visible   (Figure 4(b) and Figure 4(c)). This cement type is observed within grainstone facies with moldic porosity. Dissolution of unstable fragments created moldic porosity. The materials that are produced from dissoluted particle formed
Figure 4. Main detected diagenetic processes in the Ilam formation.
Equant Cement and fill the space between of grains. This type of cement mainly observed in form of moderate to large size crystals, and is able to develop in various diagenesis environments such as meteoric to burial  .
5.1.3. Poikilotopic Cement
This type of cement mainly exists in large size because it contains several crystals (Figure 4(d)). Ahmad and Bhat (2006)  believed, the formation of this cement type is based on low-speed grow of calcite crystals during burial environment and can be filled free space  . This pervasive and widespread cement has filled inter-crystal and inter-grain’s porosity and caused of significant decline in porosity and permeability  . The mentioned cements are coarse in size due to low-speed growth of crystal  .
5.1.4. Drusy Cement
In this particular cement, dimension increment occurred toward the inter-grain and between grains centers. It also created in both burial and meteoric environments near the earth level  . This cement fills the inter-grain and between-grains spaces and caused of a significant decline in porosity and permeability in grain supported facies. However, it is formed in shallow burial which mostly have close relationship with meteoric diagenesis  .
5.1.5. Syntaxial Cement
One of the most important cements in the studied interval is Syntaxial cement, and is a witness to the free-water involvement  . The cement accumulates around a singular grain mostly echinoderm with composition of HMC. Vados, metoric and burial environments are favorable for formation of this cement. This cement with dark and cloudy appearance is created in marine diagenesis setting and cements with clear and bright crystal are related to burial and meteoric beds  .
During this process, primary crystals consumed and will replace by new crystals, which have a similar chemical combination  -  . As  represents, the process mainly occurs in present of water, diagenesis and wet environments, where dissolution and sedimentation are happening. In the other words, neomorphism is related to the growing some of crystals costs due to losing some other crystals and according to the  required calcium carbonate will provided from dissolution of small crystals and inter-grain’s water. In the studied interval neomorphism has important effect on porosity and subsequent reservoir quality.
 states that this process can take place as a result of consecutive activities of microorganisms such as, alga, fungus and cyanobacteria, that are exist on allocems surface. In Ilam Formation, micritization process was developed around the skeleton fragments. Although in some cases, the internal structure was loosed as a result of the intensity of micritization, which can be due to the high activities of microorganisms and inadequate water rotation   . In many cases, intensive bacterial activities destroyed entire internal structures of lime fragments (Figure 4(b) and Figure 4(d)).
Stylolite is rough and ridged surface contact between the rock mass which mainly created in more 500 meter depth  . In the some space with high content of clay mineral, stylolite and dissolution seams are more observable   . Stylolitization process can have negative or positive effects on reservoir quality. Stylolite acted as a blocking force against fluid flow and results declining of reservoir quality. In other hand, stylolite can enhance the reservoir quality on the some condition with passing of the flow and consequently in fact increasing permeability by connection of pore type  .
Dolomite occurred in Ilam Formation in singular and scattered rhombus forms but with low frequency. However, they have centralized around the stylolites and dissolution seams. Dolomite crystal centralization around the stylolites and disbandment channels demonstrates that they are generated from fluid with pressure dissolution in last stage of burial process see  . In fact, pressure dissolution caused of required magnesium concentration  .  recognized that dolomitization is more likely to have positive effect on the quality of reservoir. According to the  and  studies, the lime mud-which composed of more thickness of interval-Played a role as main supply of magnesium for dolomitization. Scattered dolomite crystals along the stylolite and pressure dissolution seams can be seen in some sample; therefore it is justified for formation term of burial model dolomitization  . The dolomites crystal with size ranged from tiny to middle presumed probably origin from hydrocarbon dissolution during migrating fluid when stylolite were active channels for flows see   .
Dissolution is an effective diagenesis process which mainly caused of secondary porosity such as moldic, vuggy and channel type. Dissolution process can take place in meteoric, marine and burial environments  . Microscopic dissolution and its consequence porosity will be discussed in details in porosity section.
Different types of visual porosity observed in Ilam Formation are described in ongoing part.
5.7.1. Moldic Porosity
Selective dissolution happened among aragonite grains of limestone and remains empty space in the grain place which are similar to the primary grains shape   (Figure 4(b)). This porosity is fabric selective, which mostly resulted during the meteoric and burial diagenesis environments   .
5.7.2. Intraparticle Porosity
Intraparticle porosity is type of porosity which has low frequency in the studied interval. In some sample the inter-grain pores may not be connected together appropriately.  represented which often this porosity formed during meteoric diagenesis.
5.7.3. Intercrystaline Porosity
Mainly, intercrystaline porosity not preserved among carbonate rocks, because these are considerably sensitive to the burial pressure and primary porosity was destroyed  . But in some cases this porosity was remained, especially when primary cementations are occurred.
5.7.4. Vuggy Porosity
This type of porosity observed in large forms (more than 1.16 mm). Vuggy porosity is nonselective fabric and has larger that grains size and mostly results of meteoric dissolution  . The most common porosity among the studied wells is vuggy porosity (Figure 4(e)).
5.7.5. Fracture Porosity
Tectonic forces are main agent for occurring of this porosity. Fracture porosity in carbonate rocks is a common feature and happened during burial diagenetic stage  .
6. Reservoir Characterization and Flow Units
Different parts of Ilam Formation are studied by petrography and petrophysics data in order to investigate the Reservoir characterization and flow unit determination. Analyzing the well logging that is achieved from petrophysic technique are used in order to determine the quality of reservoir in different zones, and additionally for determining the best zones for developing the oil fields  . Reservoir quality determined by calculating the porosity, and verifying the lithology type  . According to the  classification, Ilam’s reservoir is a diagenetic reservoir. Since, dissolution is the main key in increasing the quality of diagenetic reservoir. Identifying the reservoir zones from non-reservoir zones (reservoir unit determination) and subsequent flow unit determination is an important step for the field developing plans. There are different parameters for identifying the reservoir zones from non-reservoir zones such as, analyzing the common petrophysics logs and using petrology features which all was used in this study and represented in Figure 5. Total porosity of these zones evaluated via Geolog software and presented in Table 2 which are according to the  standard. Also, porosity distribution over the studied field is depicted by Surfer software (Figure 6). Two and three-dimension models of distribution demonstrates that highest amount of porosity occurred in Northwest of the field (Figure 6).
Detailed petrographic analyses, have led to identification of 10 micro-facies which are represented as a carbonate ramp depositional model. According to the studies, Ilam Formation endured three different diagenetic environments: meteoric, marine and burial diagenesis. Various digenesis processes have been recorded for this interval,
Figure 5. Flow unit determination across the studied field.
Figure 6. (a) Two-dimension model of porosity distribution over studied field; (b) Three-dimension model of porosity distribution over studied field.
Table 2. Average porosity and reservoir quality of determined flow unit.
such as dissolution, cementation, micritization, neomorphism and dolomitization. From these processes, only dissolution created reservoir quality by producing different types of porosity, especially vuggy types. This study revealed that Ilam’s reservoir is a diagenetic reservoir in the studied field. In order to achieve more reliable view from reservoir quality evaluations and also determination of reservoir unit using petrography and petrophysics data, the mentioned reservoir was divided into five units. Evidence shows that in unit-I1 and I2 there are significant numbers of mud supported facies, approximately no diagenesis processes recorded in these zones, which lead to the low porosity and poor reservoir qualities. Unit-I3 includes clear lime and has the highest amount of porosity presenting good quality, which endures more dissolution. Unit-I4 and I5 have moderate porosities and their reservoir qualities are fair to good. Also, according to the two- and three-dimension models of porosity distribution over the 5 units in horizontal extension, it can be stated that more percentage of porosity occurred in northwest of field.
 Sharland, P.R., Archer, R., Casey, D.M., Davies, R.B., Hall, S.H., Heward, A.P., Horbury, A.D. and Simmons, M.D. (2001) Arabian Plate Sequence Stratigraphy: GeoArabia. Special Publication, 2, 371 p.
 Talebian, M. and Jackson, J. (2002) Offset on the Main Recent Fault of NW Iran and Implication for the Late Cenozoic Tectonics of the Arabia-Eurasia Collision Zone. Geophysical Journal International, 150, 422-439.
 Mazzullo, S.J. (1994) Diagenesis in a Sequence-Stratigraphic Setting: Porosity Evaluation in Periplatform Carbonate Reservoirs, Permian Basin, Texas and New Mexico. Journal of Petroleum Science and Engineering, 11, 311-322.
 Heydari, E. (1997) The Role of Burial Diagenesis in Hydrocarbon Destruction and H2S Accumulation, Upper Jurassic Smackover Formation, Black Creek Field, Mississippi. American Association of Petroleum Geologists Bulletin, 81, 26-45.
 Zhang, J., Qin, L. and Zhang, Z. (2008) Depositional Facies, Diagenesis and Their Impact on the Reservoir Quality of Silurian Sandstones from Tazhong Area in Central Tarim Basin, Western China. Journal of Asian Earth Sciences, 33, 42-60.
 Rahimpour-Bonab, H., Asadi-Eskandar, A. and Sonei, A. (2009) Controls of Permian-Triassic Boundary over Reservoir Characteristics of South Pars Gas Field, Persian Gulf. Geological Journal, 44, 341-364.
 Baron, M., Parnell, J., Mark, D., Carr, A., Przyjalgowski, M. and Feely, M. (2008) Evolution of Hydrocarbon Migrationstyle in a Fractured Reservoir Deduced from Fluid Inclusion Data, Clair Field, West of Shetland, UK. Marine and Petroleum Geology, 25, 153-172.
 Cerepi, A., Barde, J.P. and Labat, N. (2003) High-Resolution Characterization and Integrated Study of a Reservoir Formation: The Danian Carbonate Platform in the Aquitaine Basin (France). Marine Petroleum Geology, 20, 1161-1183.
 Rahimpour-Bonab, H., Mehrabi, H., Enayati-Bidgoli, A.H. and Omidvar, M. (2012) Coupled Imprints of Tropical Climate and Recurring Emergence on Reservoir Evolution of a Mid-Cretaceous Carbonate Ramp, Zagros Basin, Southwest Iran. Cretaceous Research, 37, 15-34.
 Al-Laboun, A., Al-Quraishi, A., Zaman, H. and Benaafi, M. (2014) Reservoir Characterization of the Burqan Formation Sandstone from Midyan Basin, Northwestern Saudi Arabia. Turkish Journal of Earth Sciences, 23, 204-214.
 Bahroudi, A. and Koyi, H.A. (2004) Tectono-Sedimentary Framework of the Gachsaran Formation in the Zagros Foreland Basin. Marine and Petroleum Geology, 21, 1295-1310.
 Bordenave, M.L. and Hegre, J.A. (2005) The Influence of Tectonics on the Entrapment of oil in the Dezful Embayment, Zagros Fold Belt, Iran. Journal of Petroleum Geology, 28, 339-368.
 Alsharhan, A.S., Nairn, A.E.M. and Mohammed, A.A. (1993) Late Palaeozoic Glacial Sediments of the Southern Arabian Peninsula: Their Lithofacies and Hydrocarbon Potential. Marine and Petroleum Geology, 10, 71-78.
 Ghabeishavi, A., Vaziri-Moghaddam, H. and Taheri, A. (2009) Facies Distribution and Sequence Stratigraphy of the Coniacian-Santonian Succession of the Bangestan Palaeo-High in the Bangestan Anticline, SW Iran. Facies, 55, 243-257.
 Holliss, C. (2011) Diagenetic Controls on Reservoir Properties of Carbonate Successions within the Albian-Turonian of the Arabian Plate. Petroleum Geoscience, 17, 223-241.
 Koop, W. and Stoneley, R. (1982) Subsidence History of the Middle East Zagros Basin, Permian to Recent. Philosophical Transactions of the Royal Society A, A305, 149-168.
 Aqrawi, A.A.M., Mahdi, T.A., Sherwani, G.H. and Horbury, A.D. (2010) Characterization of the Mid-Cretaceous Mishrif Reservoir of the Southern Mesopotamian Basin, Iraq. AAPG GEO Middle Eastn Geoscience Conference, Bahrain, 7-10 March 2010, AAPG Search and Discovery Article 50264.
 Dunham, R.J. (1962) Classification of Carbonate Rocks according to Depositional Texture. In: Ham, W.E., Ed, Classification of Carbonate Rocks, Vol. 1, American Association of Petroleum Geologists Memoir, 108-121.
 Ahmad, A.H.M., Bhat, G.M. and Azim Khan, M.H. (2006) Depositional Environments and Diagenesis of the Kuldharand Keera Dome Carbonates (Late Bathonian-Early Callovian) of Western India. Journal of Asian Earth Sciences, 27, 765-778.
 El-Saiy, A.K. and Jordan, B.R. (2007) Diagenetic Aspects of Tertiary Carbonates West of the Northern Oman Mountains, United Arab Emirates. Journal of Asian Earth Sciences, 31, 35-43.
 Kirmaci, M.Z. and Akdag, K. (2005) Origin of Dolomite in the Late Cretaceous-Paleocene Limestone Turbidites. Eastern Pontides, Turkey. Sedimentary Geology, 181, 39-57.
 Ehrenberg, S.N., Nadeau, P.H. and Aqrawi, A.A.M. (2007) A Comparison of Khuff and Arab Reservoir Potential throughout the Middle East. AAPG Bulletin, 91, 275-286.
 .Smith, J.V. (2000) Three-Dimensional Morphology and Connectivity of Stylolite Shape Reactivated during Veining. Journal of Structural Geology, 22, 59-64.
 Hood S.D., Nelson C.S. and Kamp P.J.J. (2004) Burial Dolomitisation in a Non-Tropical Carbonate Petroleum Reservoir: The Oligocene Tikorangi Formation, Taranaki Basin, New Zealand. Sedimentary Geology, 172, 117-138.