Spe 128354
Spe 128354
Spe 128354
Pre - drill pore pressure prediction from 3-D seismic data in parts of the
Onshore Niger Delta basin.
were normally pressured, while nine were over Delta (Akata- Agbada) Petroleum System(see
pressured. Similarly, table 1b shows some of the figure1). The primary source rock is the upper
deepest wells drilled by Shell Petroleum Akata Formation, the marine-shale facies of the
Development Company of Nigeria Limited. All of delta, with possible contribution from
these deep exploration wells drilled by SPDC interbedded marine shale of the lowermost
were overpressured. In most of these cases, the Agbada Formation. Oil is produced from
wells were either abandoned, drilled without sandstone facies within the Agbada Formation,
reaching the objective sequence, or drilling however, turbidite sand in the upper Akata
prolonged unnecessary to sometimes as long as Formation is a potential target in deep water
seven months leading to astronomical rise in offshore and possibly beneath currently
drilling cost. It should be noted that these producing intervals onshore(Turtle,1999).
problems were encountered despite the fact that
the latest drilling practices were applied to most
of the wells (Nwaufa et al,2006).
Safe exploitation of these deep
prospects and plays therefore would largely
depend on our ability to understand the controls
on top seal strength, overpressure generation/
distribution, and trap integrity, and then
confidently incorporating this knowledge into
prospect evaluation and well design. This is very
important considering the current campaign for
deeper wells, as the economic consequences of
exploitation in areas with an unspecified risk of
abnormal pressure profiles may range from
increased drilling costs due to hazards to
unrealized prospect potential. Optimized
planning practices will therefore impact not only Figure1: Regional structural elements of the
on the costs of drilling but also on the quality of Niger Delta basin.
the reservoir evaluation and productivity
assessment.
This work is aimed at predicting pre-drill The Tertiary section of the Niger delta is divided
pore pressures in OML X, Onshore Niger delta into three formations representing prograding
basin using 3-D seismic velocity data. depositional facies that are distinguished mostly
on the basis of sand-shale ratios. Sedimentary
Background geology deposits in the basin have been divided into
The onshore portion of the Niger delta three large-scale lithostratigraphic units: (1) the
province is delineated by the geology of basal Paleocene to Recent pro-delta facies of
southern Nigeria and south-western Cameroon. the Akata Formation, (2) Eocene to Recent,
The northern boundary is the Benin flank - an paralic facies of the Agbada Formation, and (3)
east-northeast trending hinge line south of the Oligocene-Recent, fluvial facies of the Benin
West African basement massif.The northeastern Formation(Evamyetal,1978;ShortandStauble,19
boundary is defined by outcrops of the 67;Whiteman,1982).These formation become
Cretaceous on the Abakaliki High and further progressively younger farther into the basin,
east-southeast by the Calabar flank - a hinge recording the long - term progradation of prog-
line bordering the adjacent Precambrian. The
province covers an area of 300000 km2 and
includes the geologic extent of the Tertiary Niger
3 Pre-Drill Pore Pressure Prediction From 3 – D Seismic Data SPE 128354
Table 1a: shows some of the deepest exploration wells drilled by Nigerian Agip Oil Company
between 1970- 2005 (After Nwaufa et al, 2006)
radation of depositional environments of the strata, changed local depositional slopes, and
Niger delta onto the Atlantic Ocean passive complicated sediment transport paths into the
margin. The stratigraphy of the Niger delta is basin.
complicated by the syn-depositional collapse of
the clastic wedge as shales of the Akata Three major depositional cycles have
Formation mobilized under the load of been identified within Tertiary Niger delta
prograding deltaic Agbada and fluvial Benin deposits (Short and Stauble, 1967; Doust and
Formation deposits. A series of large-scale, Omatsola, 1990). The first two, involving mainly
basinward-dipping listric normal faults formed as marine deposition, began with a middle
underlying shales diapired upward. Blocks down Cretaceous marine incursion and ended in a
dropped across these faults filled with growth major Paleocene marine transgression. The
second of these two cycles, starting in late
4 A.I. Opara and K.M. Onuoha SPE 128354
Table1b: Some deep exploration wells (“HPHT wells”) drilled by SPDC in the Niger Delta
TVD(FT) GRADIENT(PSI/FT)
a) Reservoir continuity
Then tau is related to effective stress by
Reservoir continuity can be assessed by
the relationship.
considering the variation in overpressure. This is
………… (2)
because for a continuous aquifer in a given
where = initial travel time, =vertical
pressure cell, the pressures at different locations
effective stress, = tau, A and B are parameters are expected to lie along a hydrostatic gradient
determined from the behaviour of the sonic log from one another, and are hence equally
data in the normally pressured zone where overpressured. Thus, if a horizon exhibits a
effective stresses are known and, are lithology constant overpressure, this generally indicates a
dependent. hydraulically continuous unit such as a
connected reservoir. In contrast, a horizon which
6 A.I. Opara and K.M. Onuoha SPE 128354
exhibits large changes in overpressure laterally of the hydrocarbon column can be calculated
is expected to be shaly. from the predicted aquifer pressure. This
requires an estimate of the hydrocarbon
b) Maximum hydrocarbon column pressure gradient (i.e. hydrocarbon density), and
the depths of the top structure and spill point for
The height of a possible hydrocarbon column is each reservoir level. The hydrocarbon column
most commonly determined by the depth pressure should leave the predicted aquifer
difference between the top structure and the spill pressure curve and follow the hydrocarbon
point (spill point is the lowest closing contour on gradient up to the spill point. If the pressure at
a hydrocarbon trap capable of holding spill point is less than the minimum horizontal
hydrocarbons under gravitational equilibrium if stress (MHS) trend (taking into account the error
the formation is permeable). Figure3 shows the bars on the determination of the stress trend)
standard procedure for estimating maximum then the trap integrity is expected to be OK. If
hydrocarbon column from pressure depth plot. the hydrocarbon column pressure crosses or
However, in some overpressured areas, the comes very close to (within a couple of hundred
extra pressure at the top of a hydrocarbon psi or 20bar) the MHS trend, there is a very high
column is enough to exceed the fracture trap integrity risk.
strength of the top seal, and cause breaching
and leakage from the trap. Hydrofracturing of the d) Fault Sealing
top seal also has important implications for Analysis of interpreted faults in terms of
supply of hydrocarbons to the shallower pressure is best carried out using vertical
reservoirs in some areas in the Niger delta effective stress (VES) or overpressure plots
basin. Similarly, added pressure from the (where changes in pressure across a fault will
hydrocarbon column may leave a very small be more obvious).It is only possible to suggest
drilling margin. that a fault is sealing if there is a change in
pressure across the fault: a fault in a
hydrostatically pressured area will exhibit the
same pressure on either side of the fault,
whether or not it is sealing. In interpreting
sealing faults, it is important to bear in mind the
sealing mechanism. If the fault seals by
juxtaposition of sands against shales, then the
pressure difference should be correct. But a
cataclastic or clay smear fault seal can only
withstand a certain pressure difference – thus
not all faults which are associated with a change
in the pressure predicted from seismic are
necessarily sealing.
e) Pressure prediction
Pore pressure prediction is important
Figure2: Diagram for estimation of the for both the exploration and exploitation of
maximum hydrocarbon column from hydrocarbons. In exploration for hydrocarbons,
pressure-depth plot (after Indrelid, 1997). knowledge of the pressure distribution is of vital
importance for the development of fluid
c) Trap Integrity migration models, to study the effectiveness of
To assess the risk of trap breaching and seals and also to rank prospects. Similarly,
hydrocarbon leakage, and also the likely during drilling, pressure prediction allows the
pressure on drilling into the reservoir, the effect planning of the well in real time to control
7 Pre-Drill Pore Pressure Prediction From 3 – D Seismic Data SPE 128354
The 3-D seismic velocity data were across a fault indicates a higher risk of fault seal
picked down to 6seconds at an average of 250 breakdown. Note that there are no visible
millisecond interval on a grid of about 24m x changes across the faults in the normally
24m.The seismic section revealed dense pressured section. Similarly, In figure 9d, the
faulting and structural deformation. Structural predicted MES at prospect A and B have the
interpretation of the seismic data of the area same value of 9523- 11428psi while prospect C
revealed that the structural style is dominated by has an MES of 3809– 5238psi (27.5 – 37.4Mpa).
growth faulted rollover anticlines with footwall This shows that prospects A and B have a much
and hanging wall closures and collapsed crest lower risk of trap breaching than prospect C.
faults, associated with shale diapirs. The root Similarly, this risk in prospect C is increased
mean square (RMS) velocity data was converted even more by considering the effects of a
to interval velocity using appropriate (medium) potential hydrocarbon column. The structures of
smoothing parameters as shown in figure 6 prospects A and B are low and hence a
below. Note that the interval velocities are not relatively short hydrocarbon column (hence
erratic but some significant velocity variations small hydrocarbon buoyancy pressure) is
can be observed. There is a major reversal in expected in contrast to prospect C which has a
velocity associated with top of overpressures higher relief. Filling the prospect with
resulting from the influence of shale diapirs hydrocarbons would result in an even smaller
below 2.5seconds. The analyzed seismic minimum effective stress, and in all likelihood a
velocity data which has been depth converted blown trap. Hence, prospect B is ranked highest
were extracted at proposed well locations and followed by A, while prospect C has the least
compared with the checkshot and sonic ranking.
velocities of offset wells for the purposes of Fluid pressures predicted from the
identifying the presence of anisotropy and bad seismic velocity cube and extracted at well
seismic data(Figure 7).However, it appears that locations are shown in figure 10 below. In Ah-
the data quality is good enough since no field (figure 10a), +500psi above hydrostatic
considerable anisotropy effect was noticed. The pressure (overpressure) is at about 6000ftss and
calibration parameters (A&B which relates the is taken to be the top of overpressure in the
vertical effective stress at any depth to the area. The top of overpressure can be seen to
porosity and density at the same depth) together vary from 8000ftss in Oo-field, to 9000ftss at Og-
with the seismic velocities corrected for field and then 1100ftss at Is-field.
anisotropy were used to predict the pore Finally, maximum hydrocarbon column
pressure field of the area. Figures 8 shows the estimated from the predicted pressure depth
pore pressure cubes of the study area plots in figure 10 vary between 35ft to 350ft with
generated from seismic data giving a 3-D an average column of 55ft. This relatively short
visualization of the pressure field. hydrocarbon column is associated with the fact
The predicted pressure cubes were then that most of the reservoirs in the study area are
interpreted as a backdrop to seismic data as not filled to their synclinal spill points.
shown in figure 9 below. The purpose of this is
to predict the effects of trap integrity, reservoir Interpretation and discussion
continuity, sealing faults and maximum Pore pressure prediction and
hydrocarbon column. Figures 9a-9d all exhibited interpretation in OML X, Onshore Niger delta
a very wide variation in predicted pressure basin revealed that the onset of mild
below 2.5seconds which indicates absence of overpressure(<0.71psi/ft) in the area to be about
connected reservoirs in the overpressured 6000ftss in Oo -field, 8000ftss around Ah -field,
sections. Similarly, figure 9c showed distinct 9000ftss in Og-field and about 11000ftss at Is -
changes across the faults in overpressured field. Analysis of seismic data revealed that very
section. This shows that the faults are sealing. A high (hard overpressures nearing lithostatic
bigger difference in the predicted pressure
9 Pre-Drill Pore Pressure Prediction From 3 – D Seismic Data SPE 128354
Figure 7: Seismic velocity data at well locations calibrated with checkshot and sonic velocity data
Figure 8a:3 –D Pore Pressure Cube Generated Figure 8b:3 -D Vertical Effective Stress Cube
Using Conventional Velocities Generated From Seismic Data Using DMO
Velocity
Figure 8: Pore pressure cubes from seismic data generated using velocity to pressure transform.
10 A.I. Opara and K.M. Onuoha SPE 128354
Figure 9: Seismic generated pressure cubes interpreted as back-drop to seismic data for
interpreting reservoir continuity, maximum hydrocarbon column, sealing faults and trap integrity
Figure10: Seismic predicted fluid pressure- depth plots in the study area.
11 Pre-Drill Pore Pressure Prediction From 3 – D Seismic Data SPE 128354
Distribution of overpressures shows a well-defined the Akata Formation. However, previous studies
trend with depth to top of overpressures increasing by Ichara and Avbovbo(1985); Obah(1989) and
towards the central part of the basin at a maximum Weber and Daukoru,1975,have shown that
depth of about 13000ft (3940m). This variation in overpressures in the Niger delta often occur at
the depth of top of overpressures within the area is depths shallower than Akata Formation.
believed to be related to faulting and shale Explaining this phenomenon, Weber and Daukoru
diapirism with top of overpressure becoming (1975), revealed that overpressures are
shallower with shale diapirism and deep with encountered in the Tertiary Niger delta as a result
sedimentation. Similarly, data acquired in deep of rapid loading of the undercompacted shales of
wells from across the concessions in the Niger the Akata Formation by the sandy Agbada and
delta show no macro-structural trend.Rather, Benin Formations. The Akata Formation is in
overpressures were observed in areas with simple contact with the sandy paralic Agbada sediment in
rollovers, especially at the hanging walls, while three different ways. In the first place, there is the
normal formation pressures were observed in vertical transition from continuous marine shale
areas with k-faulting pattern and collapsed crest into paralic sediments. Secondly, there are lateral
structures. It is most likely that these fault patterns facies transitions and interfingering of sand and
have effect on the formation pressure, providing clay and thirdly Akata shale is in many places in
relief to a potential pressure build-up (Nwaufa, et juxtaposition with Agbada paralic sediments
al, 2006).Hanging wall fault closures are most across faults. In each of these cases, fluids
often sealing and therefore retains significant expelled from the overpressured Akata shales
columns of hydrocarbons often leading to high may inflate (charge) the pressures in the adjacent
minimum effective stress and fluid charging. sands. Similarly, most of the mild overpressures
Similarly, late hydrocarbon generation and shale within the Agbada Formation in the Niger delta are
diagenesis often leads to overpressures within as the result of undercompaction (possibly
upthrown shales(Evamy,etal,1978).The chemical compaction disequilibrium) of the
subsurface of the Niger delta basin is extensively interbedded marine shales of the lowermost
deformed by growth fault structures and roll over Agbada Formation. Consequently, overpressures
anticlines (Weber and Daukoru, 1975; Evamy etal, are often encountered before the Akata shale is
1978; Merki, 1972 Xiao and Suppe, 1992). reached.
Hanging wall rollover anticlines developed
because of listric-fault geometry and differential Conclusion
loading of deltaic sediments above ductile shales. Finally, overpressures observed in rocks
These growth faults are characterized by thicker within the Niger Delta owe its distribution not only
deposits in the downthrown (hanging wall) block to the mechanisms but also to the re-distribution of
relative to the upthrown block. Here, fluids during and after the generation of
sedimentation is rapid on the hanging wall with overpressures, as a result of structure and
respect to the footwall leading to a thicker stratigraphy. Changes in formation relief, geometry
sediment pile at the hanging wall. This rapid and faulting play major roles in pressure
sedimentation often leads to overpressuring distribution and re-distribution within the basin.
effects. Overpressure in the Niger delta is not caused by
In addition, it was revealed that majority of mechanisms associated with undercompaction
the overpressure situations occur between the alone; there are other dominant mechanisms
depth intervals of 6000ft to 13000ft(1820 and especially thermal expansion mechanism (where
3940m) falling largely within the Agbada fluids charging resulting from aquathermal,
Formation. In the Niger delta basin, the hydrocarbon generation and clay diagenesis all
occurrence of overpressures is largely believed to working together) occurs. These mechanisms
be associated with the undercompacted shales of seem to correlate with high pressures at great
12 A.I. Opara and K.M. Onuoha SPE 128354