Destiny Project
Destiny Project
Introduction
The Akata Formation, nestled within the extensive Niger Delta Basin in West
significance. This section furnishes the essential backdrop for the research,
The Niger Delta Basin, also known as the Niger Delta province, is located in the
Niger Delta and the Gulf of Guinea along the western coast of Nigeria. This
extensional rift basin is complex and holds significant economic value due to its
highly productive petroleum system. It ranks among the largest subaerial basins in
Africa, covering an extensive area. The basin's sediment fill is voluminous and
provide insights into its formation and regional tectonics. The Niger Delta Basin is
part of an extensional basin network within the broader tectonic structure known as
the Benue Trough, with the Cameroon Volcanic Line and a passive continental
margin bordering it.The formation of the Niger Delta Basin can be attributed to the
unsuccessful separation of the South American plate and the African plate, which
occurred during the initial stages of the South Atlantic opening. This rifting process
commenced in the late Jurassic period and concluded in the middle of the
with many of them taking on the form of thrust faults. Concurrently, sediment
deposition occurred in the late Cretaceous, initially in the form of syn-rift sands
leading to the formation of high-angle normal faults and the rotation of fault
blocks. Transitioning into the Paleocene era, a notable shoreline transgression took
place.
During the Paleocene period, the Akata Formation was deposited, followed by the
Agbada Formation in the Eocene. This deposition process caused the underlying
Subsequently, in the Oligocene, the Benin Formation began to accumulate, and this
deposition process is still ongoing today. The entire basin can be categorized into
distinct zones based on its tectonic structure. These zones consist of an extensional
zone, located on the continental shelf and resulting from crustal thickening, a
transition zone, and a contraction zone situated in the deep-sea region of the basin.
Fig 1.1 Geologic map of the Niger Delta Basin and the Benue trough, and the oil fields in the
region.
1.2 Exploration Significance:
hydrocarbon reservoirs in the Niger Delta Basin. Exploration efforts often target
areas where the Akata Formation acts as both a source and a seal for hydrocarbons
primary source rock. Its organic-rich shale layers have generated substantial
hydrocarbon reserves over geological time, making it a focal point for assessment
1.2.3 Structural Analysis: The subsidence and diapirism associated with the
a crucial role in both the economic development of the region and its
and natural gas. Its clay-rich sediments have acted as both a source rock and a seal,
generating and trapping vast quantities of hydrocarbons beneath the Niger Delta
(Tuttle et al., 2015). This has made the Niger Delta one of the world's major oil-
the Niger Delta, including the Akata Formation, have raised environmental
concerns. Spills and pollution associated with oil and gas activities can harm local
ecosystems, aquatic life, and vegetation. These impacts have been a source of
1.2.7 Wetland Preservation: The Niger Delta is home to rich and diverse
including those related to the Akata Formation, are essential for safeguarding these
Given these challenges, there is a clear need for dedicated research efforts to
Formation. This research should particularly focus on AVO, OWC, and brightspot
analysis, aiming to bridge the knowledge gap and enhance the precision of
1.What are the seismic characteristics and attributes of the Akata Formation in the
2. How can AVO analysis be applied to predict hydrocarbon reservoirs within the
Akata Formation?
4. What are the key indicators and methodologies for identifying brightspots within
The objectives of this study are the guiding beacons that illuminate our path toward
1. To investigate and analyze the seismic characteristics and attributes of the Akata
3. To develop methods for determining the Oil-Water Contact (OWC) within the
Akata Formation.
The significance of this study extends far beyond the boundaries of academic
Akata Formation's seismic attributes, and in doing so, it holds the promise of
and foremost, this study contributes to the reservoir of geological knowledge. The
Akata Formation, with its enigmatic alternating shale and sandstone layers, has
structures and geological processes that have shaped the formation over millennia.
Yet, the implications of this study extend well beyond the academic sphere. The
Basin and beyond. By applying advanced Amplitude Versus Offset (AVO) analysis
techniques, the study seeks to provide explorers with a refined toolkit for
predicting the presence of hydrocarbon reservoirs. This not only enhances the
precision of exploration strategies but also mitigates risks associated with drilling
determining the Oil-Water Contact (OWC) within the Akata Formation carries
immense practical weight. Accurate OWC determination is a linchpin of effective
formation, implications for hydrocarbon exploration AVO, OWC, and Brights pot.
Chapter 2
role. The Akata Formation, situated in various offshore basins worldwide, has
geological setting for exploration, marked by its distinct seismic characteristics and
Contact (OWC), and Bright Spots have emerged as critical tools for deciphering
formations like the Akata. These attributes are derived from seismic data, offering
valuable insights into the subsurface, and have become integral to the toolkit of
This literature review aims to explore the seismic characteristics and attributes of
the Akata Formation within the context of hydrocarbon exploration. It delves into
underpin our understanding of AVO, OWC, and Bright Spots in this geological
setting. By synthesizing existing knowledge and research, this review seeks to shed
light on the relevance and implications of these seismic attributes for optimizing
This introduction sets the stage for the theoretical framework that will guide our
application of these attributes to formations like the Akata Formation stems from
their ability to transform seismic wave reflections and refractions into valuable
extracting or deriving a quantity from seismic data that can be analysed in order to
delineation of bright spots, flat spots, polarity reversals and dim spots, coherency,
and amplitude versus offset. Attributes that can show the presence of hydrocarbons
accumulation in the subsurface. Types of DHI are bright spots, flat spots, ‘dim
Based on the this study we will be explaining the following seismic attributes
and geophysics. It involves studying how the amplitude of reflected seismic waves
changes concerning variations in the offset, which refers to the distance between
the seismic source and the receiver. This technique is employed to glean valuable
The foundation of AVO analysis lies in the fact that different rock layers exhibit
reservoir and surrounding geological formations, part of the incident wave energy
gets reflected back to the surface while the remainder continues to propagate
downwards.
This dependency of reflected seismic wave amplitude on the offset angle at which
it's recorded has been widely recognized since the early 20th century when Karl
describe the relationship between the reflection coefficient and the angle of
incidence. AVO analysis primarily focuses on how variations in the offset angle
Fig 2.1. Diagram showing how the layout of sources and receivers affects the angle
of incidence.
The basic concept behind the AVO technique is the analysis of the reflection
process. we take advantage of the fact that the measured amplitude is related to the
strength of the
reflection (reflection coefficient), and that the reflection coefficient depends on
three
parameters:
1. Change in P-wave velocity across the interface
2. Change in S-wave velocity across the interface
3. Change in density across the interface
Vp1, Vs1 , ρ1 - P-wave velocity, S-wave velocity, and density for layer 1.
Vp2, Vs2, ρ2 - P-wave velocity, S-wave velocity, and density for layer 2.
At normal incidence the reflection coefficients depend only on and
The Akata Formation, part of the sedimentary fill in the Niger Delta Basin,
al., 2017). This formation is situated above the older syn-rift sediments of the late
Cretaceous period and is overlain by the younger Agbada Formation, which was
The formation's clay-rich nature reflects the fine-grained nature of the sediments,
which were derived from the erosion and weathering of pre-existing rocks in the
The Akata Formation is a prolific hydrocarbon source rock in the Niger Delta.
Organic-rich shale layers within the formation have undergone thermal maturation,
natural gas, migrate from the Akata Formation to reservoir rocks within the basin
(Tuttle et al., 2015). This makes the Akata Formation a critical component of the
primary source rock and a potential target in deep water offshore and possibly
been generated or are capable of being generated. They form one of the necessary
elements of a working hydrocarbon system. They are organic rich sediments that
Growth faults are the dominant structural features in the Niger delta which are
load and gravitational instability of the Agbada sediment pile accumulating on the
mobile undercompacted Akata shales. Toe thrusting at the delta front, lateral flow
and extrusion of the Akata pro- delta shales during growth faulting and related
extension also account for the diapiric structures on the continental slope of the
Niger delta in front of the prograding depocenter with paralic sediments (Doust
and Omatsola, 1990; Reijers, 1996). Three major depositional cycles have been
identified within Tertiary Niger Delta deposits (Short and Stauble, 1967; Doust and
Omatsola, 1990). The second of these three cycles, starting in late Paleocene to
Eocene time, reflects the progradation of a “true” delta, with an arcuate, wave and
tide dominated coastline. These sediments range in age from Eocene in the north to
Quaternary in the south (Doust and Omatsola, 1990). Deposits of the last
depositional cycle have been divided into a series of five depobelts (Doust and
and most post sedimentary faulting ceased with the abandoned depobelt. A
depobelt therefore, forms the structurally and depositionally most active portion of
Research Methodology
The seismic reflection technique measures the time taken for a seismic wave to
travel from a source (at a known location at or near the surface) down into the
ground where it is reflected back to the surface and then detected at a receiver,
which is also at or near the surface at a known position as shown in Figure 2.5.
This time is known as the two way travel time (TWT). The seismic method gives
materials present in the subsurface. The seismic method is divided into three
processes:
Data acquisition
Data processing
finding the depths of reflecting surfaces and the seismic velocities of subsurface
rock layers. The Principle involves, a seismic signal e.g. an explosion is produced
at a known place and time and the echoes reflected from the boundaries between
rock layers with different seismic velocities and densities are recorded and
analyzed
3.2 Well Data
Six wells were made available in the data set 4 of the wells are vertical (KB-1, KB-
2, KB-4 and KB-5) while 2 are deviated (KB-3 and KB-6). The wells are presented
Data Loading/ Troubleshooting Data QA/QC Well data evaluation(Pick tops and
correlate across wells) Well-to-seismic tie Fault mapping Horizon mapping Time
importing the seismic data which was in SGY format, then the wells were loaded
one after the other. The well data were in LAS File. Directional data for the
deviated wells (ASC File) were loaded thereafter. Check shot data for the various
wells were loaded then loaded. The check shots format was in Text document.
After data loading, data quality check was carried out. The assessment of the data
was to know the quality of the data and to have knowledge of the availability of
some data such as the well logs. This step gives the interpreter an idea of how to go
Reservoirs were identified and correlated across the six wells using the Gamma ray
logs and Resistivity logs as seen in Figure 3.3. The lithologies penetrated by the
wells were identified using the Gamma ray log. A shale base line was established.
Deflection of the Gamma ray log signature to the right of the shale base line was
interpreted as shale (non-reservoir lithology) and deflection of the Gamma ray log
signature to the left side of the shale base line was interpreted as sand (reservoir
lithology). For the resistivity log, deflections to the left were interpreted as low
resistivity (high conductivity) and deflections to the right were identified as high
resistivity (low conductivity). Hydrocarbon bearing reservoirs are characterized by
net reservoir thickness, gross reservoir thickness, porosity, water saturation and
hydrocarbon saturations were calculated. A composite line was drawn across the
base map which assisted in the order in which the wells were arranged.
The stacking patterns of the lithology signature from the logs play a vital role in
predicting the environment in which sediments were deposited. In this study, due
deposition was based on the shape of the gamma ray log only. The principal shapes
observed were the bell, the funnel and the cylinder. The gamma ray log was used
because of its ability to give greater variety of shapes, greater definition and has
more „character‟ than other logs. Some of the environments investigated for in this
mostly associated with a bell-shaped gamma ray log signature which indicates a
associated with a serrated funnel-shaped gamma ray log signature which indicates
shaped gamma ray log signature which indicates a transgressive marine shelf and a
funnel shaped gamma ray log signature which indicates a prograding marine shelf.
Deep Marine Environments. The deep marine environments includes the slope
channel, the inner fan channel, the middle fan channel, the supra-fan channel and
3.3.5 Petrophysics
In this study, petrophysics will be carried out to have a better understanding of the
reservoir, interconnection of the pore spaces and how they affect the migration and
lithology, porosity, water saturation, density and permeability. Well logs are often
Well-to-seismic tie is the process of tying your well to your seismic data with
perform the well-to- 40 seismic tie, a synthetic seismogram was generated and
matched to a real seismic trace thereafter features from the well were then
correlated to the seismic data. The procedures for this process are:
1. Editing and calibrating the sonic and density logs. This process is simply referred
choosing the appropriate reflection series (usually p-waves) and constructing this
series in two-way-time.
3. Performing the match. In doing this, the best match location is determined.
4. After the well-to-seismic tie is done, the quality of the tie can be adjusted by
· Bulk shift
· Squeeze
· Stretch
The sonic log calibration is the process of bringing the sonic log into agreement
with the seismic times. Practical and theoretical analysis of the factors that
influence the accuracy of well ties shows that timing is paramount (White 1997).
calibration started with picking checkshot times. Thereafter knee points were
chosen. Knee points have to be chosen at major jumps in the sonic log so as to
avoid introducing artificial or abnormal reflection coefficients. Note that
character of the logs and knee points allow direct control of the calibration. Drift
At the start of every seismic interpretation, there must be an attempt to tie seismic
seismic data interpretation is to ensure that the well seismic and the surface seismic
at the well trajectory look as similar as possible. After this is done, events in the
well log and on the surface seismic can be linked thereby subsequently correlating
in the Earth‟s crust are as a result of plate tectonic forces. In this study, it was
essential to map out fault networks on the seismic data. Faults were mapped along
the in-lines before been traced out on the cross-lines. The faults were mapped at a
line spacing of 25m along the in-lines. Faults were identified on the seismic
sections by observing the following seismic behaviors which serve as an indication
of faults:
The horizon mapping was done after the faults were mapped. Before mapping
horizons, seismicto-well tie must be carried out using the generated synthetic well
After the mapping of the various horizons, the time map for each horizon was
generated. The time map is characterized by contour lines which show points of
equal elevation. Points of equal elevation are also characterized by the same
colour.
3.4.2 Time-Depth Conversion
The time-depth conversion was carried out using Microsoft Excel 2013. The
Checkshot data was used. A Velocity model was plotted with depth (ft) against time
equation for the trendline was generated which was used in the time map to depth
The depth map was generated using the 2nd order polynomial which was generated
in the time depth conversion. The depth map was used in prospects identification
and evaluation. The depth map is characterized by contour lines which show points
of equal elevations.
Isopach maps were generated to show the thickness between each of the horizons
that were interpreted. The Isopach maps were characterized by contour lines
joining areas of equal elevation and colours which show points of equal elevations.
phase, frequency, polarity and velocity are useful for hydrocarbon indication.
spots and flat spots. Seismic attributes are useful in geological interpretations for
temperature and polarity of seismic. Sheriff (1980), Chambers and Yarus (2002).
The seismic attributes analyzed in this study includes RMS Amplitude, sum of
values and deciding its configuration, whether a four way closure, fault assisted or
attributes.
using the maps and the petrophysical parameters of the particular reservoir.
Volumetric estimated was carried out to find the Oil originally in place (OOIP),
original gas in place (OGIP) and Stock tank oil initially in place (STOIIP).
Chapter 4
In this chapter, the reservoirs observed in the various wells, their depositional
The geological and geophysical properties of the field are presented including the
various maps and models. All the maps are presented, tested and prospects are
A lithostratigraphic correlation was carried out across all the wells which was done
by mapping out the top and base of each reservoir unit. 8 reservoirs units were
correlated which includes Sand A, Sand B, Sand C, Sand D, Sand E, Sand F, Sand
composite line from west to east linking close wells together as seen. The
correlated stratigraphy shows almost uniform layer thicknesses across the field for
most of the levels. However, sand F indicates marked variation in thicknesses with
the thickness part centred around KB-5 (well 5). The thickest part of sand F is
erosion of layer.
4.1.2 Depositional Environments
The depositional environments for all the sands were investigated and the results
Sand A: Sand A has an average thickness of about 400ft across all wells. The
gamma ray signature in figure shows a funnel shape succession. The funnel shape
Sand B: Sand B has an average thickness about 200ft across all wells. The gamma
Sand C: Sand C has an average thickness of about 290ft across all wells. The
aggradational stacking pattern. This stacking pattern is most likely associated with
Sand D: Sand D has an average thickness of about 260ft across all wells. The
gamma ray log indicates a blocky slightly serrated profile. This profile indicates an
aggradational stacking pattern which is associated with a channel-fill. The top of
Sand E: Sand E has an average thickness of about 300ft across all wells. The
gamma ray log indicates a funnel-shaped succession. The funnel shape indicates a
Sand F: Sand F has an average thickness of about 1000ft across all wells. The
gamma ray log indicates a blocky, highly serrated profile. This indicates an
aggradational stacking pattern. This stacking pattern is most likely associated with
a channel-fill.
Sand G: Sand G has an average thickness of about 400ft across all wells. There is
a discontinuity of Sand G in well TMB-05. This can be due to faulting. The gamma
ray response shows a funnel-shaped succession which is very serrated. The funnel-
environment.
Sand H: Sand H has an average thickness of about 190ft across all wells. The
gamma ray log indicates a funnel- shaped succession which is highly serrated. The
The faults were mapped on every 16 seismic section. Faults were noticeable mostly
on the in-line because they were acquired in a direction that is perpendicular to the
regional fault trend. Figure 4.3(a) shows the un-interpreted seismic section while
Figures 4.3(b), 4.3(c), 4.3(d) and 4.3(e) show the interpreted faults. 10 faults
denoted as F1, F2, F3, F4, F5, F6, F7, F8, F9 and F10 were interpreted from the
seismic data (Table 4). Faults F1 and F4 are major faults extending across the
whole field known as major regional growth faults, and they are interpreted to be
the major faults that bounds the coastal swamp depo-belts from the offshore and
central swamp depo-belts (Doust and Omatsola, 1990). Synthetic and antithetic
Six Horizons were interpreted across the field which includes sand A, sand B,
sand E, sand F, sand G and sand H as shown in Figure 4.4. Time maps, depths
maps and seismic attributes maps such as RMS amplitude, sum of amplitudes, sum
of energies and average instantaneous phase were generated for each of these
horizons.
4.2 Time Maps
The structural time map for sand A as shown in Figure 4.5(a) is characterized with
a contour interval of 5ft. The colors show difference in elevation of the horizon.
Four faults (F1, F3, F4 and F6) cut through this horizon.
a contour interval of 5ft. The colors show difference in elevation of the horizon.
Four faults (F1, F3, F4 and F6) cut through this horizon.
The structural time map for sand E as shown in Figure 4.5(c) is characterized with
a contour interval of 10ft. The colors show difference in elevation of the horizon.
Four faults (F1, F3, F4 and F6) cut through this horizon.
The structural time map for sand F as shown in Figure 4.5(d) is characterized with
a contour interval of 10ft. The colors show difference in elevation of the horizon.
Four faults (F1, F3, F4 and F6) cut through this horizon.
The structural time map for sand G as shown in Figure 4.5(e) is characterized with
a contour interval of 10ft. The colors show difference in elevation of the horizon.
Four faults (F1, F3, F4 and F6) cut through this horizon.
The structural time map for sand H as shown in Figure 4.5(f) is characterized with
a contour interval of 10ft. The colors show difference in elevation of the horizon.
Four faults (F1, F3, F4 and F6) cut through this horizon.
4.2.7 Time-Depth Conversion
Time-depth plot of checkshot data is shown in Figure 4.6. The second order
polynomial “y=0.0005 + 2.6331 + 154.07” was generated from the time-depth plot.
It was used for the time-depth conversion, where y=depth from the well data and
The structural depth map for sand A as shown in Figure 4.7(a) is characterized
with a contour interval of 20ft. The colors show difference in elevation of the
horizon. Four faults (F1, F3, F4 and F6) cut through this horizon. A drilled and
tested area is identified on the map alongside two undrilled and prospect areas. The
tested area is seen to be a fault assisted closure being assisted by fault F3. The
undrilled and prospect areas are both 4-way closures. Figure 4.7(b) shows the sand
5.1 Conclusion
3D seismic interpretation was carried out on the KB Field using well logs and
Seismic data volume. The six wells used in this study were correlated and eight
Ten faults were interpreted from the seismic inlines denoted F1, F2, F3, F4, F5, F6,
F7, F8, F9 and F10. All interpreted faults are normal faults – growth faults (F1 and
F4), synthetic faults (F2, F6, F7 and F9) and antithetic faults (F3, F5, F8 and F10).
Within the study area, there are neither thrust nor strike-slip faults. North dipping
hydrocarbons in wells KB-1, KB-3, KB4 and KB-6. Structural highs such as
rollover anticlines were interpreted to be part of the geology of the study area. Five
selected horizons were mapped out which include; H1, H2, H3, H4 and H5. These
horizons were used for the location of the hydrocarbon zones in the field. From the
structural maps, it was observed that the hydrocarbon accumulations are associated
with anticlines, fault assisted closures and fault dependent closures. The tested
areas from all sands were supported by the seismic attributes maps that were
analyzed. Only prospects in sand A, sand G and sand H were supported by the
seismic attributes maps. The petrophysical values which include the porosity, Net
to Gross, Water saturation, hydrocarbon saturation have values that are almost ideal
for the Niger Delta reservoir sands with average porosity value of 0.24, water
saturation of 0.52 and average Net to gross value of 0.6. The lower the water
saturation, the higher the hydrocarbon saturation in the reservoir sand. 113 From
the volumetric estimations, sand G is seen to have the largest hydrocarbon zone
5.2 Recommendations
I recommend that the prospects should be tested to improve the viability of the
assessment
References
Abam T.K.S (2016). Engineering Geology of the Niger Delta, Journal of Earth
Adegoke O.S., Oyebamiji A.S., Edet J.J., Osterloff P.L., Ulu O.K., (2017).
Ayolabi E.A and Adigun A.O., (2013). The use of seismic attributes to enhance
research, 2(2).
Shaw
analysis of Z Field, eastern Persian Gulf. Article in the Leading Edge 2017.
Doust, H. and Omatsola, E. (1990) Niger-Delta. In: Edwards, J.D. and Santogrossi,
Three dimensional seismic, well log and structural analysis of Igbobi Field,
Edwin A.N.,Agyingyi C.M., Nton M.E. and Oladunjoye M.A., (2018). Seismic
Ejedavwe J., Fatumbi A., Ladipo K. and Stone K., (2002). Pan-Nigeria exploration
well lookback (post-drill well analysis). Shell Petroleum Development
Etuk N.O., Aka M.U., Agbasi O.E. and Ibuot J.C., (2020). Evaluation of seismic
attributes for reservoir characterization over Edi field, Niger delta, Nigeria
168-172.
Fajana A.O., Ayuk M.A., Enikanselu P.A. and Oyebamiji A.R., (2018). Seismic
C., Lutgens F.K., Tarbuck E.J. and Tasa D., (2012). Essentials of geology,
11th edition, Boston: Prentice Hall, page 32. Marfurt, K., (2005). 3D
Characterization.
Okeke H.C., Okoro A.U., Ezediegwu P.C., Chinwuko A.I., Onuigbo E.N. and
Vol.8, No.2 pp. 65- 78. Olatunbosun L.G. (2015). 3D Seismic Reflection
3(6): 178-184.
Opara A.I., Anyiam U.O. and Nduka A.V.,(2011). 3D Seismic Interpretation and
Satti I.A., Deva Ghosh and Wan Yusoff W.I., (2014). Integrated Seismic