Electro-sequence analysis and natural resources potential of a transitional environment in central swamp depobelt, Niger delta, Nigeria

Sequence stratigraphy and Electrosequence of the Don-Field Central Swamp depobelt of the Niger Delta Basin, was carried out using well logs, biostratigraphic and core data. The aim of this study is to re-characterise and revalidate the reservoirs within the study area and provide an optimum framework for predicting reservoir character, development and reservoir quality distribution. Detailed petrophysical evaluation reveals proficient hydrocarbon reservoir units with variable parameters (i.e. porosity, permeability and shale volume) that can accommodate productive hydrocarbon reserve. Core data analyses revealed heavy minerals like pyrites and glauconites that can be of economic interest. The gamma wireline signatures of the wells show reservoir units that are predominantly blocky in shape indicative of fluvial channels. Fluvial channels are typically associated with coastal and/or marginal marine environments. More so, the aforementioned observations were validated using palaeoenvironmental interpretation of strata penetrated by Well-1, Well-2, and Well-3 via the sidewall core analyses. Sequence stratigraphic surfaces such as maximum flooding surfaces (MFS) and sequence boundary (SB) were observed to show good correlation with palynological and foraminiferal zones. The stacking patterns (i.e. progradation, retrogradation and aggradation) encountered gave acumen into systems tracts (i.e. lowstand system tract (LST), transgressıve system tract (TST) and highstand system tract (HST) interpretation. These system tracts constitute a depositional sequence set of middle—late Miocene age bounded chronologically by MFS of 9.5 Ma, 12.8 Ma, 15.0 Ma, 15.9 Ma and SB of 12.1 Ma, 13.1 Ma, 15.5 Ma, 15.9 Ma. The SB and MFS were delineated based on Shell Petroleum Development Company (SPDC) biozonation scheme: F-zone and P-zone were within P820-P680 and F9500-F9600/F9700 as contained in the biofacies interpreted data of index fossils (Globorotalia foshi S.I; Globorotalia praefohsi; Globorotalia peripheroacuta; Orbulina spp; Praeorbulina sicana; Paragloborotalia mayeri; Neogloboquadrina acostaensis) which corresponded to the alphanumeric code N8–N16 planktonic foraminiferal zones of Blow.


Introduction
Over the years, the global demand for oil and gas products has continued to increase and this situation has intensified exploration activities in search for new hydrocarbon reserves. Based on the aforementioned reasons, the Niger Delta Basin has been extensively studied to unravel subtle stratigraphic and structural traps (Knox and Omatsola 1987;Weber and Daukoru 1975;Evamy et al. 1978;Ejedawe 1981;Urang et al. 2020). In order to meet the increased global demand for oil and gas, different geological, geophysical and geochemical approaches have been employed in hydrocarbon exploration and exploitation. However, these methods tend to be complex, expensive, and pose enormous challenge during formation evaluation and decision-making (Ladipo 1992;Stacher 1995). Concurrently, it is important to take into account the risks associated with the environment, economy and technology uncertainties during oil and gas exploration. Formation evaluation provides maximum framework for reservoir assessment geared towards successful exploration, development and exploitation of commercially viable hydrocarbon reservoirs. The sedimentological and stratigraphic facets of Africa's leading oil province have been opened-up in recent time by the works of Reijer (2011) and Mode et al. (2015) based on the pioneer researches of Short and Stauble (1967), Weber (1971), Knox and Omatsola (1987), Ladipo (1992), Stacher (1995) and Durand (1995). Previous research works on sedimentological, biostratigraphic and sequence-stratigraphic styles clarified the collective influence of eustatic cyclicity and local tectonics (Ladipo 1992;Stacher 1995;Reijers et al. 1997;Owoyemi and Willis 2006;Magbagbeoloa and Willis 2007;Anyiam et al., 2019). Thus, the genetic sequence of depositional belt is marked at the top by the maximum flooding surface and at the base by a transgressive surface. Galloway (1989a, b), based on reports of Frazier (1974) proposed that maximum flooding surfaces rather than sub-aerial unconformities be used as sequence boundaries. Thus, strata bounded by maximum flooding surfaces within transgressive shales are more readily identifiable in the Niger Delta. Individual sea-level cycles are reflected in the Niger Delta in various sedimentary sequences (Haq et al. 1988). Interferences of these cycles with different periods resulted in the formation of sequence sets that are chronostratigraphically confined and sedimentologically characterised. The aim of this study is to reflect, revalidate and establish a cost effective way by which sequence sets can provide optimum framework for predicting reservoir character and quality distribution in the central swamp depobelt, Niger Delta-Nigeria.

Location and geology of the study area
The study area is part of the Cenozoic oil prolific Niger Delta Basin according to Merki (1972). The area "Donfield'' is located in the Central Swamp Depobelt of the Niger Delta within Oil Mining Lease-12 (OML-12) block with an area of about 91,803.278 square feet. It is located between Longitudes 7°0.22 I E-7°0.32 I E and Latitudes 4°0.81 I N-4°0.85 I N (Fig. 1).
The Niger Delta is a terrigenous epicratonic embayment which lies at the edge of the Africa continental crust (Weber and Daukoru 1975;Evamy et al. 1978;Ejedawe 1981). The sedimentary basin is composed of five depositional units known as depobelts that are considered as separate entities based on the diversity in stratigraphic, structural framework and hydrocarbon distribution (Hospers 1965;Stoneley 1966;Burke et al. 1972;Doust 1987;Knox and Omatsola 1989). These depobelts are: (i) the Northern Delta of Late Eocene-Middle Miocene, (ii) Greater Ughelli of Oligocene-Early Miocene, (iii) Central Swamp of Middle-Late Miocene, (iv) Coastal Swamp of Late Miocene-Pliocene and (v) Offshore of Middle-Late Pliocene (Weber 1971;Weber and Daukoru 1975;Doust and Omatsola 1990). The deltaic deposit are composed predominantly of regressive clastic successions mainly of sedimentary origin and this regressive packages consists of deltaic deposits that are defined by three main lithostratigraphic units, i.e. Akata, Agbada and Benin formations (Fig. 2). The marine shales of the Akata formation at the base consist mainly of overpressured marine shales that are about 19,685.04 feet thick. The Akata formation is overlain by the Agbada formation that is composed of paralic successions of interbedded sands and shales, about 13,123.36 feet thick in the central parts, but however, thins seawards andtowards the delta margin. The non-marine alluvial continental depositsof the Benin formation is composed mainly of fluviatile gravels and sands that thickens up to 6561.68 feet in the delta centre and thins towards the delta margin. The Niger Delta basin currently produces hydrocarbons mainly from rocks of Eocene-Pliocene age, (Weber and Daukoru 1975;Evamy et al. 1978;Ekweozor and Daukoru 1994;Ekweozor and Okoye 1980;Haack et al. 2000). However, due to the advancement in exploration technology, investigation in the basin is gradually shifting to the delta toe and flanks of the delta, where Cretaceous prospects are expected.

Materials and methods
Formation evaluation, involves the estimation of lithological properties, reservoir potential and performance based on parameters such as porosity, permeability, grain size, degree of compaction, paleoenvironmental indices, amount of matrix and cement components. The aforementioned parameters were determined through analysis of wireline log suites, sidewall core samples and biostratigraphics data. The flow chart of the methodology is as presented in Fig. 3.

Extraction of petrophysical parameters from wireline logs
The well log data from Don-Field was made available for research by Shell Petroleum Development Company, Nigeria. An indispensable step in formation evaluation procedure is the estimation of amount of shale present in the prospective sandstone. (Ulasi et al. 2012;Urang et al. 2020). This aided the exact estimation of formation porosity and fluid content contained within the pay zone. Wireline logs were 1 3 used to ascertain depths of deposition, thickness of hydrocarbon productive intervals, specific reservoir parameters and to distinguish between fluids in the reservoir encountered using derivative models of Archie's Eqs. (1942) Suites of wireline logs ranging from gamma ray (GR), resistivity (LLD), sonic (DT), density (RHOB) and neutron (NPHI) logs were used in this study. This study considers the case of tertiary unconsolidated siliciclastic rocks where composition, pore geometry, interbedded thin shale layers, effect of invasion, pore fluid types and saturation levels, fluids distribution pattern, cracks and their orientation, grain packing, pore networks are possible sources of confidence in quantitative reservoir parameter evaluation. Owning to the aforementioned effectiveness in petrophysical evaluation their significance can sway or converse depending upon the area, geology, compaction and depositional environment. Volume of shale (Vsh) is calculated and estimated through gamma ray index (IGR) from GR log to approximate clay contents in the reservoir rock using Larionov Eqs. (1969). The porosity of the reservoir interval is calculated and corrected by using density and neutron logs derivatives. The water saturation within reservoir pores was estimated from resistivity log using Archie's Eq. (1942) for clean sand intervals. The hydrocarbon saturation (Sh) (the percentage Fig. 1 a Base map of the Don-field, central depobelt, Niger delta b Niger delta basin and its stage development.adapted from Short and Stauble, 1967 c Regional basins map of Nigeria showing Niger delta (Adapted from Corredor et.al., 2005) 1 3 of pore volume in a formation occupied by hydrocarbons) was obtained by subtracting the value obtained for water saturation from 100%. The permeability, (K) (which is the ability of a rock to transmit fluid) was calculated by Wyllie and Rose (1950) log derived permeability equation. The ratio of water saturation to flushed zone water saturation gives the amount of hydrocarbons which have been moved by the invasion process. This ratio inferred as the moveable hydrocarbon index (MHI). Similarly, the work by Akaerue, et. al. (2018) in the same study area gave comprehensive mathematical derivatives and equations which were used in the current work to compute the petrophysical parameters. Gamma ray log was used to delineate different lithologies. The qualitative identification of reservoir zones was done using the log signatures from the studied packages. Strata with high GR, high neutron values, low resistivity and low density values indicated shale and was thus eliminated. Well log patterns of the sedimentary facies were then identified and interpreted through the gamma ray and resistivity log motifs in conjunction with biofacies interpretation, consequently providing information on rock properties and depositional environment as modelled by Busch (1975), Schlumberger (1985), Mode et al. (2015) and Anyiam, et al. (2019) (Fig. 4a, b).

Core data
Core data are usually generated during or after drilling operations, especially from exploration or production wells for the purpose of obtaining different information related to subsurface geology (Keelan 1982). This provided information on the physical rock properties and facies analysis (lithology, porosity and fluid evaluation). The rock fabric and degree of compaction were considered as these affect porosity and permeability which are fundamental properties of the reservoir in determining quality and effectiveness. Shale volume was also estimated from the core analysis. Paleoenvironment of deposition of the rock and reservoir quality were determined to complement wireline log signatures. Total 65 core sidewall samples were taken in the interval of 4446 feetto 10,110 feetand analysed to measure uncertainty in petrophysical and rock physics parameters. Proper sedimentological analysis was carried out to ascertain the degree of heterogeneity in the sandstones as this could greatly affect quality.

Biostratigraphic reports
Biostratigraphic analysis of fossil material derived from core and cuttings delivers information on the age of sediments, the correlation surfaces between wells, and the depositional environment. The biostratigraphic data supplied by SPDC Nigeria Limited of Well-1 and Well-2 were used for the reservoir characterisation (Tables 1, 2, 3 and 4). A total of 105 composited samples from depth of 3340 to 9855 feet and 3700 to 10,640 feet, respectively, were processed and analysed for micropaleontological and palynological studies. Interpreted palynological (P) zone and Foraminifera (F) zone based on SPDC bio-zonation scheme were calibrated and depth matched with corresponding wireline logs. The population and diversity of the benthic and planktonic foraminifera integrated with the wireline motifs were used for the stratigraphic surfaces interpretation and ages assigned. (Fig. 4b, 5, 6).

Petrophysical analysis of Well-1 reservoirs
By applying a cutoff value of 70.1 API for the GR log, five reservoir units were delineated in Well-1 within the Agbada formation. The five reservoir units, i.e. A1 to A5 have thicknesses of 320, 160, 190, 52 and 25 feet, respectively ( Table 5).The thickness of the individual net sand units ranged from 15 to 110 feet indicating that these reservoirs are not completely clean but contain some inter-fingerings of shale (a typical characteristic of paralic succession such as the Agbada formation). Average porosity values in the reservoirs range from 16.07 to 21.51% suggesting, good to very good porosity reservoirs for the well (Rider 1986;Etuefeotor 1997). The average volume of shale (Vsh) in reservoir A1 is 5.37% suggesting clean sand reservoir (Davies 2002). This means that the average volume of shale (Vsh) value in reservoir A1 may not affect its water saturation (Sw) significantly, unlike in other reservoirs (A2 to A5) where average volume of shale was more than 10 to 15% which is the permitted threshold according to Hilchie, (1978). Average resistivity and permeability values of 26.54 and 33.87 ohms-metre, and 131.08 and 1920.22 md, respectively, of reservoirs A5 and A3 suggested that these reservoir units contain hydrocarbon with excellent permeability. Reservoirs A1, A2 and A4 have average resistivity values of 1.5, 1.96 and 1.14 Ωm, respectively. This suggests that the fluid within the reservoirs may be formation water. Also, these reservoirs have average water saturation percentage values of 79.41, 30.12 and 40.6%, respectively. Hence, suggesting high water saturation zones with little or no hydrocarbon.

Petrophysical analysis of Well-2 reservoirs
Well-2 is located at the eastern part of the `Don Field'. A cutoff value of 69.0 API for the GR log was used to differentiate reservoirs from non-reservoirs. The reservoirs identified in Well-2 (B1 to B7), ranged in thickness from 30 feet in reservoirs B1 to 269 feet in B3 (Table 6). Reservoirs have very good to excellent average porosity values which ranges from 20.71% in reservoir B7 to 22.97% in reservoir B2. Reduction in average porosity values from reservoirs B1 to B7 is probably due to compaction resulting from overlying sedimentation. Average volumes of shale recorded are between 13.79 and 29.25%, suggesting that these reservoirs are not clean. Reservoirs B2, B4 and B5 show evidence of hydrocarbon saturation with average resistivity values of 28.80, 23.19 and 13.33 Ωm, respectively. Each of the reservoirs has average Sw of 13.38, 13.82 and 20.32% which was considered poor and implies that hydrocarbon saturations is good. Based on the negative crossing of the neutron and density logs by implication the hydrocarbon in reservoir B2 is gas, and the gas-water-contact (GWC) was delineated at about 5986 feet. Reservoirs B4 and B5 are probably hydrocarbon filled with little or no formation water. This probably explains the tracking together of density and neutron logs signatures. Reservoirs B1, B3, B6 and B7 are water saturated with little or no hydrocarbon based on the low resistivity values observed. The average moveable hydrocarbon index (MHI) values of reservoirs 0.356, 0.426, 0.480 and 0.505 for the respective reservoirs, indicative of good fluid moveability.

Petrophysical analysis of Well-3 reservoirs
Well-3 is located in between Well-1 and Well-2 in the central part of the Don Field (Fig. 1a). The reservoir units' average thickness ranges between 113 feet in reservoir C7 and    1 3 211 feet in reservoir C5, while the sand thickness of these reservoirs ranged from 37 to 205 feet (Table 7). Average porosity values of 19.90% to 27.05% suggest good to excellent porosity. The average Vsh in these reservoirs ranged from 13.60 to 28.54% indicating that the reservoirs may

Depth of hydrocarbon occurrence
The depths where hydrocarbons occur in a reservoir, especially at the contact zones are marked by increase in resistivity values. Result analysis shows that the depth of hydrocarbon occurrence varies from well to well due to variations in thickness of the overlying sediments, and possibly faulting (Anyiam et. al. 2019). The sand master report of well-2 and well-3 used in the determination of hydrocarbon depth of occurrence shows the fluid contacts within the reservoir, such as, hydrocarbon water contact (HCWC), gas down to (GDT), gas up to (GUT), oil up to (OUT), oil down to (ODT), water up to (WUT) and oil water contact (OWC)  (Tables 8 and 9). In well-2, the depth of 5982 feet marks OUT indicating the depth of hydrocarbon occurrence, while in well-3, the GUT and OUT at depths of 8754 and 8934 feet, respectively, show the zones of hydrocarbon accumulation (Fig. 7).

Core data analysis
Representative core data from well-3 were gotten from the interval between 8770 and 9890 feet. Core description shows physical and chemical characteristics of the lithologic intervals Revil et al. 2002;Glover 2010). The analysis was based on the sidewall sample description (Fig. 8). The genetic facies identified were the lower shoreface, upper shoreface, marine shale, tidal channel and channel heterolithics according to Mode et al. (2015) (Fig. 4b). The depositional environments spans through incised canyons, channels, inner mid shelf, shelf margin, and slope margin. The reservoir units represent intervals that are rich in sandstones (89 to 98% of sand). Such intervals contain grey, silty, well sorted and loosely consolidated sands with presence of few heavy minerals, pyrite, rare mica flakes and weathered glauconite (Fig. 8). The sorting and poorly consolidated natures of pore throats of sand grains in the reservoir units imply that these reservoir units have good fluid moveability. The result shows that reservoir intervals evaluated have good porosity and permeability qualities due to sphericity and good sorting of the grains. However, the presence of heavy minerals and clayey particles could reduce porosity and permeability qualities of rocks, because these particles could form pore throat blockages between sand grains thereby reducing pore radius and thus porosity and permeability. The core result shows the occurrence of different lithofacies at depths of 8824, 8830, 9256 and 9259 feet. At certain depths, there are non-reservoir intervals due to protrusion towards the south most position as observed in Well-1 (Fig).

Lithologic identification
Petrophysical and core analysis for depths between 8750 and 9750 feet shows that low gamma-ray log and high gammaray log values were recorded for the sand and shale units, respectively. Gamma-ray log values recorded between 8750 and 8840 feet, 9050 to 9256 feet, 9308 to 9371 feet and 9485 to 9595 feet indicate the presence of reservoir units (Fig. 8). These depth ranges also agrees with the core description for those depths, where very high sand percentages in comparison to shale percentage were observed. The lithologic compositions within these intervals are mainly silts and well-sorted, loosely consolidated sands. On the other hand, between 8790 and 9050 feet, 9250 and 9357 feet and 9601 and 9763 feet, non-reservoir units were identified based on

Fluid type identification
The resistivity log and core description for depths between 8759 and 8840 feet and 9050 to 9256 feet show hydrocarbon reservoir units that have relatively high average resistivity values than adjacent water saturated reservoir units (Tables 5, 7 (Figs. 7 and 8). Water bearing intervals are outlined by low resistivity since resistivity increases when the formation water is replaced by oil while porosity remains constant. Also, density decreases when the water is replaced by oil in the formation with the same porosity, thus the hydrocarbon bearing intervals were not characterised at those points. (Archie, 1942).

Presence of heavy minerals
High gamma-ray responses observed at some clean sandstone intervals were as a result of the presence of heavy minerals such as weathered glauconite, pyrite and mica flakes. For better clarity, the integration of core data with well logs enhanced the determination and differentiation of these radioactive reservoirs from shale. The presence of weathered glauconite, mica flakes and pyrite generally indicates a transitional environment, i.e. coastal/marginal marine environment (Porrenga 1967;Prince 2018). In the studied intervals, glauconite is rare to common, while pyrite, shell fragments, ferruginous materials and mica flakes are very rare. Glauconite and pyrite are important accessory minerals used for environmental studies. Glauconite forms as an authigenic mineral during the early stage of diagenesis of marine sediment, and may also indicate slow rate of deposition (Smulikowski 1954;Hower 1961;McRae 1972;Odin and Matter 1981;Lewis and McConchie 1994;Amorosi 1995) and/or gap in deposition in a case of mature glauconitic mineral with mean potassium (K 2 O) content of ~ 7% (Odin and Fullagar 1988). It is extremely susceptible to subaerial weathering and it is not known as a reworked second cycle detrital mineral (Selley 1997). The presence of glauconite in sand is an indication that the sand is of marine origin. This has also been reported by other researchers within the area from integrated evidences from sedimentology, stratigraphy, mineralogy and geochemical investigation (Lewis 1962;Adegoke et al. 1971 and1976;Petters 1979 and1982;Seiglie et al. 1982;Ojo 1996;Ozumba 1997;Fadiya 1998;Okosun and Liebau 1999;Chafetz and Reid 2000;Goki et al.2007 andOpeloye 2012). More so, pyrites within the shale units suggest reducing conditions. The combination of these accessory minerals with gamma-ray responses helped in sub-dividing the environment into regressive barrier bar, tidal sand wave, and tidal channels Mode et al. (2015). According to Prince (2018), these sub-environments suggest shoreface deposits in marginal marine settings (Fig. 4,  8 and 9).

Electrosequence analysis
Electrosequence involved the interpretation of gamma ray and resistivity log for diagnostic characteristics such as trends, shapes, and abrupt breaks, which are indicative of the depositional environment and energy of deposition (Serra and Sulpice, 1975). The analysis revealed facies ranging from Marine-Shale, Heterolithic channels, upper and lower shorefaces, Shaly-Sandstone and/to Sandstone Facies (Fig. 4b). The depositional systems of Shoreface, prodelta, distributary, tidal and fluvial systems were deduced from gamma ray log signatures. The Bell-shaped log patterns on gamma ray logs indicate a fining upward trends or an upward increase in gamma ray value (typical feature of fluvial channel deposits). The Funnel-shaped log patterns indicate a coarsening upward trend, clearly showed deltaic progradation. Cylindrical (blocky or boxcar) log motif are on thick uniformly graded coarse grained sandstone unit, probably deposits of braided channel, tidal channel deposits (Fig. 4b, 9 and 11).

Sequence stratigraphic interpretation
The sequence stratigraphy analysis provided an opportunity to understand lithofacies, stacking patterns, depositional

Fig. 12
Zonal scheme/bio-zone marker species, and ages, (Modified from Blow 1969) sequences and systems tracts distribution across the study area. The erosional sequence boundaries (SB) as observed are characterised by sizeable sand unit defined in well log by abrupt and sharp bases of thick low gamma ray value intervals usually separating coarsening-upwards intervals from fining upwards intervals ( Fig. 10 and 11). The maximum flooding surfaces (MFSs) are characterised by thick and extensive shale units defined by thick high gamma ray value intervals that separate overall fining-upward interval from coarsening-upwards intervals, and high abundances and diversity of microfossils ( Fig. 10 and 11). These stratal surfaces were delineated using floral and faunal abundance and diversity combined with palynological and foramimiferal zones encountered. Consequently, the stratigraphy of Niger Delta sections has conventionally been based on biochronological interpretations of fossil spores and foraminifers. Palynological and foraminiferal biozones have proved to be reliable tools for large-scale correlations between Hydrocarbon wells (Reijers, 2011). Following the planktonic foraminiferal zonal scheme of Blow (1969) as amended by Kenneth and Srinivasan (1983) (Fig. 12). The biostratigraphic framework and alphanumeric codes based on biostratigraphic zonation customised by the Shell Petroleum Development Company of Nigeria (SPDC) (Fig. 5 and 6)   (Blow, 1969) and F9500-F9600/ F9700 (Shell zonation) in the study area (Tables 1-4). The Shell biozonation was published by Evamy et al. (1978). This zonation allows correlation across all facies types from continental to marine. These biozones are based on index planktonic and benthonic foraminifers recognised from transgressive marine shales (Petters, 1982). The regionally consistent marine marker shales, as well as some of the more local shales (Fig. 13 and Table 10), have usually been named after one of the marker species and are coded systematically (Fig. 5, 6 and 14). The generated paleobathymetric map shows generally that sediments were deposited within neritic through bathyal environments, aligning with the progradational pattern of deposition of the Niger Delta (Fig. 15) (Emery and Myers, 1996;Kendall, 2003). Three MFS were mapped and dated paleontologically; 9.5 Ma (3980 feet), 12.8 Ma (5850 feet) and 15.0 Ma (7662 feet) recognised as maximum shale peaks in while MFS 15.9 Ma (9120 feet) was recognised in Well-2. Also, three SB were mapped and dated; 12.1 Ma (4800 feet), 13.1 Ma (6250 feet) and 15.5 Ma (8450 feet) in well-1 and 15.9 Ma (9120 feet) in Well-2 (Fig. 10, 11 and Tables 1-4). The aforementioned ages indicate Langhian-Serivillian stage of Middle-late Miocene on the global geochronological scheme of Harland et al. (1990) (Fig. 14). Relative thickness of the systems tracts shows changes in rate of sedimentation caused by variations in local conditions like tectonics, effects of gravity, eustatic sea level changes and availability of accommodation space (Haq et al. 1988). Three systems tracts, namely the highstand systems tract (HST), the lowstand systems tract (LST) and the transgressive systems tract (TST) are all represented in Well-1, but in Well-2, only the highstand system tract (HST) and transgressive system tract (TST) were delineated. Sequence stratigraphic correlation across Well-1 and Well-2 shows thickening of sediment package from northeast to southwest ( Figs. 9 and 10). This reveals the evidence of structural influence on stratigraphy which have been attributed to synsedimentary and gravity tectonics within the Niger Delta basin. The evolution of the Niger Delta basin is controlled by pre-and syn-sedimentary tectonic activities described by Evamy et al. (1978), Ejedawe et al. (1984), Knox and Omatsola (1989), Stacher (1995). Mode et al. (2015) agree that the occurrence of the identified sequence stratigraphic surfaces at different depths along dip and strike directions in the 1 3 Fig. 14 Chronostratigraphic data sheet (west and east halves combined) of the Niger Delta. (Doust and Omatsola 1990) wells shows evidence of faulting within the field. Hence, the sediments are thinner in the up-dip (northeast) section and thicken down-dip (southwest) with a corresponding decrease in net-to-gross values (Fig. 10). This thickening of sediment package down-dip could be attributed to the high rate of subsidence and deposition that is influenced by structures on stratigraphy (Evamy et al. 1978). Also, flattening at various MFS reveals a shift of depositional centres from the north to the south. The depositional sequences recognised were made up of genetic units of LST, TST and HST. The sands of LST and HST constitute the reservoir packages, whereas the shales of the TST in which most of the MFS were delineated serves as seal/cap rocks for the reservoir units. The alternation of reservoir units of the LST and HST, and the seal units of the TST offer good stratigraphic configurations for entrapment, hence, was the target for hydrocarbon exploration (Dim 2017). Sequence stratigraphic analysis was not carried out in Well-3 due to non-availability of sequence stratigraphic data for the well. The summary of the identified rise or shale prone wedge with thin fine grained turbidites forms on the slope lowstand prograding wedges. These sediments were deposited when relative sea level was at or near the shelf margin and rising slowly, but were slow enough for sedimentation to take place with shelf edge deltaic process (Posamentier and Vail 1988). The lowstand prograding wedge interval in well-1 is characterised by funnel deltaic sand body that laterally prograded into proximal lagoon shales. The prograding lowstand wedge complex is identified only in well-1 at a depth range of 8350 -8560 feet. This interval corresponds to channel fill environment with progradational parasequence set as shown. (Table 11 and Figs. 9, 10 and 11).

Transgressive systems tract (TST)
Transgressive systems tract develops as a result of increase in the rate of sea-level rise which has an overall deepening upward bathymetric signature. The rate of rise is at a maximum and sediment supply is more than the available space can accommodate. At this stage, the shelf is no longer in tune with the rising sea-level. Deltaic progradation ceases and sands are trapped up-dip in the estuaries (Scholle and Spearing 1998). Transgressive systems tract is observed in Well-1 (7750 to 78,340 feet and 9252 to 9760 feet) and Well-2 (6667-6700 feet) with maximum flooding surfaces (MFS 15 Ma,15 Ma and 15.9 Ma) at the upper boundaries. The constant deepening interval above the lowstand tract in Well-1 at 8350 feet shows a proximal lagoon setting

Lowstand systems tract (LST)
Sediments associated with lowstand systems tract in well-1 occur in incised valleys fill with coarse grained low-sinuosity channel or estuarine sandstone in response to sea level that drowns the funnel lowstand wedge complexes (Kendall 2003). The proximal lagoon shales in Well-1 can serve as seal or cap rocks to the underlying channel fill reservoirs which contain hydrocarbon (Table 11 and Figs. 9, 10 and 11).

Highstand systems tract (HST)
This stage features a general decrease in sea level and it is characterised by initial aggradation of shales into intervals of shallowing-up facies. The rate of sea level rise is at a minimum and depositional rates exceed rate of sea level rise causing parasequences to build basinward in aggradational to progradational parasequence sets. Both fluvial and deltaic sands near the top of the unit prograde laterally into neritic shales. The highstand systems tract was identified in well-1 and 2, respectively, at 6900-7500 feet and 8340-8860 feet for well-1, and between 5960 and 6600 feet for well-2. This prograding highstand complex corresponds to the shallow marine depositional environment. The coarsening upward sequence with frequent deepening well log spikes gives rise to a good to excellent porosity value with occasional shale intercalations that tend to reduce porosity and permeability. Furthermore, the shore sand of the highstand systems tract have good reservoir properties as stated by Coleman and Prior (1980) and Anyiam et al. (2019) (Table 11 and Figs. 9, 10 and 11).

Conclusion
The electrosequence, reservoir character and natural resources potentials of the `Don Field' Central Swamp Depobelt of the Niger Delta were comprehensively assessed and delineated. All reservoirs units in the wells encountered were classified as good reservoirs based on their viable Petrophysical properties delineated through wirelogs and core analysis. The biostratigraphic analysis reveals sequence stratigraphic surfaces marked by MFS 9.5 Ma, 12.8 Ma, 15.0 Ma, 15.9 and SB 12.1 Ma, 13.1 Ma, 15.5 Ma and 15.9 Ma of middle-late Miocene with bio-zones range of N8-N16 or F9500-F9600/F9700 corresponding with fossil index abundance of Globorotalia foshi S.I; Globorotaliapraefohsi; Globorotaliaperipheroacuta; Orbulinaspp; Praeorbulinasicana; Paragloborotalia mayeri; Neogloboquadrina acostaensis. The electrofacies identified are lower shoreface, upper shoreface, marine shale, tidal channel, channel, and channel heterolithics. However, the presence of heavy minerals like pyrite and clay particles as noted in some of the reservoir may affect fluid moveabilty. The core and well log analysis divulges that sparse shale beds were also present, and the sandstone is composed of 89% quartz and about 11% clay and other minerals. The clay is a mixture of illite (60%) and kaolinite (40%). High resistivity values recorded at intervals where sand particles are discoloured by carbonaceous films indicate possibly hydrocarbon zones. The maximum flooding surface mark regional seals and caps the major reservoir units across the various fields in the study area. The LST and HST sand units of the basin floor fans, channel, and shoreface are potential good hydrocarbon reservoirs, whereas the marine shale units of the TST, which is capped by MFS, also provides regional seals to the reservoir units.