Comparison of the Characteristics of Low Velocity Layer (LVL) in the Mangrove Swamp and in the Upper Flood Plain Environments in the Niger Delta, using Seismic Refraction Methods

Copyright: © 2016 Uko ET, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.


Introduction
Seismic exploration involves the generation, detection, analysis and interpretation of elastic waves in the earth to study the sub-surface properties of the earth. A seismic section or profile should represent accurately the configuration of the subsurface. But due to topography and near-surface irregularities and other factors, this is not so. One of the factors is the weathered zone, otherwise called the low velocity layer (LVL) or weathered layer. This portion of Earth affects and is impacted by various human activities such as building, excavating, tunnelling, storing or accidentally releasing hazardous materials. Moreover many of the planet's minerals, groundwater, and cultural resources are located in the uppermost subsurface layers [1][2][3][4][5][6][7][8][9][10].
This layer effectively slows down the propagation of seismic waves through the earth and hence increases their travel times. This can lead to erroneous depth predictions of structures during final interpretation of seismic sections and reflections may appear deeper than they actually are [11,12]. In seismic reflection work, one very important method in eliminating this problem is the application of static corrections. But before static corrections can be determined and applied, the LVL characteristics of velocity, thickness of the LVL and its elevation must be known and thus incorporated into the processing programme of seismic reflection data. The aims of this research are, therefore, to determine the LVL characteristics of thickness, weathering velocities and the consolidated layer velocities in the Upper Flood Plain and the Mangrove Swamp areas in the Niger Delta, and to make a comparison between these characteristics in the two zones. The study is of high significance in the sense that the computed results from the LVL can be used in static, elevation and dynamic corrections in reflection seismic survey [4,13] in water exploration [2]. The LVL characteristics can also be used in foundation studies for the design and construction of roads, bridges, high-rise buildings and dams [5].

Theoretical background
The weathered layer or low-velocity-layer: The weathered layer Abstract Sixteen (16) surface-laid-geophones and fourteen (14) downhole-laid-hydrophones experiments were conducted in Upper Flood Plain and in Mangrove Swamp areas respectively in parts of the Niger Delta. The aim is to compare low-velocity-layer (LVL) characteristics in the two environments. The velocity and depth of the weathered layer and those of the consolidated layers were calculated using depth-time plots. The interpreted data showed a substantial variation of the weathered layer thickness and elevation in the two study areas. In the Upper Flood Plain, LVL thickness varies between 2.8 m and 40 m with an average of 21.46 m. In the Mangrove Swamp, the thickness varies between 2.0 m and 5.5 m with an average of 3.40 m. The weathered layer thickness in the Mangrove Swamp is fairly uniform having an average of 3.4 m. This highly variable LVL thickness indicates the necessity of correcting for this layer during seismic reflection processing. of the layer). Shots taken in this layer tend to be of low frequencies because the layer is capable of absorbing high frequency signals and releasing lower frequency ones. Since higher frequency signals contain more information on the subsurface, it is appropriate that in order to acquire good quality reflection data, shots have to be taken below the weathering layer. Because of the first factor, records from shots in this layer often are of poor quality and efforts are made to locate the shot below the LVL.
Surface-laid-detectors refraction method: Surface-laid-detectors refraction method is a technique to map near-surface geologic structure by using head waves that are detected by surface-planted geophone detectors. The source/shot point is close to the surface. Head waves from the source enter a high-velocity medium refractor at critical angle and travel in the high-velocity medium nearly parallel to the refractor surface and are detected by the surface-laid geophones. The objective is to determine the arrival times of the head waves to the refractors [7,8]. The velocity of seismic waves within overburden and rock layers is controlled by the important parameters of elasticity and density. The ground is generally assumed homogenous in its elastic properties, and its velocity will therefore vary both vertically (with depth) and laterally.
For a two-layer medium in which the energy is located at a depth of D s in the first weathered layer and sub weathered layer [D 1 ] (Figure 1), the thickness of the first weathered layer (D o ) is given by Dobrin, Uko et al., [2,15]: is the topmost layer of the earth's surface. It should not be confused with geologic weathering although it is often the result of it. For weathered or low velocity layer, seismic velocities which are lower than the velocity in the underlying consolidated layer imply that gas (air or methane resulting from the decomposition of vegetation) fills at least some of the pore space. This layer is usually made up of loose unconsolidated sedimentary materials although igneous, metamorphic, and consolidated sedimentary outcrops can have a weathered layer due to exfoliation and dilatation. Moreover, this layer, which is usually 4 to 60 m thick, is characterized by low seismic velocities which are not only low (250 and 1000 m/sec) but at times highly variable [2,6]. Velocities within the weathered layer vary both vertically and laterally.
The importance of the low velocity layer is four fold [14]: (1) the absorption of seismic energy is high in this zone causing false structures. It is therefore essential that the effects of the weathered layer be eliminated from the records, (2) the low velocity and the rapid changes in velocity have a disproportionately large effect on travel-times. Shooting in the weathering layer results in low energy penetration of the sub-surface for a given source energy level and also excites a greater proportion of boundary waves. These effects degrade the signal-to noise ratio of the resulting recording. Because of this, it is preferably economical to shoot beneath the weathering layer when possible, (3) the marked velocity change at the base of the LVL sharply bends seismic rays so that their travel through the LVL is nearly vertical regardless of their direction of travel beneath the LVL, and (4) the very high impedance contrast at the base of the LVL makes it even an excellent reflector, important in the generation of multiples (at the base where t 1 and t 2 are the intercept times on the time-distance graph, and v o is the velocity of the first weathered layer. The velocities of the subweathered and consolidated layers are v 1 and v 2 respectively, and the sum of D o and D 1 is the total thickness (D w ) of the weathered layer.
For the derivation of these equations, the following assumptions were made. Seismic wave travels along straight raypaths with velocities that are constant within each layer that is isotropic and homogenous. The layers are horizontal bedding, velocity increases with depth, and only P-wave is considered. The recorded travel-time data may be interpreted by the use of graphical method.

Downhole-laid-detectors refraction methods:
The features of the down hole include source/shot point is close to the surface; array of hydrophone-detectors are down the hole [3]. The hole is normally filled with water. In the early 1980s in Nigeria, the practice was to fire shots at various depths in the borehole beginning at the bottom and continuing upwards until the shot is just below the surface of the ground. The field challenge faced ranged from borehole wall collapses after some shots are taken at depth to shot-misfires. This results in the borehole cavingin and blocking the hole. Consequently, there are downtimes and extra costs in re-drilling new holes. Now an array of hydrophone receivers is fixed in the hole as shown in Figure 1. Uko et al. [16] observed that there exists a 7-18% difference between the velocities determined by the surface-wave and downhole methods. The surface-wave method gives spatially averaged velocities along the line of traverse coincident with the geophone spread. Downhole technique is mainly used in marshy swampy terrains where geophones planting become a problem that could lead noisy or dead traces for poor planting. In the case of downhole seismic method, velocities are specifically applicable only to the borehole location. The implication is that downhole measurements can be adversely affected by the local borehole factors and conditions resulting from plastic casings and excessive collapse around the flushedzone of the borehole [16].
Measurements of arrival times of the seismic pulses to the geophones at various depths give the velocity of the propagation of pulse in the ground. Since the elastic properties of the layer are assumed homogenous generally, the velocity of the ground motion will therefore vary both with depth and horizontal distance.
For a source depth of h, an offset of X, the vertical velocity t v corrected for the source depth is [3]: where Z is the respective receiver depths, h is the source depth, X is the offset, t s is travel time for slant path (the raw time picked from the seismogram). The velocities were computed from the reciprocals of the slopes of the straight-line segments based on the thickness equation, t i is the intercept time.

Location, geology and geomorphology of study areas
The Niger delta is located at the southern end of Nigeria. It is a constructive lobate-arcuate delta situated between latitudes 4°-6°N and longitudes 4°-9°E ( Figure 2). It covers a geographical area of 259,000 km [17]. Its subaerial encompasses about 75,000 km 2 and extends 300 km from apex to mouth. It is predominantly a savannah covered lowland covering over a million square kilometres [18]. The basin is bounded on the south by the Atlantic Ocean, on the west by the Benin flank, on the north by the Anambra Embankment and the Abakaliki Anticlinorium and on the east by the Calabar Flank ( Figure 2). The geology of the Niger Delta is only known through the numerous subsurface data acquired during oil prospecting. The history and structure of the Niger Delta are relatively known through several syntheses [19][20][21][22].  It often contains sandstone lenses near the contact with the overlying (Agbada formation) and attains a minimum thickness of 1200 m [18,20]. The Agbada formation consists of upper predominantly sandy unit, minor shale intercalations and a lower unit which is thicker than the upper sandy unit. The sandy beds of the formation form prolific aquifers. The high permeability of the Benin formation, the overlying lateritic red earth and weathered top of the formation as well as the underlying clay-shale member of the Bende-Ameki and Ogwashi-Asaba series provide the hydrologic conditions favoring aquifer formation in the area.
Gravelly beds up to 9 m thick were reported by Amechi [23]. The top sediments are aerated unconsolidated sandstones and have a highly variable thickness throughout the region. This formation consists of sands, clay, lignite, peat and some granular gravel. The clay strata in this formation are of varying thickness from 3 m to 10 m. The sands and sandstone are coarse grained, gravelly and locally fine grained [24]. They range in color from white to yellowish brown. The high yielding permeable Benin formation, the overlying lateritic red earth and weathered top of the formation as well as the underlying clay-shale member of the Bende-Ameki series provide the hydrologic conditions favoring aquifer formation in the area. Generally, the sediments in the Niger Delta range in age from Paleocene to recent [18,20].
The climate of the study area is tropical and dominated by two main seasons: the dry and rainy seasons. Average mean annual rainfall is over 2400 mm and according to hydro-geological studies over 40% of this infiltrate and recharge groundwater [25]. The soil profile is remarkably uniform throughout the area. Approximately, the whole areas consist of deeply weathered and intensely leached soil. The heavy rainfall coupled with the drainage nature of the sub soil is conducive for the high infiltration of rain water.
The proximity of Nigeria to the equator is responsible for the general high temperature. A mean annual temperature of 27°C is recorded in part of the Niger Delta Region. Minimum temperatures in the coastal states are highest in February, March and April and lowest in January and August. Relative humidity near the coast is about 80% to 100% at dawn and 70 to 80% in early afternoon at maximum temperature.
The vegetation cover is thick mangrove in a very swampy environment. Generally speaking, the topography of the Benin formation is flat lying with the greatest height not exceeding 180 m above sea level. The ground level however, slopes gently towards the Niger River flood plains, being directed mainly by the Cross, Imo, Kwaibo and Sambrario River systems. The vegetation is typically tropical rain forest.

Materials and Methods
In the downhole detector array method, hydrophone spread consisting of 24 hydrophones at specified interval was lowered into a hole drilled to 66 m deep. The drilling was by manual rotary method using steel pipes. Water and drilling mud were employed during drilling PVC casings were installed after the drilling was completed. This was to prevent the collapse of the hole. The hydrophone spread was then lowered into the hole. A heavy metal weight was attached to the end of the hydrophone spread to enhance the stability of the spread in the hole and to confirm that the spread reaches the bottom of the hole.
The hydrophones were spaced from the surface and were linked  The picked travel times were corrected to account for the 5 m offset distance from the seismic source to the borehole head. This correction approximates as though the data were recorded with the seismic source placed exactly at the borehole head. The velocities and the LVL thickness were determined from the corrected travel times.
The field set-up of the survey for surface detectors, a spread of 12 geophones was planted on the earth's surface. Shots of 0.2 kg dynamite charge were exploded from offsets of 5 m. The average results for both direct and reverse shots represent the velocity structure of the location. The monitor records were first edited for dead and noisy traces. The arrival times were picked from the traces of the monitor records and used to compute the velocities, using Equation 3.

Results and Discussion
The sample monitor records for the two areas are shown in Figures  4 and 5. The computed elevations, velocities and LVL thickness are presented in Tables 1 and 2. Figures 6 and 7 show typical time-offset (T-X) profiles in Upper Flood Plain and Mangrove Swamp respectively. The velocity of the weathered zone was calculated from the reciprocals of the slopes of the various line segments of the plots using Equations 1 and 2 above. The computed thicknesses, velocities were used for interpretation. The interpreted data showed a substantial variation of the weathered layer thickness and elevation in the two study areas.      [2][3][4][5][6][7][8] who worked on LVL characteristic in different environments in the Niger Delta. The results of this work can used for static correction in seismic processing, for planning and assessing risk for engineering structures, and for groundwater exploration.

Conclusion
Based on analyzed results, the following conclusions are deduced for areas.

5.
The results of this work are very important for static correction in seismic processing, for planning and assessing risk for engineering structures, and for groundwater exploration.