Abstract

Conventional departures from Archie conditions for petrophysical attribute delineation include shaliness, fresh formation waters, thin-bed reservoirs, and combinations of these cases. If these departures are unrecognized, water saturation can be overestimated, and this can result in loss of opportunity. Wireline logs of four (4) wells from Apete field were studied to delineate petrophysical attributes of shaly-sand reservoirs in the field. Shale volume and porosities were calculated, water saturations were determined by the dual water model, and net pay was estimated using field-specific pay criteria. Ten sand units within the Agbada formation penetrated by the wells were delineated and correlated and their continuity was observed across the studied wells. The reservoirs had high volume of shale (Vcl), high hydrocarbon saturation, low water saturation, and good effective porosity ranging 12.50–46.90%, 54.00–98.39%, 1.61–46.0%, and 10.40–26.80%, respectively. The pay zones are relatively inhomogeneous reservoirs as revealed from the buckle’s plot except in Apete 05. The direction of deposition of the sands was thus inferred to be east west. Empirical relationships apply with variable levels of accuracy with observation of the porosity-depth, water saturation-depth, and water saturation-porosity trends. Core data is recommended for better characterization of these reservoirs.

1. Introduction

Shales can cause complications for the petrophysicist because they are generally conductive and may therefore mask the high resistance characteristic of hydrocarbons [1]. Several factors are to be considered when delineating petrophysical attributes for shaly-sand reservoirs because clay minerals add conductivity to the formation especially at low water saturations.

Clay minerals attract water that is adsorbed onto the surface, as well as cations (e.g., sodium) that are themselves surrounded by hydration water. This gives rise to an excess conductivity compared with rock, in which clay crystals are not present, and this space might otherwise be filled with hydrocarbon. Using Archie’s equation in shaly sands results in very high water saturation values and may lead to potentially hydrocarbon bearing zones being missed. Moreover, in clean sands, the irreducible water volume is a function of the surface area of the sand grains and therefore the grain size, but for shaly sands the addition of silt and clay usually decreases effective porosity due to poorer sorting and increases the irreducible water volume with the finer grain size [2]. Archie’s equation was developed for clean rocks, and it does not account for the extra conductivity caused by the clay present in shaly sands. Therefore, Archie’s equation would not provide accurate water saturation in shaly sands. In fact, water saturations obtained from Archie’s equation have a tendency to overestimate the water in shaly sands. Several models have been proposed by many researchers for shaly-sand analysis such as Juhasz model, dual water model, Indonesian model, Waxman and Smits model, and so forth.

2. Synopsis of the Geology

The stratigraphic sequence of the Niger Delta comprises three broad lithostratigraphic units, namely, a continental shallow massive sand sequence, the Benin Formation, (2) a coastal marine sequence of alternating sands and shales, the Agbada Formation, and (3) a basal marine shale unit, the Akata Formation (Figure 2). Outcrops of these units are exposed at various localities (Figure 1). The Akata Formation consists of clays and shales with minor sand intercalations. The sediments were deposited in prodelta environments. Petroleum in the Niger Delta is produced from these unconsolidated sands in the Agbada Formation. Characteristics of the reservoirs in the Agbada Formation are controlled by depositional environment and by depth of burial. The sand percentage here is generally less than 30%. The Agbada Formation consists of alternating sand and shales representing sediments of the transitional environment comprising the lower delta plain (mangrove swamps, floodplain, and marsh) and the coastal barrier and fluviomarine realms. The sand percentage within the Agbada Formation varies from 30 to 70%, which results from the large number of depositional offlap cycles [3]. A complete cycle generally consists of thin fossiliferous transgressive marine sand, followed by an offlap sequence which commences with marine shale and continues with laminated fluviomarine sediments followed by barriers and/or fluviatile sediments terminated by another transgression [4, 5].

The Benin Formation is characterized by high sand percentage (70–100%) and forms the top layer of the Niger Delta depositional sequence. The massive sands were deposited in continental environment comprising the fluvial realms (braided and meandering systems) of the upper delta plain. The Niger Delta time-stratigraphy is based on biochronological interpretations of fossil spores, foraminifera, and calcareous nonnoplankton. The current delta-wide stratigraphic framework is largely based on palynological zonations labeled with Shell’s alphanumeric codes (e.g., P630, P780, and P860). This allows correlation across all facies types from continental (Benin) to open marine (Akata). There have been concerted efforts, within the work scope of the stratigraphic committee of the Niger Delta (STRATCOM), to produce a generally acceptable delta-wide biostratigraphic framework [9] but not much again has been accomplished after several data gathering exercise by the committee. The sediments of the Niger Delta span a period of 54.6 million years during which, worldwide, some thirty-nine eustatic sea level rises have been recognized [10]. Correlation with the chart of Galloway [11] confirms the presence of nineteen of such named marine flooding surfaces in the Niger Delta. Eight of these are locally developed. Adesida et al. [10] defined eleven lithological mega sequences marked at the base by regional mappable transgressive shales (shale markers) that are traceable across depobelt boundary faults and proposed these as the genetic sequences that can be used as the basis for lithostratigraphy of the Niger Delta.

3. Methodology

Composite wireline log data from four well logs were interpreted. The basic analysis procedure used involves the following steps; each of which is described in the following sections.

3.1. Import and Well Log Data

The well data was imported into the software used and well log correlation (Figure 3) was done after which the petrophysical attributes were delineated. Well correlation helped in determining the direction of thickness of sand being mapped and the lateral continuity of these reservoirs.

3.2. Zoning and Point Selection

Zoning is of vital importance in the interpretation of well logs. The logs were split into potential reservoir zones and nonreservoir zones. Hydrocarbon bearing intervals were identified and differentiated based largely on the readings from the shallow and deep reading resistivity tools. However, hydrocarbon typing (oil and gas differentiation) was based on density-neutron logs overlay.

3.3. Compute Shale Volume from the Gamma Ray

This was derived from the gamma ray log first by determining the gamma ray index [12]: where = gamma ray index; = gamma ray reading of the formation; = minimum gamma ray reading (sand baseline); = maximum gamma ray reading (shale baseline).

For the purpose of this research work, Larionov’s [13] volume of shale formula for tertiary rocks was used: sh: volume of shale and gamma ray index.

3.4. Compute Total Porosity and Shale-Corrected (Effective) Porosity

Total and effective porosity was estimated from the density, neutron, and sonic logs using Archie’s equation.

3.5. Compute Water Saturation

Water saturation was estimated using Archie’s water saturation formula and Schlumberger’s dual water model.

3.6. Estimate Net Pay

Calculate net pay using field-specific net pay cutoffs. Cutoff criteria used are water saturation < 50%, porosity > 10%, and volume of clay < 50%.

3.7. Use of Crossplots

Pickett, buckles, and neutron-density crossplots were generated to understand reservoir properties. Pickett plot which is a resistivity-porosity plot generated was used to determine saturation values alongside Archie parameters and . Porosity is calculated from the neutron porosity and density porosity logs and is plotted against the resistivity data obtained from the deep resistivity log. Porosity is plotted on the -axis with a logarithmic scale ranging from 0.1% to 100%, while the resistivity is plotted on the -axis with a logarithmic scale ranging from 1 to 100 ohm meter. In order to properly characterize the reservoir sands delineated and correlated across the studied wells, Buckles plot, a plot of Sw versus , was generated to depict whether or not the sands are at irreducible water saturation. Porosity is plotted on the -axis with a scale ranging from 0 to 40% porosity (shown in decimals), while water saturation is plotted on the -axis with a scale ranging from 0 to 100% (shown in decimals) water saturation. The scale for bulk volume water lines (grey lines) ranges from 0.01 to 0.25 and is shown as a secondary -axis. The implicit assumption in the Buckles plot approach is that the product of irreducible water saturation and porosity is constant. Empirical relationships were also established for porosity-depth, water saturation-depth, water saturation-porosity, and permeability-depth to check the trends.

4. Results and Discussions

Petrophysical attributes as porosity (effective and total), reservoir thickness (net and gross), water saturation (Archie and dual water model), and volume of shale (Tables 1, 2, 3, 4, 5, 6, and 7) were delineated in this research work. Results from petrophysical analysis revealed reservoir Sand F to be the most viable reservoir with net thickness as high as 126.50 ft. All the ten reservoirs exhibited good petrophysical attributes with high porosity and hydrocarbon saturation. The sands are shaly sands with high volume of shale resulting in overestimated water saturation in the reservoirs; dual water model was used for estimating the water saturation for better appraisal of the reservoirs due to their shaly nature.

4.1. Crossplot

Correlation analysis was performed to determine whether the petrophysical attributes (water saturation and porosity) are interdependent. Generally, the effective porosity decreases with depth (Figure 4) with high correlation coefficient except in Apete 15 where there was increase in porosity with depth. The observed reduction in depth would likely be due to the effect of compaction resulting from overburden pressure. Water saturation generally increases with depth (Figure 5) except in Apete 16. This implies that reservoirs in Apete field occur in shallow depth hence the unavailability of reservoirs as we go deeper into the subsurface. Efforts made to delineate trends for water saturation and effective porosity were marginally efficient. From crossplot results (Figure 6), effective porosity reduces with an increase in water saturation in Apete 05 and 06. Specifically, in Apete 15 and 16, there was no correlation at all, as the correlation coefficient was extremely low which suggests that no relationship exists between the two petrophysical attributes. Buckles plot for the reservoirs in the wells were generated (Figure 7); the results from these plots reveal only Apete 06 to be at irreducible water saturation as the data points align along the bulk volume of water (BVW) trend line due to the consistency of the data points. The reservoirs zones in this well are considered to be homogenous; therefore, hydrocarbon production from Apete 06 should be water free [14]; that is, the reservoirs would have a low water cut.

Pickett plot (Figure 8) reveals the reservoirs to be somewhat shaly which is observed with saturation exponent being less than 2 in the best porosity type. This is further confirmed from the Neutron-Density crossplot (Figure 9). Pickett plot was used to determine Archie parameters, tortuosity , and cementation exponent which is approximately 1 and 2, respectively. Neutron density crossplot reveals more of laminated clay in the reservoirs; this should be taken into consideration during well planning. The shale morphology generally changed from laminated to dispersed which affects saturation mixing function hence the need to use another saturation model (Schlumberger’s dual water model) which is designed specifically for shaly sands, rather than the conventional Archie’s water saturation model. Results from both water saturation models used show a wide disparity which could not have been noticed if only the conventional Archie’s model was used which could have led to bypassing some reservoirs as well as undervaluating the reserves in this field.

5. Conclusion

In the study of well logs from the Apete field, Niger Delta, it was observed that Apete 06 is the most economic well drilled in this field; apart from having the most presence of reservoirs it also has the highest net thickness of 639.50 ft (Table 7). Generally, water saturation increases with depth and porosity reduces with depth as a result of compaction. Water saturation-porosity trends cannot be emphatically established except in Apete 05 and 06 where porosity reduces with increase in water saturation. The reservoirs are shaly sands with the shales mostly occurring as laminated clays which could act as impediment to flow during production and therefore causing reservoir compartmentalization. Shale morphology changes from laminated to dispersed, thereby affecting saturation mixing functions.

Glossary

Buckle’s plot: A plot of water saturation against porosity generated to depict whether or not the sands are at irreducible water saturation ( on -axis and on the -axis)Bulk volume of water (BVW): Percentage of the total rock volume occupied by waterCore data: Set of data derived from analysis of core (rock) samplesEustatic sea level: Sea level change which occurs on a global scaleFluviatile: Pertaining or relating to rivers, found in or near riversFossiliferous: Means containing fossilsLaminated: Composed of layers bonded togetherMarine flooding surface: A surface of deposition at the time the shoreline is at its highest landward positionMorphology: Means this refers to the description of the shape of geologic featuresOfflap: The arrangement of strata deposited on the sea floor during the progressive withdrawal of the sea from landPickett plot: Plot of Archie’s saturation parameters against resistivity of water so as to estimate the water saturation in such a reservoirTransgression: Progressive movement of the sea towards landTransitional environment: Environment situated between the continental realm and the marineWater cut: Amount of water produced with oil.