1. Introduction
Permeability is a crucial factor in the characterization of oil and gas reservoirs as it regulates how fluids flow through it [
1]. The complexity of formation continuity and, particularly, the differences in pore space-related parameters like porosity, permeability, and capillary pressure are becoming increasingly apparent to the oil industry. Among all the formation parameters that oil and gas engineers use, permeability is one of the most crucial. It is used to decide whether to complete and put a well into operation or to leave it abandoned. The initial depositional process as well as later diagenetic and tectonic modifications are reflected in these variances. For the purpose of forecasting the performance of reservoirs and creating a management plan for field production, simple models frequently fall short. Reservoir engineers are realizing more and more how critical it is to optimize recovery based on the reservoir properties. As a result, a precise understanding of the distribution of permeability both vertically and laterally is crucial.
Because it regulates production rate, permeability is a crucial factor in reservoir development and management. In general, increased porosity, larger grains, and better sorting all result in an increase in permeability. Permeability has been studied extensively in carbonate formations for many years. It has long been thought that vertical permeability in carbonate formations is problematic. In carbonates, pore connectivity is the primary factor influencing permeability [
2]. Carbonate reservoirs exhibit heterogeneity as a result of different depositional conditions and subsequent diagenetic processes. Most heterogeneous reservoirs are composed of formations that have varying pore geometry and/or contain shale. Reservoir heterogeneity is now seen as an opportunity to recover more oil rather than a concern in this instance. The reservoirs are considered to be nonuniform, isotropic, and homogeneous in early reservoir engineering research. But in more recent research, a more accurate representation of the reservoir’s complicated geological consistency has been obtained through the application of permeability–porosity relationships.
The deposited environment has a crucial role in the formation of primary porosity [
3]. Porosity and permeability are typically higher in high-energy deposits than low-energy deposits, which have lower permeable intervals. Low-energy deposits frequently have high porosity; however, permeability may also be low if the pore throat size is too tiny. Porosity is both created by and destroyed by diagenesis [
4]. For instance, early cementation may prevent compaction and retain primary porosity by reducing porosity and permeability, whereas dissolution primarily enhances porosity with the exception of specific stylolites that may form barriers. Because it is difficult to assess directly, permeability prediction presents a challenge for information appraisal and reservoir modeling. Understanding permeability is crucial for creating 3D reservoir models, comprehending oil and gas production, and, ultimately, developing a development strategy [
5].
Traditional methods of determining permeability include core analysis and well testing. These traditional techniques are very time and money consuming. The difficulty of precisely estimating permeability in diverse carbonate reservoirs is challenging [
6]. Several studies attempted to use intricate mathematical formulae to relate permeability and reservoir characteristics; nevertheless, this led to an inaccurate calculation of the formation permeability values. The most popular technique for figuring out permeability is called conventional core analysis, which forms a nonlinear relationship between porosity and permeability [
7]. By using this connection, permeability as a function of porosity can be predicted. Vertical permeability is one of the most important reservoir characteristics due to its great role in various oil field activities such as optimizing the well location, production rate, directional drilling, completion method, and perforation properties. Permeability anisotropy, the ratio of vertical to horizontal permeability, is vital in reservoir modeling due to its role in both vertical and horizontal wells.
In multilayered reservoirs, each layer has distinct vertical permeability [
8]. Thus, these are reservoirs usually classified based on their permeability anisotropy (the ratio of vertical permeability (Kv) to horizontal permeability (Kh)). In earlier reservoir engineering research, reservoirs were proposed to be nonuniform but isotropic and homogeneous. However, recent studies describe the reservoir as heterogeneous in terms of permeability and porosity [
9,
10]. Generally, the horizontal permeability is higher than the vertical permeability for rocks with nonuniform and large grains [
11,
12,
13,
14], which is the case in most hydrocarbon reservoirs. Clark [
14] also demonstrated that the permeability of the rock would be very high and nearly identical in both vertical and horizontal directions if it were primarily made up of large, uniformly rounded grains. Generally speaking, horizontal permeability is greater than vertical permeability, especially when the sand grains are tiny and irregularly shaped. This group includes a large proportion of the oil and gas reservoirs.
Various studies have been performed to estimate the absolute permeability of different carbonate and sandstone reservoirs [
15,
16,
17,
18,
19]. It is important to mention that predicting permeability and developing a straight equation for permeability as a function of petrophysical properties for carbonate reservoirs is a complicated task because of the availability of several uncertainties (e.g., anisotropy and heterogeneity) in these reservoirs. Thus, many researchers used core analysis measurements to predict permeability [
20,
21,
22]. Therefore, rock permeability is impacted by several factors, including the form, size, and shape of the grains that constitute the reservoir formations. These factors and the permeability of the rocks are determined using laboratory tests on the core samples. Laboratory tests are usually rare and expensive. In addition, these tests do not represent the whole reservoir because the core samples are taken at specific depths and there is a lack of vertical permeability measurements.
Thus, it is crucial to develop a relationship to predict the vertical permeability based on available core measurement data. To do so, a relationship between microscopic core data and microscopic-level attributes using the concept of hydraulic mean radius should be developed. Permeability–porosity variations in this concept are characterized as ().
The groundbreaking equation for predicting the permeability–porosity relationship was developed by Kozeny [
22]. Carman [
23] modified the Kozeny equation and developed the modified Kozeny–Carman equation [
24]:
where
K is permeability,
is effective porosity,
is 3D pore shape factor,
is tortuosity, and
Sgr is specific surface area per unit volume.
Equation (1) can be rearranged in terms of the hydraulic mean radius variations (
) by dividing both sides by effective porosity
and taking the square root for both sides as follows:
In addition, to simplify the prediction of the permeability–porosity relationship, the reservoir quality index (
RQI) has been proposed as a function of hydraulic mean radius variations (
). The
RQI model considers the pore, pore throat, and grain distributions with other macroscopic parameters [
25]:
Based on the
RQI model, Zahaf and Tiab [
26] developed a generalized equation for predicting vertical permeability (
Kv) as a function of hydraulic mean radius by substituting horizontal permeability (
Kh) instead of absolute permeability (
K):
where
A and
B are specific field case coefficients.
However, all previous equations for predicting vertical permeability were developed for sandstone reservoirs. Even though the uncertainties in carbonate reservoirs in terms of anisotropy and heterogeneity are well known, the prediction of the vertical permeability of carbonate reservoirs has not received sufficient attention. Thus, in this paper, new sets of correlations have been developed to predict vertical permeability for carbonate reservoirs as a function of horizontal permeability and porosity using 369 laboratory measurements.
2. Stratigraphic Description and Facies Analysis
The Arabian Plate is positioned to the north of the Mesopotamian Basin. In the west and east, respectively, this foreland basin empties into the folded zone of the Arabian stable shelf’s Salman zone. The rocks in the basin’s basement are metamorphic Precambrian rocks. The majority of the infracambrain to quaternary rocks make up the sedimentary cover, which thickens toward the east. The extensional and compressional tectonic stages are the two distinct stages of the Cretaceous tectonic growth of the Mesopotamian Basin. A passive global border sedimentary setting was formed by the South-Neo-Tethys between the Late Tithonian and the Early Turonian. The closure created a compression tectonic stage between the Late Turonian and Maastrichtian, which led to the Arabian Plate entering the foreland basin evolution stage. Folds with basement involvement are known as anticlines, and they are found in the southern and central Mesopotamian Basins [
27].
The southern Mesopotamian Basin, a transitional region between the Arabian Plate and the Zagros Mountains, is where the region of study is located. Due to various forces (e.g., climatic, eustatic, and tectonic), more than 3500 m of Cretaceous deposits have settled across this region. During the period of the Cretaceous, the region was a portion of a wider limestone layer that was situated along the Arabian Plate’s northeast passive border. The Mesopotamian basin’s Cretaceous succession consists of six supersequences. They are all separated from one another by maximum flooding surfaces and unconformities. Limestone dominated the studied carbonate formation, which is called the Mishrif Formation [
28]. The Mishrif Formation is made up of three third-order sequences stratigraphically. Each of them is made up of an upper regressive and a lower transgressive hemicycle. Either conformity or disconformity layers separate each of the third-order sequences, and mudstone and wackestone are the predominant grains on all maximum flooding surfaces. This formation is divided into three members. These members can be further divided into six zones, six subzones, and 7–8 strata. The three primary pore types created in this formation are microfracture pores, matrix pores, and dissolved vug. The formation is split into three 3rd-order sequences. There are algal, bioclastic, rudistid, and foraminiferal-rich facies in the heterogeneous carbonates of the Mishrif Formation. They are split into two permanent restrictive cycles by the facies growth and an unconformity surface at the region level. According to current biostratigraphic research, the formation is early–late Cenomanian, but the dating system places it in the middle of the Centaurus–early Turonian period. Along the Mishrif platform, the border between Iran and Iraq exhibits a high energy trend with a thickness of about 400 m [
27,
28].
Seven facies, spanning from shallow open marine Tidal Flat, are recognized that correspond to palaeo-depositional environments. The paleo-geomorphy decreased toward the north when the formation as a whole was being deposited [
29]. The paleo-geomorphical trend is controlling the decline in Mishrif’s thickness especially around the carbonate ramp and the reef-rimmed carbonate platform. Considering the specific palaeo-geographical setting throughout this basin, rudist congregations might move lateral into shoal facies or be displaced by them, which occupied the platform margin in this depositional model without forming a continuous rim. The historic slope break on this limestone platform is its most noticeable feature. The relationships between these facies are outlined below [
30,
31].
2.1. Deep Marine Facies
The micritic and argillaceous limestones that make up the units of deep marine facies range in thickness from 1 to 40 m, and they are accompanied by thin shales. Limestones are made up of fine-grained bioclasts that have been carried from nearby shallow-shelf regions as well as wackestones and mudstones that are rich in pelagic foraminifera. Deposition occurred beneath the base of the storm wave as well as the photic zone. Both the photic zone and the storm wave base were below which deposition occurred.
2.2. Shallow-Shelf Open Marine
There were six identified microfacies in this facies type of Mishrif Formation. The thick and massive units of limestone that make up the faces are made up of bioclastic wackestones and packstones. Green algae, mixed rudist, bivalves, and echinoids debris, as well as benthonic and planktonic foraminifera, are examples of bioclasts. Algae and colonial corals are also found. Shallow-shelf open-marine facies become increasingly rudist- and grain-rich as they move upward. The rudist debris was carried downslope and deposited as packstones and floatstones, suggesting that it originated from adjacent rudist biostromes. This facies association contains high porosity, grain-rich bioclastic limestones.
2.3. Rudist Biostromes
Throughout rudstones and packstones, para-autochthonous and autochthonous rudists predominate in rudist biostromes. Rudist bioclasts tend to be highly compacted, usually broken but rarely rounded, and poorly arranged. Rudists in the Mishrif Formation lack characteristics of reefs or build-ups; as a result, rudist structures are called lithosomes, and dense faunal groups are called congregations.
2.4. Shoal Facies
Grain-dominated packstones and grainstones make up the majority of the shoal facies. Rudists, echinoderms, algae, pellets, and coated grains are examples of bioclasts. Small benthonic foraminifera, mollusks, and allochthonous corals are also found locally, although their poor preservation suggests lengthy transport routes.
It is believed that the shoal facies deposits contain very little clay because they were deposited above the fairweather wave base in high energy settings.
2.5. Back-Shoal Facies
This type is varying from bioclastic and foraminiferal-rich wackestones to packstones formed of shoals or rudist biostromes during comparatively calm water conditions. In addition to the less common elements such as pellets, algae, echinoderms, gastropods, and benthonic foraminifera, bioclasts consist of fragments of coarse rudist and Chondrodonta bivalve. It is believed that the back-shoal facies originated in photic environments with low energy and shallow water.
2.6. Lagoonal Facies
Lagoonal facies, which are found in specific regions of the Mishrif, are identified through a range of benthonic foraminiferal species, such as textulariids, milolids, and alveolinids. Spicules of sponge, rudist pieces, gastropods, and ostracods are found locally and exhibit minimal breakage, suggesting short transport distances. Deposits of lagoon facies can have a thickness of over 30 m.
Thin layers of grainstones and peloidal and rudistid packstones may be found within the lagoonal successions. These are believed to be small-scale biostromes and shoals that developed in the lagoonal environment.
2.7. Tidal Flat Facies
The top part the Mishrif contains tidal flat facies, which are composed of lime mudstones with poor fossil preservation and Charophyta wackestones.
Deposits on tidal flat facies are frequently highly recrystallized, oxidized, and bleached. Although facies are difficult to distinguish from lagoonal deposits, thicknesses vary from 1 to 15 m.
3. Materials and Methods
A total of 369 experimental samples from various depths of a carbonate formation were collected from an Iraqi southern oil field. For all 369 samples, the porosity, horizontal permeability, and vertical permeability were recorded as a function of depth. All the samples were loaded for Soxhlet cleaning using refluxing hot toluene and methanol. Samples are deemed clean when no further evidence of hydrocarbons (fluorescence under UV light) or salts (silver nitrate precipitation) is present. After cleaning, the samples were then dried in a regular oven at 116 °C. Once dried, the plug samples were stored in a desiccator and allowed to cool to room temperature prior to the analysis. The weight of each sample was determined. The direct grain volume of all plug samples was then determined using Core Laboratories Ultra-poreTM 400 (Core Laboratories, Houston, TX, USA). The grain density is determined by dividing the dry weight of the plug by the grain volume determined by the Ultra-poreTM 400. Ultra-poreTM 400 uses Boyles’s law to determine the pore or grain volume from the expansion of a known mass of helium into a calibrated sample holder. The equation used for the calculation of grain volume is derived from the basic Boyle’s law as follows:
which can be expressed in term of grain volume:
where
Vg is the grain volume of the sample,
Vm is the matrix cup volume (sample holder),
Vr is the known reference volume (filled by helium),
P1 is the initial pressure, and
P2 is the expanded pressure.
The bulk volume (
Vb) is calculated by measuring the dimensions of the sample. Then, pore volume (
Vp) and porosity (ø) are calculated as follows:
Pore volume measurements were made using the Boyles law method combined with an advanced calibration technique. The effective porosity of each sample was then calculated using the pore volume determined at the designated confining stresses by the CMSTM—300 and the direct grain volume from the Ultra porosimeter. The porosity determined by this technique has the advantage of automatically accounting for bulk volume compressibility due to pore volume reduction at overburden stress.
Based upon the sample condition, the shape and dimensions for horizontal and vertical plug samples were deemed suitable for permeability testing, and all of them were selected for CMS TM-300 measurements. The direct pore volume, permeability to air, and Klinkenberg permeability of plug samples at ambient conditions were determined using Core Laboratories CMSTM-300. CMSTM-300 is an automated computer-controlled unsteady state pressure decay permeameter and porosimeter that measures pore volume, non-reactive liquid permeability (Klinkenberg), calculates an equivalent air permeability at a specified mean pressure, uses the Forchheimer Inertial Factor (Alpha and Beta), and uses the Klinkenberg slip factor at programmable, sequential confining stresses from 800 psig (minimum confining stress) to 10,000 psig on 1 and 1 ½ inches in diameter.
A transient pressure fall-off technique is used by the instrument to evaluate permeability. Helium is poured into a tank whose capacity is precisely known. Helium is released into the environment through the hydrostatically stressed core sample. A calculation of the gas flow rate and pressure drop across the core sample at any given time is made possible by tracking the decaying tank pressure over time. This measurement’s data are used to calculate the equivalent air permeability at a given mean pressure, the Klinkenberg slip factor, and the Forchheimer inertial factors (alpha and beta).
The statistical summary of the measured porosity (
), horizontal permeability (
Kh), and vertical permeability (
Kv) of the studied samples (i.e., 369 experimental samples) is shown in
Table 1.
4. Results and Discussion
In this section, we developed new correlations to predict vertical permeability for carbonate reservoirs as a function of horizontal permeability and porosity using a nonlinear regression approach. Firstly, we developed a new correlation for predicting vertical permeability for the whole data (369 samples;
Table 1) as a function of horizontal permeability and porosity. The results indicate that the best form representing the relationship of vertical permeability with horizontal permeability and porosity is found to be as follows:
This new correlation is able to predict the vertical permeability with correlation coefficients (RSQ) of 0.853.
Figure 1 shows the vertical permeability (measured and estimated from the developed correlation; Equation (9)) as a function of horizontal permeability and porosity.
Figure 1 indicates that the new correlation estimates the vertical permeability with good accuracy.
In order to increase the accuracy of the prediction vertical permeability, the used data (369 samples) were divided into four groups depending on the anisotropy (
Kv/
Kh) values, which are less than 0.1, 0.1–1, 1–10, and more than 10. Subsequently, a new correlation for predicting the vertical permeability of each data group as a function of the horizontal permeability and porosity has been developed. The following correlations were identified for each group:
The range of input and output variables that are used to develop the new correlations (Equations (10)–(13)) is summarized in
Table 2. To check the accuracy of our established new correlations, the output data for prediction were compared with the experimental results. The correlation coefficient (RSQ) for the new correlations was 0.826 for group 1 data with (
Kv/
Kh) > 10 (Equation (10)), 0.963 for group 2 data with 1 < (
Kv/
Kh) < 10 (Equation (11)), 0.947 for group 3 data with 0.1 < (
Kv/
Kh) < 1 (Equation (12)), and 0.751 for group 4 data with (
Kv/
Kh) < 0.1 (Equation (13)). The statistical analysis above reveals high values of the correlation coefficient, indicating the high accuracy of our developed correlations (Equations (9)–(12)). Furthermore,
Figure 2,
Figure 3,
Figure 4 and
Figure 5 cross-plot the experimentally measured vertical permeability versus those determined from new correlations (Equations (9)–(13)). The highly demonstrated agreements suggest that our proposed correlations are able to predict the vertical permeability for carbonate reservoirs as a function of horizontal permeability and porosity with excellent accuracy (
Figure 2,
Figure 3,
Figure 4 and
Figure 5) and
Table 3.