This study presents results of log analysis from Mpapai well, which is located in the East Pande Block, southern coast of Tanzania. The study aimed at assessing the hydrocarbon potential of lithological units encountered during drilling of Mpapai well. To achieve the general objective, suites of wire-line logs from Mpapai well were used for the analysis. Based on wire-line logs, three types of lithology were identified which include sandstone, shale and limestone. Seven sandstone bodies marked as MpapaiA, B, C, D, E, F and G were identified with their tops and bases at the depth interval from 3004 m to 4008 m. Four zones among seven sandstones bodies marked as MpapaiB, E, F and MpapaiG were identified as reservoir zones. Computed petrophysical parameters for the four reservoir zones gave an average total porosity ranging from 14% to 21% with low permeability in the range of 3.92 mD to 13.67 mD. The low permeability indicates that the reservoir sand bodies are impermeable, that might have been affected by the geology of the area where high content of clay minerals reduces permeability due to filling in open spaces. The fluid type defined in the reservoir zones is basically water with high saturation greater than 75%, which indicates that the proportion of void space occupied by water is high, consequently low hydrocarbon saturation and production. Despite of fair to good porosity, the low permeability and high-water saturation indicate that the quality of Mpapai prospect is poor.
Southern coastal basins of Tanzania are located in the southern tip of the Somali basin, which is connected to the Natal basin in South Africa through the Mozambique Channel [
The East Pande Block belongs to Pande Formation (Late Cretaceous to Early Oligocene) of the Kilwa Group that consists of both onshore and offshore parts of the Southern coastal basins of Tanzania. Other formations forming the Kilwa group are the Nangurukuru (Late Cretaceous to Paleocene), Kivinje (Paleocene to Lower Eocene) and Masoko (Middle Eocene) formations [
to a total depth of 4153 m. The preliminary petrophysical analysis of well logs indicated that gas traces were encountered in the upper strata of Mpapai well. The presence of water and gas traces encountered at 781 m of water depth suggests that the East Pande Block have possible petroleum system. This suggests that the petroleum system in the East Pande Block contains matured source rock somewhere in the basin. Despite the presence of hydrocarbon traces within the Mpapai well, little is documented about petrophysical properties of reservoir rocks. Lithological units that constitute the East Pande Block are not well known. In addition, the sequence and thickness of sedimentary units, which are potential for petroleum system, are also not well established.
Petrophysical analysis of well logs is one among the most useful and important tools for reservoir rock characterization. It helps to define physical characteristics of rocks such as lithology, porosity, permeability and fluid saturation. The analysis is also useful in identifying potential reservoir intervals, distinguishing the type of fluid in a reservoir and estimating hydrocarbon reserves [
The East Pande Exploration Block is one among the south-eastern Tanzania coastal basin located along the passive continental margins of western Indian Ocean developed on the Precambrian Pan-African basement [
The bedrock of the Kilwa Peninsula area is composed of Upper Cretaceous to Lower Miocene claystone and clays, where the Late Campanian to Middle Eocene aged clays are assigned to the Kilwa Group [
The East Pande formation (Late Cretaceous to Early Oligocene) (
The geology of the area is strongly affected by two-stage break-up of Gondwana,
which began at about 300 Ma [
The rifting and spreading was accompanied by a major marine transgression onto the passive continental margins of both West and East Gondwanan blocks. In Tanzania, this resulted in thick Upper Jurassic and Lower Cretaceous sedimentation in a series of marginal basins, which include the Mandawa and Ruvuma Basins. By the end of the Early Cretaceous, the axis of spreading had jumped east of Madagascar which resulted in a period of stabilization along the East African margin. A Late Cretaceous transgression resulted in widespread deep-water clay and mud facies accumulating in offshore basins, and similar conditions continued into the Paleogene. Thick clay deposits are particularly well developed along coastal Tanzania. At about 35 Ma the current East African Rift system had been initiated and was followed by renewed tectonic movement along the Davie and Mozambique submarine ridges [
The present sequence stratigraphic framework and architecture of offshore Tanzania basins, spanning the post break-up Jurassic to Neogene has been well established by (
The stratigraphy of the East Pande exploration Block is defined by Megasequence K1 to Pg2 in which, Megasequence K1 (Berriasian-Aptian) is not well resolved, this is due to limited well penetrations offshore for seismic calibration, and poor biostratigraphic dating. The K2 Megasequence (Early Albian-Early Campanian) records the development of a thick, muddy slope characterized by
the presence of long-lived slope valley systems in the architecture of Tanzania. MS K3 records the complete infilling of canyons in Mid-Late Campanian times, by back stepping, ponded sheet-systems, comprising high-density turbidity and hybrid flow deposits. The Base Tertiary unconformity marks another major change in depositional style at the base of Megasequence Pg1 (Paleocene-Eocene). MS Pg1 is divided into three parts based on changes in depositional style and unconformities recorded in offshore well data. Ng1 Megasequence is characterized by the onset of major slope failures, which result in plugging, and disruption of the slope channel systems. A renewed phase of slope failure and extensive fan deposition represents the Pliocene-Recent (Ng2 Megasequence) and numerous deep canyons still characterize the slope today.
The study accomplished the objectives using well log data from Mpapai well provided by Tanzania Petroleum Development Company (TPDC). Well logs include gamma ray log, photeolectric factors (PEF), resistivity logs which include shallow, medium and deep resistivity, density log, neutron-porosity log and sonic log. Well logs from Mpapai well were used to analyze the petrophysical properties of the reservoir rocks. Both qualitative and quantitative interpretations were conducted using different well log data. The wireline logs used were gamma ray log, resistivity logs, photoelectric effect log (PEF), sonic log, neutron and density log. Log interpretation was conducted using Techlog software version 2013.4.0.1.
The identification of lithology is fundamental to all reservoir characterization because the physical and chemical properties of the rock that holds hydrocarbons and/or water affect the response of every tool used to measure reservoir properties [
Identification of lithology from gamma ray log was done by reading API values on the gamma ray curve which ranges from 0 to 150 API where the sand/shale baseline was inserted at 75 API in order to differentiate the two lithologies (sand and shale). Therefore, all formations with gamma ray values less than 75 API were classified as sandstone, while those with gamma ray values greater than 75 API were classified as shale.
In this study, classification of lithology was based on the values given in
Neutron and density logs were sketched together in the same track for doing some comparison and the correlation between the two curves leads to better lithology identification. When both the density and neutron logs show lower value, it indicates sandstone formation whereby overlay of the two log curves shows limestone lithology and the increase of neutron and density value indicates shale lithology.
Lithology | Sandstone | Dolomite | Shale | Limestone | Anhydrite |
---|---|---|---|---|---|
PEF value barns/electron | 1.8 - 2 | 3.14 - 4 | 1.8 - 6 | 5 | 5 |
Neutron-density cross plots are easy to use for clean (non-shaly) reservoir rocks. The plots are entered with a bulk density and apparent neutron porosity. A rock type (sandstone, limestone, or dolomite) and a corrected porosity can be read from the cross plot. The neutron-density cross plot was used to determine lithology of Mpapai well by plotting together neutron and density logs and use gamma ray log for scale. Points were observed to fall into different lithological regions and the interpretation of the lithology (sandstone, carbonate and shale) was based on color legend bar, which is the intensity of gamma rays. The gamma ray intensity represents the amount of concentration of the radioactive elements present in the neutron-density logs, thus giving different types of lithology.
Reservoir is the only zone, which is potential for economic interest because it contains storage space for fluid (hydrocarbon or water) to accumulate. Thus, must first of all be identified in order to evaluate important parameters suitable for hydrocarbon exploration. Reservoir identification was conducted after the interpretation of various lithology of Mpapai well. The clean (non-shale) formation was marked in different zones as reservoir rocks. Reservoir rocks are defined as subsurface pool of hydrocarbon or water contained in porous or fractured rock formation [
Porosity and permeability are the most important physical properties of the reservoir rocks. For rock to be named a reservoir has to be porous and permeable. Sandstone (which covers 62% of the petroleum reservoir) and Limestone are two sedimentary rocks, which are used as reservoir rocks [
Gamma ray log was used in the identification of reservoir rock based on the fact that, sandstone reservoir exhibits very low radioactivity because of low content of radioactive elements hence have low gamma ray value and the log deflect to the left of shale/sand baseline [
It is very important to identify the type of fluid in a reservoir rock, because reservoirs may contain hydrocarbon (oil and gas), non-hydrocarbon fluid (water) or both. For a reservoir to contain hydrocarbons the zone should be porous with resistivity values higher than those of water-bearing zones [
The analysis of shale volume was conducted in order to determine the amount of shale or radioactive minerals, contained in the reservoir rock in order to delineate zones of interest. In this study, the volume of shale was estimated using Gamma ray logs and neutron-density methods. The method that gave the minimum value of shale volume was chosen for further estimation of porosity and water saturation in order to minimize the effect of reducing reservoir quality caused by high content of shale volume.
Gamma-ray log is one of the best tools used for identifying and determining the shale volume. This is principally due to its sensitive response on the radioactive materials, which normally concentrate in the shaly rocks. In this study, the first step to estimate the reservoir’s shale volume of Mpapai well from gamma ray logs was to calculate the Gamma Ray Index by using the Dresser Atlas 1979 formula in equation below:
I G R = G R log − G R min G R max − G R min (1)
where I G R is the gamma ray index, G R log is picked log value for each zone, G R min and G R max indicate gamma ray values picked in sand and shale base lines respectively.
The gamma ray log has two empirical responses. The linear and non-linear responses, the nonlinear responses are based on geographic area or formation age, which includes the Larionov model [
V S H ( 1969 ) = 0.083 ∗ ( 2 ( 2 ∗ 3.7 I G R ) − 1 ) (2)
The separation between neutron and density porosity is a common method for calculating shale volume. This method is accurate only when the shaly sand contains pure quartz plus clay minerals. Most sandstones are not pure quartz, but rather may contain other minerals like Micas, siderite, ankerite, and volcanic rock fragments. These minerals are heavier than quartz, which cause excess separation by reducing density porosity and increasing neutron porosity [
V s h = ϕ n − ϕ d ϕ N s h − ϕ D s h (3)
where V s h = Volume of shale, ϕ n = Neutron porosity in sand, ϕ d = Density porosity in sand, ϕ N s h = Neutron porosity in shale portion or adjacent shale and ϕ D s h = Density porosity in shale portion or adjacent shale.
Porosity is very important parameter of the reservoir rock as it is used to describe the amount of open space filled with fluid (hydrocarbon or water). In this study, porosity was calculated from two methods, density log and density-neutron combination logs. The criteria used in classifying porosity are given in
Density log is a good method for determining either total or effective porosity in single or multiple mineral fluid-filled reservoirs. The method of estimating porosity from the density requires determining the matrix density (ρma), the density log reading (ρb), and the fluid density (ρfl) at the depth of interest. The matrix density is determined by the lithology. Normally, sandstone is 2.65 g/cm3, limestone is 2.71 g/cm3, and dolomite is 2.87 g/cm3 [
Average Porosity (ϕ) values | ϕ < 0.05 | 0.05 < ϕ < 0.1 | 0.1 < ϕ < 0.15 | 0.15 < ϕ < 0.25 | 0.25 < ϕ < 0.30 | ϕ > 0.30 |
---|---|---|---|---|---|---|
Qualitative description | Negligible | Poor | Fair | Good | Very Good | Excellent |
In this study, the porosity was also determined from density log using the formula below:
ϕ = ρ m a − ρ b ρ m a − ρ f (4)
where ϕ = Total porosity, ρ m a = Matrix density (or grain) density (2.65 g/cm3), ρ b = Formation bulk density from log and ρ f = Fluid density (1.1 g/cm3).
The combination of neutron and density logs provides a good source of porosity data, especially in formation of complex lithology. Better estimates of porosity are possible with this method than using other tool separately such as density and sonic because inferences about lithology and fluid content can be made. Porosity from Neutron-Density log can be calculated mathematically using the following equation [
ϕ N − D = ( ϕ N 2 − ϕ D 2 2 ) 1 2 (5)
where ϕ N − D = Neutron-density porosity, ϕ N = Neutron porosity and ϕ D = Density porosity.
In this study, two models (the Archie’s and Indonesian models) were used to calculate water saturation of the reservoir rocks and results were compared. The Archie’s model works well in homogeneous or clean sand reservoir [
Archie’s equation is most famous method for calculating water saturation in clean and Shale free formation. Archie formula is based on the fact that the only conductive material in the formation is salt water; but in a sandy shale formation, ions that are released along with shale are also responsible for conducting electrical current [
S w = ( a ∗ R w ϕ m ∗ R t ) 1 n (6)
where S w = Water saturation, R w = Formation water resistivity, ϕ = Total porosity, R t = True formation resistivity, a = Tortuosity factor, m = Cementation exponent and n = Saturation exponent.
Water saturation ( S w ) results from the formula are comparatively easy to calculate and because it is not a quadratic equation, it gives results that are always greater than zero. Calculation of water saturation by using this method depends on porosity, shale volume and resistivity, water and deep resistivity. Water saturation from Indonesian model is given by the following formula:
S w = { [ ( V s h 2 − V s h R s h ) 1 2 + ( ϕ e m R w ) 1 2 ] 2 R t } − 1 n (7)
where V s h = Volume of shale, R s h = Resistivity of shale, S w = Water saturation, ϕ e = Effective porosity, R w = Water resistivity of formation, m = Cementation coefficient, n = Saturation capacity and R t = Real resistance.
Hydrocarbon saturation S h is the percentage of pore volume in a formation occupied by hydrocarbon. In this study, the hydrocarbon saturation was determined by subtracting the value of water saturation from 100%, as illustrated in the equation below:
S h = ( 100 − S w ) % (8)
Well log is one among the methods used to estimate permeability of a reservoir rock. There are two principal log measurements used to estimate permeability, these are resistivity and porosity logs. Porosity log is frequently preferred than resistivity log because it is strongly correlated to permeability [
K = 0.136 ϕ 4.4 S w i r 2 (9)
where K = Permeability in mD, ϕ = Porosity and S w i r = Irreducible water saturation Irreducible water saturation was estimated from Crain’s method [
S w i r r = ϕ × S w ϕ e f f (10)
where ϕ × S w represents the bulk volume of water (BVW).
The bulk volume of water is used to indicate if the reservoir rock is at irreducible water saturation or not. If values for bulk volume of water obtained at different depths are constant, they indicate that the reservoir zone is at irreducible water saturation otherwise the reservoir zone is not at irreducible water saturation [
Petrophysical analysis of well logs and cross plot reveals that Mpapai well consists of three types of lithology: sand, shale and carbonate (limestone). The well section was divided into four depth intervals with their corresponding log curves for easy description of lithology type. In connection to the four depth intervals used to describe the three types of lithology, the same intervals were used to form four zones that consist of alternating sand and shale formations (Figures 4-7;
From the four chosen depth intervals shown in the previous section, seven clean sand bodies were identified across Mpapai well, which were named as MpapaiA, B, C, D, E, F and MpapaiG. Generally, based on visual observation of well logs, zones that showed low gamma ray values, relatively high porosity and high resistivity values were identified as reservoir zones. Therefore, only four sandstone bodies named MpapaiB, E, F and MpapaiG among seven identified zones were marked as reservoir zones (
The neutron-density logs combination and resistivity logs were used for the identification and characterization of various fluids in the reservoir zone. Based on visual observation of these logs four selected reservoir zones named MpapaiB, E, F and MpapaiG among seven selected reservoir zones were identified as gas bearing zones. This is due to the presence of neutron-density crossover and
Average k values (mD) | 1.0 - 15 | 15 - 50 | 50 - 250 | 250 - 1000 | >1000 |
---|---|---|---|---|---|
Qualitative description | Poor to fair | Moderate | Good | Very Good | Excellent |
Depth | Sandstone Bodies | Reservoir Zones | Lithology |
---|---|---|---|
3003.66 - 3008.04 | MpapaiA | Sand | |
3014.99 - 3018.79 | MpapaiB | MpapaiB | Sand |
3079.88 - 3085.2 | MpapaiC | Sand | |
3119.53 - 3125.8 | MpapaiD | Shaly Sand | |
3199.37 - 3235.95 | MpapaiE | MpapaiE | Shaly Sand |
3246.74 - 3252.29 | MpapaiF | MpapaiF | Sand |
3968.57 - 4007.26 | MpapaiG | MpapaiG | Sand |
high resistivity values. Neutron and density crossover were observed in some intervals and marked by yellow color as shown in the log curves for MpapaiB, E, F and MpapaiG respectively (see Figures 10-13). In resistivity logs, values in the reservoir zones were observed to be relatively higher, which also indicate the presence of hydrocarbon. As described in methodology part, resistivity logs are commonly used to differentiate types of hydrocarbon fluids in the sense that liquid hydrocarbon normally display higher resistivity values compared to gas zones. Based on these observations, the type of hydrocarbon fluid that could be found in these reservoirs is gas.
Four reservoir zones have been selected in Mpapai well with average thicknesses ranging from 3.8 m in MpapaiB to 38.69 m in MpapaiG. The computed
average shale volume content of the reservoirs ranges between 0 and 0.079 v/v (
Zones | Top (m) | Bottom (m) | Gross (m) | Net (m) | NTG (v/v) | BVW | Av-Vsh (v/v) | Av-Por % | Av-Sw % | Av-K (mD) |
---|---|---|---|---|---|---|---|---|---|---|
MpapaiB | 3014.99 | 3018.79 | 3.8 | 2.896 | 0.762 | 0.522 | 0.055 | 20.9 | 86.3 | 13.68 |
MpapaiE | 3199.37 | 3235.95 | 36.58 | 5.181 | 0.143 | 0.552 | 0.077 | 14.2 | 75.3 | 7.068 |
MpapaiF | 3246.74 | 3252.29 | 5.55 | 1.372 | 0.247 | 0.136 | 0 | 15.5 | 64.1 | 11.16 |
MpapaiG | 3968.57 | 4007.26 | 38.69 | 26.52 | 0.685 | 4.298 | 0.079 | 16.5 | 98.4 | 3.929 |
75.3%, 64.1% and 98.4% for MpapaiB, E, F and MpapaiG respectively while that of Archie model was observed to be 84.8%, 74%, 63.2% and 97.4% for MpapaiB, E, F and MpapaiG respectively (see Figures 10-13). The water saturation values suggest that reservoir zones are water bearing with low hydrocarbon saturation (1.6% to 35.9%). The average permeability was observed to be fair to moderate with average values ranging from 3.9 mD to 11.2 mD (see
The petrophysical analysis of Mpapai well on identification of lithology indicates that the well consists of three types of lithology, which are shale, sand and very little carbonate. The most dominant lithological unit encountered at Mpapai well was shale formation with thickness ranging from 7 m to 716 m. These shale units as described in the study by Nicholas et al. [
Based on well log analysis, a total of seven clean sand formations were identified, which were named as MpapaiA, B, C, D, E, F and MpapaiG. Of the seven sandstone bodies, four reservoir zones with high resistivity values and the presence of neutron-density crossover were identified. These include MpapaiB (3.8 m thick), MpapaiE (36.58 m thick), MpapaiF (5.55 m thick) and MpapaiG (38.69 m thick). The reason for categorizing the four sand units as different reservoir zones was based on eliminating thick shale beds between reservoirs so as to reduce the effect of increasing shale volume when computing other petrophysical parameters. The shale formation between the identified reservoirs could thus be interpreted as source rock when located below reservoir rock and as a seal rock when located above the reservoir rock.
The four selected reservoir zones were analyzed quantitatively to estimate the values of shale volume, porosity and water saturation by using empirical formulas as described in the methodology part. After applying cutoff values of 0.5 shale volume, 9% porosity and 50% of water saturation, the net pay thicknesses for the four selected reservoir zones of MpapaiB, E, F and MpapaiG were found to be 2.90 m, 5.18 m, 1.37 m and 26.52 m respectively (
Generally the quality of reservoir zones of Mpapai well is strongly affected by both local and regional geology of the area. Locally, the East Pande Block is highly composed of clays with very fine grains [
Generally, by considering all parameters such as reservoir thickness, shale volume, porosity, permeability, water saturation and hydrocarbon saturation from the log analysis performed in this study, the quality of the reservoir sand units of Mpapai well is poor. Some reservoir zones (MpapaiB and F) are very thin, which also reduce the quality of reservoir zones if they do not extend lateral. Selected reservoir zones have an average porosity ranging from fair to good. MpapaiB has good porosity of 21% while MpapaiE, F and G have fair porosity of 14%, 15% and 16% respectively. Permeability obtained from this analysis is fair indicating that the reservoir sand bodies are impermeable that might have been affected by the geology of the area where high degree of compaction, cementation (presence of carbonate cement in sandstone formation) and high content of clay minerals reduces permeability due to filling in open spaces. The quality of the reservoir zones is also poor due to low hydrocarbon saturation in which more than 75% of the reservoir zones are filled with water. The water saturation indicates that the proportion of void space occupied by water is high consequently low hydrocarbon saturation and low hydrocarbon production.
The authors are grateful to Statoil Companyin collaboration with the University of Dar es Salaam in Tanzania for funding this research as part of MSc. study. We thank Tanzania Petroleum Development Company (TPDC) for providing the data used in this study. We also thank many individuals from the University of Dar es Salaam and TPDC for their assistance in completion of this work. We also thank anonymous reviewers for helpful comments.
Lyaka, A.L. and Mulibo, G.D. (2018) Petrophysical Analysis of the Mpapai Well Logs in the East Pande Exploration Block, Southern Coast of Tanzania: Geological Implication on the Hydrocarbon Potential. Open Journal of Geology, 8, 781-802. https://doi.org/10.4236/ojg.2018.88046