Rock typing is an important tool in evaluation and performance prediction of reservoirs. Different techniques such as flow zone indicator (FZI), FZI* and Winland methods are used to categorize reservoir rocks into distinct rock types. Generally, these methods are applied to petrophysical data that are measured at a pressure other than reservoir pressure. Since the pressure changes the pore structure of rock, the effect of overburden pressure on rock typing should be considered. In this study, porosity and permeability of 113 core samples were measured at five different pressures. To investigate the effect of pressure on determination of rock types, FZI, FZI* and Winland methods were applied. Results indicated that although most of the samples remain in the same rock type when pressure changes, some of them show different trends. These are related to the mineralogy and changes in pore system of the samples due to pressure change. Additionally, the number of rock types increases with increasing pressure. Furthermore, the effect of overburden pressure on determination of rock types is more clearly observed in the Winland and FZI* methods. Also, results revealed that a more precise reservoir dynamic simulation can be obtained by considering the reservoir rock typing process at reservoir conditions.
Keywords: Overburden pressure; Carbonate rocks; Rock type; Reservoir quality index; Flow zone indicator; Winland method
Edited by Yan-Hua Sun
Classification of reservoir rocks into different rock types, called reservoir rock typing, is an essential tool in drilling, production and especially reservoir studies. A petrophysical rock type is presented as a group of rock samples that have similar petrophysical and geological properties that influence fluid flow (Stolz and Graves [
Generally, rock typing approaches are performed at a pressure other than reservoir pressure, especially at atmospheric pressure. However, changes in pressure can alter the pore structure of the rock. When pressure is applied to different rocks, they respond differently, consequently a rock sample in a rock type that was determined at the atmospheric pressure may shift to another rock type when pressure changes (see Fig. 1). It should be noted that the effect of overburden pressure on rock is considered in commercial simulators by the rock compressibility (C
Graph: Fig. 1 Schematic of effect of pressure on rock typing process (RT1: Rock type 1, RT2: rock type 2 and RT3: rock type 3). Each circle represents a rock sample
In this work, 113 core samples from one of the carbonate oil reservoirs in the Middle East have been categorized into different rock types using RQI/FZI, Winland and FZI* methods at five different pressures to investigate the effects of pressure on the process of rock type determination. It is worth mentioning that the effect of pressure in the rock typing process has not yet been investigated. In other words, most of the research examines rock type determination at a specific pressure.
In this study, first reservoir rock typing is defined. Then, three selected methods of rock typing are applied to classify studied rock samples at different overburden pressures, and finally, the results of the three methods are discussed thoroughly.
Various techniques have been suggested for classification of reservoir rocks into rock types such as the J-function method, RQI/FZI technique, capillary pressure approach and the Winland method (Soleymanzadeh et al. [
RQI/FZI has been derived from Kozeny–Carman equation which is based on assuming a porous medium as a bundle of capillary tubes. It is obtained by combining Poiseuille's equation and Darcy's law (Zhao et al. [
Graph
where k is permeability (mD), ϕ is porosity, F
Rearrangement of Eq. (
2
Graph
Amaefule and Altunbay ([
3
Graph
Also, RQI is defined as follows:
4
Graph
where k is permeability in mD. Normalized porosity (ϕ
5
Graph
Substituting Eqs. (
6
Graph
Taking logarithms of both sides of Eq. (
7
Graph
where RQI and FZI are in µm, and ϕ
The rock typing methods were frequently used to classify reservoir rocks at atmospheric pressure. It is obvious that values of porosity and permeability in reservoir conditions are different from their values at atmospheric pressure. Therefore, using data at atmospheric pressure may result in an incorrect rock typing process and inaccurate reservoir performance prediction. A proper solution for considering pressure effect on rock type is to perform the RQI/FZI method at reservoir pressure.
Winland performed mercury injection experiments on a large set of sandstone and carbonate rock samples to correlate porosity, permeability and the size of the pore throats. His multiple regression analysis for various mercury saturations revealed that the best correlation coefficient (R
8
Graph
where r35 is in µm, k is uncorrected air permeability in mD, ϕ is porosity in percentage.
r35 can be used as a basis to classify a reservoir into different rock types. All rock samples with similar r35 constitute a single rock type and lie on an iso-pore throat curve.
The base form of Kozeny–Carman equation is obtained by combining Poiseuille's equation and Darcy's law, as noted in Sect. 2.1. This form of Kozeny-Carman equation is as follows:
9
Graph
where r
10
Graph
Herein, FZI* is in µm which can be calculated for each sample from measurement of its porosity and permeability. Hence, rocks with the same FZI* lie within an individual group. The fluid flow behavior of this group is assumed to be the same. Taking logarithm from both sides of Eq. (
11
Graph
It is inferred from Eq. (
In this work, permeability–porosity data related to 113 carbonate core samples from a carbonate reservoir were used. These data have been obtained at different confining pressures (atmospheric pressure, 2000, 4000, 5000 and 6000 psia). Porosity and permeability of these rock samples at atmospheric conditions are illustrated in Fig. 2.
Graph: Fig. 2 Porosity and permeability of 113 core samples at atmospheric conditions
The value of r35 at atmospheric conditions is calculated for all samples from Eq. (
Graph: Fig. 3 r35 of all samples at atmospheric pressure
Table 1 summarizes the average (Ave) and median (Med) of permeability, porosity and r35, FZI and FZI
Average and median of porosity, permeability and r35 of the rock samples at different pressures
Pressure Permeability, mD FZI, µm FZI*, µm Ave Med Ave Med Ave Med Ave Med Ave Med 14.7 10.16 0.98 0.169 0.153 0.984 0.534 0.612 0.465 0.138 0.083 2000 9.54 0.7 0.153 0.139 0.947 0.538 0.626 0.460 0.131 0.080 4000 9.26 3.51 0.15 0.135 0.920 0.505 0.611 0.461 0.127 0.073 5000 9.05 0.41 0.149 0.135 0.902 0.455 0.60 0.456 0.125 0.071 6000 8.98 0.45 0.148 0.131 0.845 0.444 0.595 0.450 0.0123 0.070
The classical approach to reservoir rock typing, a semi-log plot of permeability versus porosity (Abbaszadeh et al. [
Graph: Fig. 4 Classical method of rock typing: log K versus ϕ
The first step of the rock typing process is data clustering. There are different clustering techniques can be used in rock typing processes, such as discrete rock type (DRT), histogram, parabolic plots and global hydraulic element (Abbaszadeh et al. [
In order to investigate the effect of pressure on rock type determination, the RQI/FZI method was applied to determine rock types at different pressures. These rock types are illustrated in Fig. 5. Comparing rock types at different pressures reveals that the rock type of core samples changes in various ways:
Graph: Fig. 5 RQI/FZI methods at different overburden pressures
- Increasing trend (shift from lower rock type to upper one): such as core No. 51 which has been denoted by symbol in Fig. 5.
- Decreasing trend (shift from upper rock type to lower one): for example, core No. 47 which has been shown by symbol in Fig. 5.
- Fluctuating trend: such as core No. 19 which has been indicated by symbol in Fig. 5.
- No change: major part of studied samples remained in the same RQI/FZI rock type.
Table 2 presents the number of samples for each mentioned trend.
Four different trends due to pressure change based on the RQI/FZI method
Trend Percentage, % Remarks No change 60 58% dolomite, 52% vuggy and 28% anhydrite Increasing 10 Decreasing 23 76% limestone, 88% vuggy Fluctuating 7
In order to clarify the abovementioned trends, for each trend, some samples were selected and their FZI values were plotted versus pressure in Fig. 6. In fact, each trend in Table 2 was named according to the change in FZI value versus pressure as it is shown in Fig. 6.
Graph: Fig. 6 Change in FZI during increasing pressure for different trends
Figure 7 shows the number of rock samples in each rock type. This figure reveals that the number of samples in the rock types with low values of FZI (EX1, RT7, RT8 and RT9) increased by increasing pressure. It should be emphasized that rock types EX1 and EX2 did not exist at atmospheric pressure and were added to the other rock types when pressure was increased. It means that by increasing pressure the number of rock types increases.
Graph: Fig. 7 Number of samples in each rock types based on the RQI/FZI method
Since, permeability mostly depends on pore throat size rather than pore size, the authors believe that using the Winland method which contains pore throat size (r35) leads to a clearer description of the effect of pressure changes on the rock type determination. Whereas the RQI/FZI method is based on the Kozeny–Carman model in which the pore radius and pore throat are considered equal. In order to investigate the effect of pressure on the rock type determination by the Winland method, this method was applied to rock samples at five different pressures: 14.7, 2000, 4000, 5000 and 6000 psia (see Figs. 8, 9, 10, 11, 12).
Graph: Fig. 8 Flow units based on the Winland method at 14.7 psia
Graph: Fig. 9 Rock types based on the Winland method at 2000 psia
Graph: Fig. 10 Rock types based on the Winland method at 4000 psia
Graph: Fig. 11 Rock types based on the Winland method at 5000 psia
Graph: Fig. 12 Rock types based on the Winland method at 6000 psia
Figures 8, 9, 10, 11 and 12 show that most of the rock samples shift to the left and downward simultaneously. In other words, this leads to change in the number of rock types and also changing a rock sample from one rock type to another one. In addition, these figures indicate that the number of data in the low k–ϕ zone (blue circle) increases with an increase in pressure.
Figure 13 depicts the number of rock samples in each of the rock types at different pressures. Three points are inferred from this figure: (
Graph: Fig. 13 Frequency of rock types based on the Winland method
The shift of the rock samples between different rock types (based on the Winland method) during pressure changes was examined, and results are reflected in Table 3. Indeed, Table 3 reveals that 37% of rock samples jump from one rock type to another one due to change in pressure. This means that ignoring the effect of pressure on the determination of rock types and considering k-ϕ at atmospheric pressure, make large errors in subsequent processes in a reservoir study.
Four different trends due to pressure change based on the Winland method
Trend Percentage, % Remarks No change 63 60% dolomite, 54% vuggy and 30% anhydrite Increasing 2 Decreasing 31 82% limestone, 81% vuggy Fluctuating 4 80% anhydrite
Further investigations imply that 60% of rock samples which had remained in the same rock type during changes in pressure are dolomitic. This may be due to lower compressibility of dolomite rock samples with respect to limestone samples. Also, 82% of rock samples which shift from upper curves to lower curves are limestones. It should be noted that most of these samples contain vugs. It seems that the high compressibility of these vuggy limestone samples is the main reason of this trend of Table 3. A few samples (2%) jump from lower curves to upper curves which may be related to generation of fractures in the pore structure of these samples due to an increase in pressure. The fluctuating trend in Table 3 can be attributed to the generation of induced fractures and closeness of some pores in successive steps of pressure changes. It is worth mentioning that 80% of samples with a fluctuating trend contain anhydrite. Further investigation is required to explain the effect of anhydrite content on the fluctuating trend of a rock sample. Figures 14, 15, 16 and 17 illustrate four trends of Table 3: no change, decreasing, increasing and fluctuating, respectively. In these four figures, arrow direction indicates the path of change of rock types during pressure changes.
Graph: Fig. 14 Samples remain in the same rock type during pressure changes (no change trend)
Graph: Fig. 15 Samples jump from upper curves to lower curves during pressure changes (decreasing trend)
Graph: Fig. 16 Samples shift from lower curves to upper curves during pressure changes (increasing trend)
Graph: Fig. 17 Samples fluctuate between different rock types during pressure changes (fluctuating trend)
The value of r35 at different pressures was used to explain the observed trends in Fig. 18. This figure shows the value of r35 at different pressures for four different trends. (Each part of the figure is related to one trend in Table 3.) It is inferred from this figure that, all trends in Table 3 can be interpreted based on the change in r35 during pressure changes.
Graph: Fig. 18 r35 changes due to an increase in pressure for four different trends of Table 3
Using the FZI* method, the number of rock types increased from six to eight with an increase in overburden pressure (see Fig. 19). Rock samples move to the left and downward simultaneously, which obviously implies that the quality of the majority of the rocks reduces with an increase in pressure. Comparison of rock types at five different pressures shows that rock types change in two trends: decreasing and fluctuating. Core 19 (symbol ) and Core 65 (symbol ) represent decreasing and fluctuating trends, respectively.
Graph: Fig. 19 Rock typing at different pressures using the FZI* method
The frequencies of FZI* are demonstrated in Fig. 20 which confirms the results obtained from the other two methods. Increasing pressure causes rock types EX1 and EX2 to be added to the rock types at atmospheric pressure.
Graph: Fig. 20 Frequency of rock type based on the FZI* method
Table 4 presents the effect of pressure on the rock typing process by the FZI* method and details of observed trends due to pressure change. This table shows that, similar to RQI/FZI and Winland methods, more than 50% of studied rock samples have remained in their rock types at atmospheric pressure. Furthermore, a decreasing trend is the most common trend and vuggy limestone samples are majority of the rocks which fall within this trend, as observed in RQI/FZI and Winland methods.
Four different trends due to pressure change based on FZI* method
Trend Percentage, % Remarks No change 55 61% dolomite, 52% vuggy and 28% anhydrite Increasing 0 Decreasing 42 72% limestone, 79% vuggy Fluctuating 3 67% vuggy and 67% anhydrite
Finally, it is noted that having a clearer picture of the rock pore structure, such as from micro-computed tomography (Micro-CT) scans, improves the analysis of the effect of pressure on the determination of rock types.
The following conclusions arose from this work:
- Studied samples were classified into different rock types using RQI/FZI, FZI* and Winland methods at five different pressures. Different behavior was observed for rock samples during changes in pressure. The majority of the samples remained in the same rock type during pressure increases. Some of the samples shifted from an upper curve to a lower curve, and a few samples change from a lower curve to an upper one. In addition, several of the rock samples showed fluctuating trends. These four different trends can be attributed to the mineralogy and change in pore structure of the studied samples.
- Most of the rock samples which remained in the same rock type during pressure changes are dolomitic. It seems that this is related to the lower compressibility or higher density of this type of rock. In contrary, the drastic changes in rock types occur in the limestone rock samples which contain vuggy porosity. The higher compressibility of these samples is the main reason of this behavior.
- In RQI/FZI, FZI* and Winland methods, it is observed that the number of rock types increases with an increase in pressure. Also, the number of rock samples in the lower rock types (the lower quality rock types) increases. Furthermore, generally, at pressures greater than a specific pressure (in this study, 4000 psia), the number of rock samples in the higher rock types (the higher quality rock types) remains constant.
- The Winland method gives a clearer picture of changing rock samples between different rock types. This is because the Winland method has been developed based on the size of pore throats (r35).
- The effect of pressure on the rock type determination implies that the process of reservoir rock typing should be performed at the reservoir conditions.
By Aboozar Soleymanzadeh; Saeed Parvin and Shahin Kord
Reported by Author; Author; Author