and the quality of k4 reservoirs was studied through analyzing petrophysical data. In fact, the aim of this article is to study sedimentary formations, types of porosity and diagenetic control on evolution of porosity using petrophysical data and to investigate their effect on the quality of k9 reservoir zone in Kangan Formation, in one of the hydrocarbon fields in the south of Iran. To study more dimensions of this subject, the following hypotheses are presented.
1.1. The Main Hypothesis
It is feasible to study the effect of sedimentary facies, types of porosity and diagenetic evolution control on the quality of reservoirs using petro physical data in Kangan Formation in one of hydrocarbon fields in the south of
• It is feasible to study the effect of sedimentary facies on the quality of the reservoir in Kangan formation in one of hydrocarbon fields in the south of Iran using petrophysical data.
• It is feasible to study the effect of types of porosity on the quality of the reservoir in Kangan formation in one of hydrocarbon fields in the south of Iran using petrophysical data.
• It is feasible to study the effect of diagenetic evolution control on the quality of the reservoir in Kangan formation in one of hydrocarbon fields in the south of
So far, many studies have been carried out on sedimentary facies, diagenetic processes, type and geometry of pores in different formations such as Kangan Formation k4   . No comprehensive research has been done on the effect of sedimentary facies, types of porosity and diagenetic evolution controls on the reservoir quality especially in Kangan k4 using petrophysical data. So this article presents an innovative research. Kangan Formation is located in South Pars field of
Figure 1. The location of South Pars hydrocarbon Field in the Persian Gulf.
1.2. Theoretical Foundations
Due to the subject of this article, the following terms are explained:
Facies: This term was first introduced by reference  . In this definition, the facies are referred to the specific characteristics of a sedimentary unit. These characteristics include lithologic and fossil features. These properties include: color, layering, composition, texture, sedimentary structures and fossil elements. As there are significant differences between the fossil and lithologic characteristics, two types of facies i.e. lithofacies and biofacies are recognized  .
Diagenetic processes in Kangan Foramtion: Different diagenetic processes associated with atmospheric diagenetic environment, seafloor and burial diagenetic environment are seen in facies of Kangan Formation. The most important processes are dissolution, cementation, micritization by microorganisms, compaction, recrystallization, dolomitization and anhydritization  .
Reservoir Quality: By determining the hydraulic flows, it is possible to separate the reservoir sections from non-reservoirs. Each of these hydraulic flow units is defined based on the relationship between porosity and permeability of different reservoir sections and indicates the quality of the reservoir in that section  .
Sedimentary environment: The sedimentary environment is part of the earth where sedimentary materials are deposited. It is distinguished physically (depth, sedimentation rate, direction of water flow, specific weight, and water temperature), chemically (oxidation and reduction, pH of the environment, salinity) and biologically (existence or absence of organisms or fossils) from the adjacent environment  .
2. Literature Review
Reference  studied the facies, sedimentary environment and diagenetic evolution control on the quality of reservoir in Upper Permian Upper Dalan Formation in Kish Gas Field (Zagros offshore basin). Kish gas field in Zagros off shore basin is located between Iran and Qatar and is based on a 222 m thick continuous core. The quality of the reservoir was severely affected by the changes in the rock fabrics and subsequent diagenetic alterations. The types of porosity are reinter-particle, moldic and connected vug. The interparticle relationship between porosity and permeability shows that different crystal sizes do not affect the reservoir quality. That study shows that pervasive pore-filling anhydrite mineralization decrease significantly the porosity and permeability, but the size of pores is slightly affected by it. In 2016, the effect of sediment and diagenetic control on reservoir properties was examined in the Permian Triassic sequence in the west of the Persian Gulf. This study provides information on the profile of facies and diagenesis and the properties of reservoirs in the Golshan Field in the Western Persian Gulf. The results show the effect of sedimentary processes on the quality of reservoirs. Among the diagenetic processes, cementation, calcite and anhydrite have a significant destructive effect on reservoir quality in prefabricated porous facies. The results also suggest that there is a close relationship between sedimentary facies and diagenetic processes in the Kangan Dallansequence  . Reference  studied the diagenesis of k4 Kangan formation in south Pars field with an emphasis on the dolomitization and anhydritization. Diagenesis includes changes that have occurred since the sediment was deposited before the onset of metamorphism. K4 Kangan formation (Lower Triassic) is the reservoir rock of South Pars Gas Field which is located in Dehram Group. Dolomitization, anhydritization, micritization, dissolution and cementation are effective diagenetic processes in Kangan formation. Identified Dolomites in this formation include fine grained dolomites, medium grained dolomites, medium to coarse grained dolomites, coarse grained dolomites or cement dolomites and saddle dolomites. Poikilotopicanhydrite, anhydrite nodules, pore-filling and pervasive anhydrite cement, layered anhydrite, fracture-filling anhydrite and anhydrite veins are of various types of anhydrites in Kangan k4 formation. The porosities in this sequence include intra-particle porosity, inter-particle porosity, moldic porosity, vuggy porosity, crystal porosity, shelter porosity and fenestral porosity. In another study, the diagenesis and reservoir characteristics of k4 Kangan formation and the upper section of Dalan formation in one of gas fields of Persian Gulf was investigated. The gas field is located in Persian Gulf and in south east of port of Bushehr. To investigate iagenetic processes, petrographic analyses were conducted on 162 thin sections obtained from cores of Kangan formation and the upper part of Dalan formation in well A of X field. Dissolution, dolomitization, cementation, compaction and anhydritization were important diagenetic processes which affected the reservoir characteristics of under studied facies. Among those processes compaction, cementation and anhydritization reduced the quality of the reservoir but dissolution and dolomitization increased its quality  .
3. Statistical Population and Research Samples
In this study, the statistical population is one of the hydrocarbon fields in south of Iran. To do the study, 166 thin sections of two wells from core of Kangan k4 formation were used as diagenetic samples, 85 of which were used in experimental studies. After preparing the rock samples and microscopic thin sections, the plaques which were made from Kangan formation were studied.
In this research, the data were collected from reports, published theses, related texts, geographic and topographic maps, internet websites, library of Geological organization and Islamic Azad University. The data were analyzed using GIS software. By investigating the microscopic and macroscopic characteristics of thin sections of Kangan Formation (k4), the facies and sedimentary environment affecting the reservoir quality of the formation in two wells of one hydrocarbon field in south of Iran was studied. In microscopic study, 166 thin sections of two wells (3 × 2) with a core of one of the hydrocarbon fields of southern Iran are used. Using polarized microscopy, allochems, the major complications of diagenesis and facies were investigated and using Sibley and Gregg classification (1987), dolomite rock textures were classified. The main goal of petrophysical researches is to study the properties of rocks and their relation with fluids inside them. Petrophysical assessment is the science of interpreting the information collected from well logs. The first step in assessing petrophysical data is collecting digital data. The digital data is entered into the software and the logs are drawn. Then according to Schlumberger log interpretation charts book environment a land mud invasion corrections are done  . To interpret logs and to choose different petrophysical techniques, it is necessary to have information about the lithology of the formation, properties of the drilling fluid and the characteristics of fluid inside the rock. The aim of this article is to assess reservoir properties such as petrophysical parameters like porosity, permeability and saturating fluids, lithological variations and properties of hydrocarbon such as useful and useless thickness and useful hydrocarbon column. The quality of k4 Kangan reservoir is studied through analyzing petrophysical data. The main petrophysical charts such as neutron, density and sonic log data were used to determine porosity.
5.1. Determining Lithology
One of the main procedures in assessing reservoir properties is lithology identification which helps separate zones with reservoir properties from zones without these properties. Lithology identification is fulfilled using four crossplot techniques. One of the applications of Dipole Sonic Imager is lithology identification. In current Dipole Sonic Logging, two components are considered for solid parts. The first component is matrix which contains main aggregates and cement and the second is shale components. In lithology, some properties such as mineral texture, shale volume and liquid content may influence on the response of the Imagers. Texture contains parameters of texture components related to matrix and the characteristics of bed and the structure of the gap. Data obtained from cores, drilling and Well Logging Imagers helps to estimate the lithology of rocks exactly.
Shale volume calculation:
Equation (1) shows shale volume which was calculated using corrected gamma ray (CGR)  .
With this technique the minimum and maximum of CGR in understudied wells were calculated. In Equation (1) CGRmin is the minimum value at non-shale formations (clean) and CGRmax is the value that is read at a 100% shale formation. This equation does not have any unit. In this equation, CGR (corrected gamma ray) relates to the desired depth. Shale volume calculated in A and B wells are 1.88 and 1.12 percent respectively and the shale volume average in both wells is 1.54 percent. The desired well is considered non-shale and clean.
Petrophysical properties of reservoirs such as porosity are the most important controllers of quality and quantity performance of the reservoir. The calculation of porosity is done using neutron, density and sonic porosity logging. Neutron porosity logs are used to calculate porosity. This tool shows hydrogen frequency or index. Neutron porosity logs response to the volume of water filled the pores. Neutron porosity log is calibrated for limestone, so in calcareous lithology the value of neutron porosity log is equal to the amount of porosity but dolomitic sandstone lithology should be corrected for matrix. The amount of porosity is obtained by Equation (2) in which f is real porosity, a and b are experimental constant values and N is the value of used tools  .
In non-shale formations (clean), porosity is calculated by density logging through Equation (3)  . According to this equation, matrix and density of the fluid determine porosity.
Based on this equation, matrix density of formation was 2.9 gram per cubic centimeter and fluid density was 1 gram per cubic centimeter. In sonic porosity logging, difference in sound speed causes separation and recognition of rock properties in a way that the less density of the filling fluid, the less speed of sound. Figure 2 shows sonic graphs by which sonic porosity was calculated. In these graphs, solid lines on the chart are made using curves and Wiley formula related to Rimmer equation.
Figure 2. Sonic graphs for calculating porosity.
Porosity can be calculated using sonic chart and Equation (4)  (Fertle et al., 1987).
where DTma: the time in which the wave passes through formation matrix, 45 (ms/ft).
DTfl: the time in which the wave passes through the fluid, 185 (ms/ft) (Table 1).
Table 1. Calculating porosity by sonic logging.
In resistance logging technique, Archie equation was used to calculate porosity. This technique is used when the sample is saturated with water and porosity is calculated using Equation (5)  . In the following equation m and a are experimental constant values, R is the resistance of filter mud and Rxo is the washed resistance zone:
Determining lithology using Crossplot technique:
Usually four crossplot techniques are used to determine lithology:
1) Neutron-density crossplot.
2) Sonic-neutron crossplot.
3) M-N plot technique.
4) MID plot technique.
Among double charts, neutron-density crossplot is the best technique for separating different minerals. In this study, lithology, permeability and porosity of Kangan formation were calculated using this crossplot. Neutron-density cross- plot has maximum resolution and it can best determine porosity in different lithologies. In Neutron-density crossplot, three lines are drawn which are related to sandstone, limestone and dolomite. To calculate porosity using this chart, density response is plotted against neutron porosity. Based on the distance between the point and the matrix lines, the position of plotting shows the percentage of the lithology. Figure 3 and Figure 4 show neutron-density crossplot that are used to determine lithology in the wells. The horizontal axis shows porosity calculated through neutron logging in calcareous (limestone) lithology and the vertical axis shows porosity calculated by density logging in calcareous lithology. As neutron logging is used, density logging in limestone is read from the straight line which is related to limestone figure. The numbers on lithology lines show the porosity.
Figure 3. Neutron-density cross-plot of well A.
Figure 4. Neutron-density cross-plot of well A.
In sonic-neutron crossplot, the resolution between different lithologies is close to neutron-density crossplot. In graphic logging, it is enough to draw sonic plots against neutron. The plot position is a point which shows the distance from the point to matrix line and presents lithology percentage. These figures have low resolution in the presence of volatile minerals. Figure 5 and Figure 6 indicate sonic-neutron crossplot of understudied wells. The results confirm the lithology and porosities obtained from neutron-density crossplot.
Figure 5. Sonic-neutron crossplot in well A.
Figure 6. Sonic-neutron crossplot in well B.
For determining lithology, M-N crossplot is used with three porosity logs. M-N crossplot is actually used to eliminate the effect of porosity and to determine triple mineralogical composition. To draw this crossplot, M and N are drawn against each other. The value of M and N is calculated using  the following equation:
In this equation, two diaphragms are used for calculating M or N so the effect of porosity and lithology is eliminated to a large extent. Thus M and N are just the lithological functions. One of the main applications of M-N crossplot is determining secondary porosity. Figure 7 shows the M-N crossplot of under studied well. In M-N crossplot, when the points are plotted above calcite-do- lomite line, and the formation contains gas, it shows that there is a secondary porosity. In Kangan formation, the samples are above calcite and dolomite points (Northeast of crossplot), which indicates the formation contains gas. There is no photoelectric factor in well B, so no M-N crossplot is drawn for well B.
Figure 7. M-N crossplot of well A.
In MID crossplot, in order to determine the lithology, first the nominal values of matrix i.e. (∆tma), (fNma) and (ρma) should be determined. Neutron-density crossplot is used to calculate the nominal values of (ρma) matrix, and Sonic-neu- tron crossplot is used for calculating (∆tma). After determining these two
Figure 8. MID-plot for well A.
parameters, MID plot was drawn. Figure 8 shows MID plot for well A.
Using neutron-density crossplot, most of the under studied formation was evaluated and lithological components such as limestone and dolomite were recognized. In some parts of the formation, a low amount of shale was found, as well.
5.2. Calculating Porosity Using Double Logging
This technique does not need any primary information and the type of porosity and lithology are determined and calculated using neutron-density crossplot, sonic-neutron crossplot and sonic-density crossplot. In calculating through neutron-density crossplot, the porosity is determined based on the values of neutron log. In Figure 9, neutron log is drawn against the density log. As is indicated in Figure 9, the gas source reduced the value of neutron log and the density of the rock which resulted in the reduction of density value. Comparing neutron with density double logging, the result indicates that the presence of gas can separate these two logs which are called the effect of rotor. Using Equation (8), the porosity is calculated through neutron-density crossplot  (Fertle et al., 1987).
In the porosity calculation method, using the sonic-neutron crossplot, porosity is determined based on the neutron values versus sonic log. Using this cross- plot both porosity and lithology are measured (Figure 10). Using sonic-density
Figure 9. Neutron-density crossplot.
Figure 10. Sonic-neutron crossplot.
Figure 11. Sonic-density crossplot.
curve, porosity is calculated through drawing the values of density logging versus sonic logging (Figure 11).
5.3. Final Results of Petrophysical Assessment
Final results of petrophysical assessment and loggings of A and B wells are indicated in Figure 12 and Figure 13 with a scale of 1/1000. Based on the results shown in Figure 12, the first right column shows the amount of porosity and the second column shows components of the formation lithology. As is shown in the figure, the dominant lithology is calcareous in the formation (the blue parts) and in some distances there is a little amount of dolomites (the purple parts). These results are correlated with geological reports as well as drilling operations. The average of effective porosity and the total porosity calculated in both A and B wells are shown in Table 2.
5.4. Facies of Kangan Foundation
Kangan formation consists of following three facies:
1) Clean carbonates facies,
2) Clay and shalebaseline, and
Figure 12. The results of petrophysical assessment in Well A (the first column from left shows Calliper, SGR, CGR and IBS loggings; the second column is related to depth and the next columns are related to NPHI, DT, RHOB loggings which are used to calculate the porosity as well as resistance logging).
Figure 13. The results of petrophysical assessment in Well B (the first column from left shows Calliper, SGR, CGR and IBS loggings; the second column is related to depth and the next columns are related to NPHI, DT, RHOB loggings which are used to calculate the porosity as well as resistance logging).
Table 2. The average of effective porosity and total porosity calculated in Wells A and B.
5.5. The Effect of Sedimentary Facies on the Quality of Kangan Formation
Investigation on the effect of sedimentary facies on the quality of Kangan formation indicates that these facies have been created in sabkha environment/tensile zones, ponds, underwater bumps, damns and open seas are related to a carbonate ramp. The most important aggregates of carbonate facies are Ooids, Peloid, Intraclasts, Oncoid and skeletal fragments such as bivalves, gastropods and ostracods. Furthermore, in facies of Kangan formation, different diagenetic processes are observed which are related to atmospheric, marine and burial diagenetic environments. The most important diagenetic processes are dissolution, cementation, micritization by microorganisms, compression, recrystallization, dolomitization and anhydritization. The investigation into core samples from thin sections of carbonate-evaporite facies shows that in Kangan formation anhydrate exists in different thickness and forms which had been formed simultaneously with depositing or in different diagenetic procedures. Anhydrate is widely seen with dolomites. It is also found in the form of independent facies or with various fabrics inside the carbonate facies of this formation. In the Kangan Formation (k4), the anhydrite facies are in laminate, layer and bulk form, with a thickness of several millimeters to several centimeters and sometimes more than one meter. Sulfate minerals are susceptible to various diagenetic processes, especially dissolution, recrystallization and replacement. The diversity of primary and secondary environments for sulfate mineral formation as well as the diversity of diagenetic processes in them has led to the absence of the same classification for the mentioned minerals.
5.6. Anhydrite Textures and Fabrics in Kangan Formation
The main anhydrite textures and fabrics in Kangan formation are:
1) Single nodules: These anhydrate nodules are found in various sizes, ranging from several millimeters to several centimeters. They have changeable forms but they mostly are in spherical and elliptical shapes. Microcrystals to macro crystals are found in them and sometimes they are stretched.
2) Chicken wire netting (poultry netting): The compressed arrangement of anhydrite nodules has led to the emergence of this kind of texture. Anhydrate nodules are surrounded by thin veins of host rock sediment. The host sediment in most samples of this texture is microcrystal dolomite.
3) Enterolithic: In this texture anhydrate nodules are interwoven and have created irregular and folded laminates or layers.
4) Lens: This fabric is a good sign of initial gypsum minerals which are turned to anhydrate through diagenetic processes and usually it is less than one centimeter.
5) Tabular/blade shaped: This fabric is mostly found with anhydrate lens fabric and lens fabric is the result of the conversion of gypsum.
6) Rose: This fabric is like bundles of roses which normally have several stalks with lengths of a few millimeters. This fabric is mostly found with other fabrics such as tabular and lens.
7) Spherulites/fan shape: Aggregation of anhydrate with stretched crystals makes spherulitic/fan shaped form. This fabric is found less than other fabrics.
8) Filler cement: This fabric is one of the most abundant fabrics in facies of Kangan formation especially in facies with greenstone texture. Usually anhydrate cement is spathic and forms between the grains, crystals, inside skeletalgrains, fenestrals, mud cracks, veins and molds due to the dissolution of orogonite grains such as Gastropods, Bivalves and Ooids. In some samples, anhydrate cement is found in the form of poikilotopic cement. This type of cement is especially found in ooid facies in a way that a single anhydrate cement crystal encloses some ooids.
In general, studies show that most textures and fabrics in sedimentary facies are formed simultaneously with deposition processes or at initial stages of diagenesis procedures.
5.7. Determining the Type of Porosity
Using 3 modes of velocity-deviation log (Positive, negative and zero velocity deviations), different types of porosity in k4 Kangan carbonate formation are determined.
5.8. Velocity-Deviation Log
Velocity-deviation log which is created by combining sonic-neutron crossplot and neutron-density crossplot, is an appropriate tool for collecting information about the dominant type of porosity in a well. There is always little difference between calculated porosities using these two techniques. The difference is regarded as secondary porosity which is related to fracture and vuggy porosity.
5.9. The Way of İnterpreting Velocity-Deviation Log
1) Positive velocity deviation (near to 1000)
In this range, the measured sound speed is greater than the calculated velocity. This mode shows porosities whose bed is a strong and cemented rock, and pores which are not linked well and do not have proper permeablity. In this range usually the dominant porosities are reinter-fossil, mouldic and vuggy porosities.
2) ±500 velocity-deviation
±500 velocity-deviation (near to zero) indicates inter-crystal, inter-grain or micro porosity. This range, except for cases where sucrose dolomite makes inter crystal porosity due to diagenesis, shows mild effect of diagenesis. This group has a better permeability than the previous group.
3) Negative velocity-deviation
The range that shows negative velocity-deviation can be interpreted in three ways:
a) Collapsing horizons in well can cause differences in calculated velocity and measured one. In this range, the measured sound speed is less than actual sound speed.
b) Fracture porosity can make negative velocity-deviation. Contrary to the fact that fracture porosity is regarded the same as secondary porosity and in the range of high velocity with positive velocity-deviation. The results of previous studies show that fracture in small and large scale reduces sound velocity so fracture porosity is in the range of negative velocity-deviation.
c) Negative velocity-deviation can be due to free gases in the reservoir as free gases reduce sound velocity.
5.10. Method of Drawing Negative Velocity-Deviation Diagram
First, environmental correction of neutron and density logs is done using Geolog software. Then the volume of shale (Vsh) is calculated through GR diagram and the following equations. The effect of shale on the amount of porosity that is determined by neutron and density logs, is corrected, too  .
According to the above equation, GRmin is the minimum of GR figure which is related to intervals free from shale and GRmax is the value that is read from GR diagram of shale layer. So the volume of shale in both wells is calculated using the following equation  (Qasim al-Askari et al., 2010):
According to above equation, the effect of shale on porosity determined through GR diagram is corrected and is the corrected value of porosity determined through neutron log; is the porosity read from neutron log and is the porosity read from neutron log in shale layer. Thus according to the following equations, the nominal porosity determined through density log is calculated and its shale effect is eliminated  .
fDa is the nominal porosity read from density log, ρma is matrix density and ρfl is the fluid density. In the next step, nominal porosity (fDa) affected by shale is corrected using the following equation  .
In the next step, using the following equation, porosity determined with neutron lg(fN) is combined with density lg(fD) and neutron-density porosity is calculated  :
It should be mentioned that in the above equation the corrected values of neutron and density diagrams are used.
Studying textures and fabrics existing in facies of Kangan formation indicates that most of those textures and fabrics are formed simultaneously with depositing or at the beginning stage of diagenesis. Some nodules are related to dry and hot environment. These nodules are made of inter-hole concentrated fluids through embedded or displacement procedures and in the capillaries of upper part of ferrite zone under Sabkha. These crystals which grow inside deposits are found in matrix of sediments above tensile zone and above inter-tensile zone. Continuing calcium sulfate saturation inside inter-hole fluids, gypsum and anhydrate nodules grow and join together. As a result, chicken wire netting or enterolithic wrinkles are created. Petro physical data of Kangan formation shows that based on the lithology of the studied zone, the dominant sediment is limestone while in some distances it is dolomite. The results indicate that anhydrate nodules mainly expand in upper parts of peritidal parasequences of Kangan facies. Anhydritization may primarily happen due to the increase of salinity of salt water or as a result of an increase in the temperature which occurs when the sea level decreases. This may show that some sedimentary parts of Kangan formation temporarily falls out of water. This phenomenon is very important in stratigraphic studies. In general, the results show that k4 reservoir of Kangan formation is in good conditions regarding saturation, porosity and lithology.
7. Research Suggestions
Due to the importance of prediction and investigation of petro physical properties of carbonate reservoirs and their porosity and permeability, it is suggested that future studies use modern techniques based on statistical probabilities. Multivariate statistical methods (clustering system, variance analysis, artificial neural networks, multidimensional fit, standard correlation analysis, etc.) are among these modern techniques which are based on statistical probabilities and are used in analyzing the characteristics of hydrocarbon wells. It is also suggested that techniques of determining thin zones with high porosity be used in future studies.