An Integrated Rock Typing Approach for Unraveling the Reservoir Heterogeneity of Tight Sands in the Whicher Range Field of Perth Basin , Western Australia

Tight gas sands in Whicher Range Field of Perth Basin show large heterogeneity in reservoir characteristics and production behavior related to depositional and diagenetic features. Diagenetic events (compaction and cementation) have severely affected the pore system. In order to investigate the petrophysical characteristics, reservoir sandstone facies were correlated with core porosity and permeability and their equivalent well log responses to describe hydraulic flow units and electrofacies, respectively. Thus, very tight, tight, and sub-tight sands were differentiated. To reveal the relationship between pore system properties and depositional and diagenetic characteristics in each sand type, reservoir rock types were extracted. The identified reservoir rock types are in fact a reflection of internal reservoir heterogeneity related to pore system properties. All reservoir rock types are characterized by a compacted fabric and cemented framework. But distribution and dominance of diagenetic products in each of them depend on primary depositional composition and texture. The results show that reservoir rock typing based on three aspects of reservoir sandstones (depositional properties, diagenetic features and petrophysical characteristics) is a suitable technique for depiction of reservoir heterogeneity, recognition of reservoir units and identifying factors controlling reservoir quality of tight sandstones. This methodology can be used for the other tight reservoirs. Corresponding author. R. K. Ilkhchi et al. 374


Introduction
Tight gas sands, as a type of unconventional reservoirs, have been studied and recognized in several countries such as United States of America, Canada, Australia and China.The importance of these reservoirs has been increasingly recognized in recent years because these reservoirs, in spite of their low porosity and especially low permeability, have large reserves of hydrocarbon that is obtainable if they are characterized exactly and extracted by new techniques.Generally, tight gas sand is defined as a low permeability reservoir that need specific treatments such as advanced drilling technology, correct way of completion and stimulation jobs for having a commercial gas production [1].
There are many studies on tight gas focusing on different aspects of these reservoirs both geologically and petrophysically.Stroker et al. [2] studied diagenesis of tight gas sand reservoir of Upper Cretaceous Mesaverde Group in the Picean Basin of Colorado, and they concluded porosity to be a complicated function of compaction, cementation and dissolution.Prodanovic et al. [3] in their study of diagenetic alteration by numerical simulation in tight gas sandstones showed that cementation has the main effect on gas trapping in these reservoirs.Rezaee et al. [1] estimated permeability of the tight gas sands in Western Australia from mercury injection capillary pressure and nuclear magnetic resonance data.They showed that pore throat radius at 10% mercury saturation (r 10 ) yields the best correlation with permeability.Bahrami et al. [4] studied hydraulically fractured tight sand gas reservoirs of the Perth Basin in Western Australia and introduced the importance of formation damage by water blocking.Ziarani and Aguilera [5] estimated pore-throat radius and tortuosity from formation resistivity data for tight-gas sandstone reservoirs.In recent years, there are also many studies on tight gas reservoirs in China that have considered different aspects of their reservoir characterization (e.g.[6]- [10]).
Reservoir sands of the Whicher Range gas field in the Perth Basin of Western Australia are mainly considered as tight gas sands due to their low porosity and low permeability.Production from these sands is not stable and cost-effective.On the other hand, uses of certain techniques such as hydraulic fracturing to stimulate the reservoir and improve the production, have not been successful until now.It is, therefore, the target of this study to characterize these sands related to their internal heterogeneity and to determine main factors controlling reservoir behavior related to pore system properties.In this respect, we have applied reservoir rock typing based on the integration of the depositional properties and diagenetic features associated with petrophysical characteristics as a comprehensive method to characterize reservoir sandstones.The results show the effect of both depositional and diagenetic properties on pore system of tight sands in the Whicher Range gas field.

Geological Setting
Perth Basin is an N-S trending, long and narrow and asymmetrical rift basin (Figure 1(a)) located in the southwest of Western Australia [11].The eastern boundary is the Darling Fault, and the basin extends offshore to the continental-oceanic boundary.The basin contains up to 12 km of Ordovician and mid-Carboniferous to Cretaceous, primarily non-marine to shallow-marine strata, below extensive thin, mostly non-marine, Cenozoic cover [12].Perth Basin hosts some tight sand reservoirs situated in both southern and northern sectors.The Whicher Range field, the target of this study, with gas trapped in deep burial (~4 km) tight sands is one of the most important fields of the basin.The field is located in a large anticline structure in the Bunbury Trough, and five wells have been drilled in this structure so far.The sandstones of Willespie Formation of Late Permian age (Figure 1(b)) constitute the main reservoir rock of the field [13].

Data and Methodology
In order to characterize tight sands of Whicher Range field, we used well log and core data from four drilled wells in the field.Facies and diagenetic features were recognized and interpreted based on core, cutting and microscopic (polarizing, SEM) studies.In order to study and classify reservoir sandstones in relation to their reservoir  to Triassic Lithostratigraphy column of the southern Perth Basin [13].In this column, Late Permian Willespie Formation from Sue Group, as the main reservoir rock of the Whicher Range field, has been highlighted.
characteristics, pterophysical rock typing was used in terms of hydraulic flow units (HFU) and electrofacies concepts.Core porosity and permeability data were used for recognition and differentiation of HFUs.Reservoir electrofacies were determined based on well log data clustering using the cluster tree method by using Matlab software.In an effort to find a reasonable and clear relationship between recognized petrophysical rock types with their depositional and diagenetic characteristics, we described reservoir rock types.Each rock type is characterized by three main parameters including depositional facies (texture and composition), diagenetic facies (diagenetic features), and petrophysical characteristics (pore system properties).

Detrital Composition
Willespie Formation is mainly composed of poorly sorted feldspathic sandstones (feldspathic arenite to subfeldspathic arenite according to Folk et al. [14], Figure 2), with subordinate conglomerate, siltstone, shale and sporadic lenticular coal seams [13].Sandstones are generally texturally and compositionally immature to submature.Texturally, sandstones show a wide range in size from very fine to very coarse and granular.Framework grains are mostly sub-angular to sub-rounded.Sorting is generally poor, but is better in fine grained facies.Monocrystalline quartz is the dominant grain type.Feldspar (potassium feldspar and plagioclase) is usually the second most abundant framework grain.Less abundant framework grains are rock fragments (metamorphic, igneous and sedimentary lithics), micas (biotite and muscovite), and heavy minerals (garnet, zircon, tourmaline, and epidot).In some facies, the abundance of micas and garnet significantly increases.Based on the regional geology, the existing core interpretation and vertical facies variations, a fluvial-paludal system is considered for the Willespie Formation.

Diagenetic Overprints
Diagenetic events including compaction, cementation and dissolution during early to burial stages of diagenesis have severely affected the reservoir sandstones of the Whicher Range field.The effects of these processes are clearly seen in modification of the primary pore system properties which, in turn, control the reservoir quality and production behavior of the field.Photomicrographs of common diagenetic features within the reservoir sandstones have been shown in Figure 2.

Pore System in Whicher Range Sandstones
The pore system in Whicher Range sandstones contains both primary and secondary pores.Both pore types have been cemented by clays, quartz and carbonates.The secondary pores are connected through the residue of the strongly cemented primary pore network, and its permeability will fundamentally determine the reservoir quality of the sandstones.In addition, micro-porosity related to authigenic clays such as kaolinite and illite/smectite with large surface area supports high irreducible water saturation.All of these features define the reservoir sandstones as a type of unconventional reservoirs named as tight sand or tight gas sands.In Figure 3, some SEM photographs of pores system of the reservoir sandstones affected by authigenic clay minerals have been demonstrated.

Relationship between Diagenetic Products and Tightness in Reservoir Sanstones
In tight gas sandstones of Whicher Range, as already discussed, diagenetic processes have severely modified pore system quantitatively and qualitatively and finally resulted in sandstones with compacted fabric and low porosity and low interconnected pore spaces.The results from core studies, petrography interpretation and SEM photographs show that diagenesis has close relationship with primary depositional characteristics, both texture and composition, of Whicher Range sands.In such deep sandstones, compaction has a significant effect on pore systems in the first stages of diagenesis, especially in fine grained facies and in sandstones with a large proportion of ductile framework grains.This can be evidenced by physical deformation and bending of soft grains such as micas (Figure 2(a)) and some rock fragments (Figure 2(d)) around harder framework grains.Precipitation of quartz overgrowths, carbonate cements and authigenic clay minerals mostly occurred after significant physical compaction of the reservoir sandstones.Quartz overgrowths generally have sparse and limited distribution in the reservoir although they are abundant in some quartz rich facies.These cements apparently have more important effect on porosity than permeability.Authigenic clay minerals have extensive distribution and as they are occluding the pore spaces and pore throats they are the main diagenetic cements in relation to reducing the permeability within the reservoir sandstones.Based on microscopic evidence, authigenic clays mostly derived from the breakdown of feldspar and detrital mica, and they infill the interstices between the framework grains.So, distribution of authigenic clay minerals as grain coating and pore lining (i.e.chlorite, mixed layer illite-smectite) as well as pore filling (mainly kaolinite) has intimate relationship with primary mineralogical composition of the sandstones.Calcite cementation has had an important role on porosity and permeability reduction, but it has a local distribution in the reservoir.Dissolution of unstable grains such as plagioclases is considered as a positive process on reservoir properties and to some extent compensates the effect of compaction and cementation, but the resultant secondary pores are isolated and are often partly to completely filled with authigenic cements, especially clay minerals.So, it can be concluded that tightness in Whicher Range sandstones is a complex function of interaction between primary depositional properties and post-depositional diagenetic processes.In this respect, compaction and cementation are two main factors controlling pore system characteristics of the Whicher Range sandstones.As a result, tight reservoir sandstones with compacted fabric, highly cemented and low interconnected pore systems have been produced.

Petrophysical Characteristics of Reservoir Sandstones
The results from core studies and microscopic analysis of Willespie Formation sandstones show an internal heterogeneity related to depositional properties and diagenetic features within the reservoir.This, in turn, is reflected in petrophysical characteristics of the sandstones.In this regard, reservoir behavior can be interpreted based on pore system properties controlled by depositional and post-depositional characteristics of the reservoir sandstones.Based on sedimentary rock typing, reservoir units (sand and silt) of Whicher Range sandstones were distinguished from non reservoir units (shale) using clustering of well log data and verification of results with petrographic data [16].In this study, in order to characterize the reservoir sandstones and to recognize the main factors controlling reservoir quality, reservoir units (sandy and silty intervals) first studied petrophysically in terms of hydraulic flow unit and electrofacies concepts resulted in separation of petrophysical rock types.Then, these rock types based on pore system properties associated with depositional and diagenetic attributes were interpreted, resulted in differentiation of reservoir rock types.

Petrophysical Rock Typing Based on Hydraulic Flow Units
Rock typing based on the hydraulic flow unit (HFU) concept is a technique by which the reservoir rocks are petrophysically classified into certain reservoir units related to pore properties within each unit.In the literature, a hydraulic flow unit is considered as a part of reservoir rock within which geological and petrophysical properties, controlling the fluid flow, are consistent and predictably different from properties of other parts of reservoir rocks ( [17] [18]).Amaefule et al. [19] defined Hydraulic Flow Unit as the elements of the reservoir rock which have similar averages of petrophysical rock properties that affect fluid flow.In this study, we use the method of Amaefule et al. [19] for differentiation of hydraulic flow units in the reservoir that has been recognized and used as a practical and common method by many researchers (e.g., [20]- [23]).The stages are summarized as follows: 1) Calculation of reservoir quality index based on porosity and permeability data and according to Amaefule et al. [19] method: where RQI is reservoir quality index, µm; k is permeability, mD; and ϕ is porosity in fraction.
2) Calculation of pore to matrix ratio (PMR) using the following equation: 3) Extraction of flow zone indicator (FZI) by the following correlation: Some properties of different HFUs associated with cut off values for their differentiation are given in Table 1.Study of identified hydraulic flow units on porosity and permeability plot (Figure 4) shows that each HFU is characterized by a certain district on the plot.This property in fact is related to pore system properties controlling fluid flow in the reservoir.Thus, each hydraulic flow unit can be indicator of a specific reservoir unit with unique pore system properties.Generally, hydraulic flow units in relation to their petrophysical characteristics (Table 1) can be classified into three main groups as follows: HFU A&B: low permeable hydraulic flow units A and B are related to very fine grained sands and silts.In these facies, although porosity may be relatively high but their permeability is very low (<0.1 mD).So, they cannot have important role on production behavior in the reservoir.In this respect they are similar to shaly intervals, and we discarded them for more investigation.Two other groups, mainly separated based on porosity value, are HFU C&D and HFU E&F with average porosity of 10% and 6%, respectively (Table 1).Based on core and petrographic evidences, sandstone facies in these two groups are very heterogeneous, and change from fine to coarse and very coarse.So, we used data clustering technique based on prtrophysical log responses (DT, NPHI, and RHOB) and pore size to extract petrophysical rock types, and unravel reservoir characteristics in each group individually.Pore size was determined by Equation ( 4) which was developed by Rezaee et al. [1] for tight gas sands.
where k is the dry gas permeability (mD), Phi is porosity (%) and r 10 is the pore throat size (micron) corresponding to the 10th percentile of mercury saturation on a cumulative mercury injection plot.

Petrophysical Rock Typing Based on Electrofacies
In this stage, we use well log data responses based on electrofacies concept to investigate and classify different petrophysical rock types within the reservoir sandstones.Electrofacies approaches are often used to determine rock types from logs.In this regard, electrofacies is applied for a set of log responses which characterizes bed/strata and allows it to be distinguished from others [24].This response against reservoir rocks is in fact a reflection of variability related to geological and petrophysical properties of the reservoir recorded continuously along the well.Thus, separation of reservoir facies based on available logs known as log facies or electrofacies, can be used as a useful, practical and economical procedure in this regard.Electrofacies in each group of defined hydraulic flow units (HFU C&D and HFU E&F), were determined individually, by hierarchical cluster analysis of well log data.The processed log data (i.e., GR, DT, RHOB, and NPHI) were loaded as an input matrix (x) in MATLAB software workspace, and we pursue the methodology according to the following stages.
Step 1: Compute the distance between data objects.We compute the Euclidean distance (Equation ( 5)) between pairs of objects as a common method.

(
)( ) where d for two clusters r and s is the various distances between the vectors x r and x s for an m by n data matrix x, which is treated as m (1-by-n) row vectors 1 2 , , , m x x x  .
Step 2: Linking between distance data.In this step, an appropriate linkage function create hierarchical cluster tree from the distance information generated in the first step.We use complete linkage function (Equation ( 6)) that uses the largest distance between objects in the two clusters, and is more compatible method with the reservoir sandstones: clusters or electrofacies has certain values of readings or responses of any used well logs that in fact reflects geological (textural and diagenetic) and petrophysical (porosity and permeability) characteristics of the reservoir rocks.
Study of extracted electrofacies related to petrophyical characteristics shows that each electrofacies has specific districts on porosity and permeability plot (Figure 6).This is mostly related to pore system properties.Thus, they can be classified into three main petrophysical groups or sand types including sub-tight (K = 1 -10 mD; EF3, EF5, and a part of EF4), tight (K = 0.1 -1 mD; EF1, EF6), very tight (K < 0.1 mD; a part of EF2) sandstones.Facies and petrophysical characteristics of each sand type are as follows.
Very tight to tight sandstones: these sandstones that are mostly related to EF2, consist of fine to medium grained sandstones, in most cases they are micaceous with numerous shale partings.The average porosity and permeability in these sandstones are 7% and 0.19 mD, respectively.
Tight sandstones: these sandstones are related to EF1 and EF6.They consists of medium to coarse and to some extent fine grained sandstones.In some facies mica is dominant with numerous shale partings.Average porosity and permeability for EF1 tight sandstones are 11% and 0.46 mD respectively and for EF6 tight sandstones they are 4% and 0.65 mD.
Sub-tight sandstones: sandstone facies in this sand type are compatible with EF3 and EF5 and a part of EF4, and consist of medium to coarse and in some case very coarse and granular sands.A minor part of porosity in these sands has remained as intergranular.Facies in EF4, petrophysically change between sub-tight and tight sandstones.Average porosity and permeability for these sand types are respectively: 19% and 2.83 mD in EF3, 10% and 2.83 mD in EF5, and 6% and 1.36 mD in EF4.

Reservoir Rock Typing of the Willespie Formation Sandstones
Study of tight sands in the Whicher Range field based on hydraulic flow units and reservoir electrofacies, as a type of rock typing, shows some aspects of internal resevoir heterogenity.But other aspects of reservoir characteristics with respect to pore system properties responsible for reservoir bahavior remains unknown.So, in an effort to recognize factors affecting pore systems and resulting reservoir properties, we consider the role of depositional and diagenetic processes on pore systems of sand types in the reservoir.In this respect, in the framwork of what we name reservoir rock typing, we study and analyse three main aspects of reservoir sandstones including depositional properties (texture and composition), diagenetic features (diagenetic facies) and petrophysical characteristics (pore system properties) with each other to show a real picture of reservoir units.Depositional properties that are related to texture and composition of reservoir sandstones, as already discussed, have some control on the pore system and also the distribution of diagenetic products.The appearence of the sandstone  facies has been changed as a result of diagenetic activities within the reservoir.This can be observed in compacted fabric and cemented framework of all sand types, and resulted in quantitative and qualitative modification of pore system in these sandstone, as they are considered as tight gas sandstones.Pore system properties of tight sandstones in the Whicher Range field is responsible for production behavior of the reservoir.Thus, knowledge of pore space properties (type, size, volume and pore connection) of the reservoir sandstones, can help in recognition of production zones and interpretation of reservoir behavior in the field.Table 2 summarizes pore system properties and descriptions in relation to different sand types and electrofacies in the reservoir.

Reservoir Rock Types
Based on the results from the integration of three aspects of the reservoir sandstones (depositional properties, diagenetic features and petrophysical characteristics) reservoir rock types are described.Investigation of different sand types in terms of depositional, diagenetic and petrophysical characteristics resulted in discrimination of eigh reservoir rock types (RRTs).Quantitative and dominant qualitative features of these RRTs are shown in Table 3 and Table 4, respectively.

Interpretation of Reservoir Rock Types in Each Sand Type
The main findings of reservoir rock typing of tights sandstones are summarized as follows.
• Compaction is a common diagenetic feature in all reservoir rock types.But the effect of these events is more striking in fine grained and mica rich sandstones and also in facies with significant amount of ductile grains and shale partings, especially in very tight to tight sandstones (RRT7).• Cementation by authigenic clays, Q-overgrowths and calcite has an important role on porosity reduction and limitation of pore connections.In this respect, pore-filling and pore-lining clay cements, mostly derived from alteration of feldspars and micas, have widely been distributed in both primary and secondary pores.• Quartz overgrowths with sparse distribution within the reservoir sandstones, have mostly affected the pore system through reduction in pore volume.But, the effect of this cement on permeability is not significant.• Local distribution of calcite cement within the reservoir has been recognized as pore occluding and also grain replacive, especially in coarse grained facies.Facies with significant calcite cement in their pore system, have very low porosity and permeability.These facies involve part of very tight to tight sandstones in the reservoir (RRT-8).• Pore system includes isolated dissolution pores and also micro pores associated with authigenic clays.However their connection is through small remaining intergranular pores which especially are preserved in coarse grained sandstones.• Study of reservoir rock types shows that all sandstone types are characterized by a compacted fabric and cemented framework, but medium to coarse/granular sandstones, mostly compatible with channel deposits, with medium to high porosity and moderately low permeability (RRT1 and RRT2) are the most promising intervals in the reservoir.

Conclusions
In this study, tight gas sandstones of the Whicher Range field were studied with respect to depositional, diagenetic and petrophysical characteristics.The study shows that tight sandstones have severely been affected by diagenetic events.Such events depend on primary depositional properties of reservoir sandstones, and have modified the pore system responsible for reservoir behavior in the field.We recognized reservoir rock typing based on three depositional, diagenetic and petrophysical aspects of the reservoir sandstones, as a comprehensive and practical method to reveal internal reservoir heterogenity and to recognize main factors controlling pore system propeties.Consequently, eight reservoir rock types were recognized.All theses reservoir rock types are characterized by compacted fabric and cemented framework.But the pore system in each reservoir rock type, depends on depositional properties and distribution of diagenetic products, resulting in different properties, which differentiate it from other reservoir units.The results show that medium to coarse/granular sandstones, mostly related to channel deposits, with medium to high porosity and moderately low permeability are the most promising intervals from reservoir quality point of view.

Figure 1 .
Figure 1.(a) Location map of Perth Basin in the southwest of Western Australia.The location of Whicher Range field in Bunbury Trough has been shown in the southern Perth Basin [15].(b) Permian to Triassic Lithostratigraphy column of the southern Perth Basin[13].In this column, Late Permian Willespie Formation from Sue Group, as the main reservoir rock of the Whicher Range field, has been highlighted.

Figure 3 .
Figure 3. SEM photographs of clay minerals in the Whicher Range sandstones.(a) General view of sample showing detrital grains (G) surrounded by clay (C) (WR2-3942.5m).(b) The pore system is coated with a well-developed film of smectite (S).Note the cellular morphology which gives the coat a large surface area (WR2-4171.1 m).(c) The primary pores are often filled with kaolinite booklets (K) and are coated with a film of iron oxide (F) (WR2-3945.51m).

Figure 4 .
Figure 4. Distribution of different hydraulic flow units on the porosity and permeability plot for the reservoir sandstones.

Step 3 : 5 . Step 4 :Figure 5 .
Figure 5. Generated dendrogram diagrams for cluster analysis of reservoir sandstones in each group of hydraulic flow units: (a) (HFU C&D) and (b) (HFU E&F).The horizontal dashed line on each diagram shows the cutoff value which based on different clusters is differentiated.

Figure 6 .
Figure 6.Distribution of different electrofacies related to HFU groups (HFU C&D and HFU E&F) on porosity and permeability plot.Electrofacies in relation to various sand types have been shown with different colors: very tight to tight (K = 0.19 mD in average) in yellow, tight (K = 0.1 -1 mD) in brown and sub-tight (K = 1 -10 mD) in red color.

Table 1 .
Cut off values based on FZI for differentiation of HFUs in the reservoir sandstones, associated with reservoir quality index and average porosity and permeability data.

Table 2 .
Pore system properties in relation to identified sand types and electrofacies in the reservoir.

Table 3 .
Quantitative parameters for each reservoir rock type of Whicher Range sandstones.

Table 4 .
Descriptive characteristics of reservoir rock types defined based on three aspects of reservoir sandstones.