Assessment of Groundwater Potential using Geophysical Techniques in Some Part of OFU Local Government Kogi State, North Central Nigeria
Assessment of Groundwater Potential using Geophysical Techniques in Some Part of OFU Local Government Kogi State, North Central Nigeria
Assessment of Groundwater Potential using Geophysical Techniques in Some Part of OFU Local Government Kogi State, North Central Nigeria
ISSN No:-2456-2165
ABSTRACT
The groundwater potential of Some part in Ofu Local government Kogi State was investigated using vertical electrical
sounding method. A total number of 20 vertical electrical sounding data were acquired using Schlumberger array
configuration. The data obtained were iterated using IPI2WIN software from which the geoelectric parameters (Resistivity,
thickness and depth) were obtained. A total of 4-5 earth modeled layers were delineated which include, Topsoil, clay, sandy
clay, sandstone and water-saturated sand layers were identified. Resistivity values for topsoil ranges from 1103-294 Ωm
with depth varying from of 0.146-1.61m; sandstone with resistivity values 2964-9655 Ωm with depth varying from 0.364-
59.86m; clay with 135-2006 Ωm and depth 0.788-11.0m; sandyclay with 927-7111 Ωm with depth 0.364-59.86m; clayeysand
of 224-2238 Ωm with depth 1.29-10.4m. Overburden thickness varied significantly with a range between 7.33-80m. Thin
overburden layers (<20m) were found at VES 13 and 15.Moderate overburden layers (20-40m). Thick overburden layers
(40-80m) and very thick layers (>80m) were identified at VES 4 and 19. The usage of VES technique has suggested future
boreholes in the area should be drilled at an expected depth of 42.8 m to 112 m depth.
TABLE OF CONTENTS
LIST OF FIGURES
Figure 1.1: Location Map of the Study Area Derives (Adopted from Q gis Version 3)- 5
Figure2.1: Geological map of the Study Area - - - - 12
Figure2.2: Geometric Arrangement of the Schlumberger Array Configuration - 13
Figure 3.4 Resistivity Map of the Study Area- - - - - 24
Figure4.1.1: VES 1 and 2 Data- - - - - - - - 24
Figure4.1.2: VES 3and 4 Data- - - - - - - 25
Figure4.2.1: VES 1Curve - - - - - - - - 26
Figure4.2.2: VES 2 Curve - - - - - - - - 27
Figure4.2.3: VES 3 Curve - - - - -- - - - 28
Figure4.2.4: VES 4 Curve - - - - - - - - - 29
Figure4.3: Estimated Table of Geologic Formation of the Study Area - - -33
Figure 4.4: Resistivity Curve Type of the Study Area- - - - - 34
Figure 4.5: Lithologic Table of the Study Area- - - - - 35
Figure 4.6: Resistivity Depth(s) of the study Area - - - - 37
Figure 4.7: Depth to Groundwater Potential- - - - - -37
Figure 4.8: Depth to Aquifer - - - - - - 39
Figure 4.9: Thickness of Overburden Map - - - - - 40
APPENDIX
CHAPTER ONE
INTRODUCTION
A. Background of the Study
The electrical resistivity method is a geophysical method that has been successfully used by several researchers for
groundwater investigations. The electrical resistivity method (non-invasive) is usually preferable because of the resistivity contrasts
obtained when the groundwater zone is reached. Research in groundwater geophysics reveals that this technique has been
successfully employed for delineating/mapping of geo-materials in the subsurface that guides exploration of groundwater resources
(Metwaly et al., 2012; Oyedele et al., 2009). The electrical resistivity method is the most versatile and the most popular method of
all the geophysical methods that are used in groundwater investigation in Basement Complex area (Olorunfemi, 2009). It is relevant
in the depth to bedrock estimation, aquifer delineation and structural mapping (Ademilua and Eluwole, 2013; Obasi et al., 2013).
Geophysical techniques have been applied to groundwater with some showing more success than others. The methods include
gravity, magnetic, seismic, electrical and electromagnetic methods (Reynolds, 1997 and Somiah, 2013). One of the methods that
have proved effective to groundwater studies are the electrical and electromagnetic. This is because many of the geological formation
properties that are critical to hydrogeology such as porosity and permeability of rocks can be correlated with electrical
Geophysical resistivity methods are based on the effective response of the earth to the flow of subsurface electrical current.
Electrical resistivity is widely used for hydrogeological studies, because the acquired data are mainly controlled by lithological
conditions of the aquifer. The method is also useful in the correlation of lithological facies between wells (Obianwu et al. 2015).
VESs, when combined with other geophysical methods, geologic mapping and available well data can greatly assist in the location
and completion of water wells in bedrock areas of complex hydrogeology. The VES method is usually considered more suitable for
the subsurface investigation of geologic environments consisting of horizontal or nearly horizontal layers, such as occur in
unconsolidated sedimentary sequences (Ojekunle et al. 2015; Badmus and Olatinsu 2012; Alile et al. 2008).
Groundwater is used for various purposes in both the developed and developing countries. Among its major purposes are for
use in municipal, agriculture and industries, often extracted through the construction of wells and boreholes. In most parts of the
countries where there is no (or limited) supply of freshwater from the surface water sources, groundwater serves the alternative
natural sources. In Nigeria, large portions of its freshwater usage come from groundwater sources. The regional distribution of
groundwater within the earth varies and this depends on the geology of the subsurface.
The reliance on groundwater is such that it is necessary to ensure that there are significant quantities of water and that the water
is of a high quality. The principal source of groundwater is meteoric water, that is, precipitation. However, two other sources are
occasionally of some consequence. These are juvenile water and connate water. The former is derived from magmatic sources,
whereas the latter represents the water in which sediments are deposited. Connate water is trapped in the pore spaces of sedimentary
rocks as they are formed and has not been expelled.
The amount of water that infiltrates into the ground depends on how precipitation is dispersed, namely, on the proportions that
are assigned to immediate run-off and to evapotranspiration, the remainder constituting the proportion allotted to
infiltration/percolation. Infiltration refers to the seepage of surface water into the ground, percolation being its subsequent
movement, under the influence of gravity, to the zone of saturation. In reality, one cannot be separated from the other. The infiltration
capacity is influenced by the rate at which rainfall occurs (which also affects the quantity of water available), the vegetation cover,
the porosity of the soils and rocks, their initial moisture content and the position of the zone of saturation. The retention of water in
soil depends on the capillary force and the molecular attraction of the particles. As the pores in the soil become thoroughly wetted,
the capillary force declines, so that gravity becomes more effective. In this way, downward percolation can continue after infiltration
has ceased but the capillarity increases in importance as the soil dries. No further percolation occurs after the capillary and gravity
forces are balanced. Thus, water percolates into the zone of saturation when the retention capacity is satisfied. This means that the
rains that occur after the deficiency of soil moisture has been catered for are those that count as far as supplementing groundwater
is concerned. The water in the zone of aeration is referred to as vadose water. This zone is divided into three subzones, those of soil
water, the intermediate belt and the capillary fringe. The uppermost or soil water belt discharges water into the atmosphere in
perceptible quantities by evapotranspiration. In the capillary fringe, this occurs immediately above the water table, water is held in
the pores by capillary action.An immediate belt occurs when the water table is far enough below the surface for the soil water belt
not to extend down to the capillary fringe. The degree of saturation decreases from the water table upwards. An aquifer is the term
given to a rock or soil mass that not only contains water but from which water can be abstracted readily in significant quantities.
The ability of an aquifer to transmit water is governed by its permeability.By contrast, a formation with permeability of less than
10-9 m is one that, in engineering terms, is regarded as impermeable and is referred to as an aquiclude. For example, clays and
shales are aquicludes. Even when such rocks are saturated, they tend to impede the flow of water through stratal sequences. An
aquitard is a formation that transmits water at a very slow rate but that, over a large area of contact, may permit the passage of large
amounts of water between adjacent aquifers that it separates.
Fig 1: Location Map of the Study Area (Adopted from Q-GIS V 3.0)
To provides valuable information about the availability, distribution, and quality of groundwater resources in the region.
To access to clean and reliable water sources is a fundamental requirement for the socioeconomic development of any region.
Identifying potential aquifers and zones with significant groundwater resources can support local communities in meeting their
water needs, which in turn promotes better health, education, and overall quality of life (Ayoade and Akintola, 2012).
CHAPTER TWO
LITERATURE REVIEW AND GEOLOGICAL SETTING
Alile et al. (2011) conducted a geoelectric investigation to identify potential groundwater resources in Obaretin - Iyanomon
locality, Edo State, Nigeria. The researchers employed the Vertical Electrical Sounding (VES) method, acquiring data from three
locations within the study area. The study revealed three to four subsurface geologic layers characterized by varying resistivity and
thicknesses and interpreted these layers as sandy topsoil, sandy clay, sand, and coarse sand or fractured bedrock, respectively. Based
on the resistivity values and layer thicknesses, the authors identified potential aquifer zones within the study area.
Adeniran et al. (n.d.) investigated the hydrogeophysical implications of geoelectric sounding at Igarra Comprehensive High
School, located in Akoko Edo Local Government Area of Edo State, Nigeria. The study aimed to determine the subsurface geology
and its impact on groundwater occurrence and movement in the area. The study used the VES method to collect geoelectric data
from various locations within the study area. The resulting data revealed the presence of three to four subsurface geoelectric layers,
characterized by distinct resistivity values and thicknesses. These layers were interpreted as topsoil, weathered/fractured basement,
and fresh basement rocks. Based on their findings, Adeniran etal. (n.d.) identified potential aquifer zones within the weathered and
fractured basement layers, emphasizing the importance of these geological features in trapping and storing groundwater.
Ojo et al. (2017) conducted a study to determine the groundwater potential in the permanent site of the University of Abuja,
located in Nigeria's Federal Capital Territory (FCT). The authors used the Vertical Electrical Sounding (VES) technique to collect
geoelectric data from ten locations within the study area. The collected data were interpreted using the Psartial Curve Matching
(PCM) method and numerical modeling software.The results revealed three to five subsurface geoelectric layers characterized by
different resistivity values and thicknesses. These layers were identified as topsoil, weathered basement, fractured basement, and
fresh basement. The authors identified potential aquifer zones in the weathered and fractured basement layers, Ojo etal. (2017)
concluded that the groundwater potential in the study area is relatively high and that the identified aquifer zones can be targeted for
sustainable groundwater exploitation.
Saleh et al. (2019) investigated the groundwater potential of Kaltungo and its environs in Northeastern Nigeria using
geoelectric techniques. The authors employed the VES method to collect geoelectric data from twenty locations across the study
area. The data were analyzed and interpreted using the PCM method and 1D-Inversion software.The interpreted results revealed
three to five subsurface geoelectric layers, including the topsoil, lateritic clay, weathered/fractured basement, and fresh basement.
Based on the resistivity values and layer thicknesses, Saleh et al. (2019) identified potential aquifer zones in the weathered and
fractured basement layers. The study concluded that the groundwater potential in Kaltungo and its environs is moderate to high and
that the identified aquifer zones can be targeted for groundwater exploitation.
Nwachukwu et al. (2013) evaluated the groundwater potentials of Orogun in the South-South region of Nigeria using the
electrical resistivity method. The study aimed to characterize the subsurface geology and identify potential aquifer zones in the
area.Nwachukwu et al. (2013) identified potential aquifers in the sand layer.
Badmus and Olatinsu (2010) carried out a research work on aquifer characteristics and groundwater recharge pattern in a
typical basement complex: A case study of Federal College of Education, Osiele, Abeokuta, South west of Nigeria. In their work,
they used vertical electrical sounding (VES) using the Schlumberger electrode array. It was revealed that Abeokuta has seven major
geological formations which are topsoil, shale or clay, sandy clay, clayey sand, sandstone, fractured basement and fresh basement.
It also was discovered that the weathered and fractured basement constitute the main aquiferous units in the area. They were able to
discover that the reasons for borehole failure and poor recharge in the area were due to inadequate geophysical investigation, the
depth at which drilling was terminated and the geological formation of the aquifers. They also engaged 3-D view to show the
overburden thickness.
Of all surface geophysical methods, the electrical resistivity method has been applied most widely for ground water
investigations out of the entire surface geophysical methods. The electrical resistivity method can be best employed to estimate the
thickness of overburden and the thickness of weathered/fractured zones with reasonable accuracy. Though both Wenner and
Schlumberger electrode configuration methods are popularly employed, the Schlumberger electrode configuration method is more
suited to the study area, ensuring better results. The method has practical, operational and interpretational advantages over other
methods such as the Wenner method of electrode arrangement (Zohdy et al., 1974). Most of the electrical resistivity techniques
require injection of electrical currents into the subsurface via a pair of electrodes planted on the ground. By measuring the resulting
variations in electrical potential at other pairs of planted electrodes, it is possible to determine the variations in resistivity. A
conventional vertical electrical sounding (VES) survey was used for quantitative interpretation where the center point of the array
remains fixed and the electrode spacing is increased for deeper penetration (Loke,1999).
In carrying out resistivity sounding survey, electrodes are distributed along a line centered about a midpoint that is considered
the location of the sounding. The electrode arrangement used in data acquisition is the Schlumberger array of electrodes. The
Schlumberger survey involves the use of two current electrodes labeled A and B, and two potential electrodes M and N placed in
line with one another and centered on some location. It is worthy to state that the potential and current electrodes are not placed
equidistant from one another. To acquire the resistivity data in the field, current is introduced into the ground through the current
electrodes and the potential electrodes are then used to quantitatively measure the voltage pattern. The apparent resistivity data
obtained from the measurements are presented on maps at various levels and they are useful in the first stage of interpretation. More
realistic sections of the earth are obtained only after interpretation of the data in terms of true variations of the resistivity distribution.
This is a very important step because it allows the estimation of the true position and depth of formations. Moreover, it is possible
to estimate the actual electrical resistivity of the region and relate it to its physical state.
D. Resistivity Method
The electrical resistivity method measures potential differences at points on the Earth’s surface that are produced by directing
current flow through the subsurface. This leads to the determination of resistivity distribution in the subsurface and to an
interpretation of. Earth materials. The movement of charges through the conducting wire is termed current.
Specifically,
Where I is current in amperes, Q is charge in coulombs, and t is time in seconds. Also, another important concept in electrical
resistivity surveying is the current density J, which is defined as the current divided by the cross-sectional area of the material
through which it is flowing
J=I/A……………………………………………………………………………………………………………………………….(2.1)
I=V/R………………………………………………………………………………………………………………………………(2.2)
V=IR…………………………………………………………………………………………………………………………………..(2.3)
Electrical resistivity method is basically conducted to measure and map the resistivity of subsurface materials. It also refers to
survey that is carried out to present the image of electrical properties of the subsurface by passing an electrical current along many
different paths and measuring the associated voltage. Resistivity method is based on the response between the earth and the flow of
electrical current. It sensitive to variations in the electrical resistivity of the subsurface measured in Ohms meter. Resistivity
measurement are conducted by inducing an electric current into the earth through two current (C1 and C2) electrodes and measuring
the resulting voltage at two potential electrodes (P1 and P2). The apparent resistivity (pa) value can be calculated based on the
current (I) and voltage (V).
………………………………………………………………………………………………………………….……(3.1)
The resistance is also proportional to the cross sectional area and the distance between the electrodes. The relationship is given
by:
……………………………………………………………………………………………………………………… (3.2)
……………………………………………………………………………………………………………………........ (3.3)
V is the voltage,
I is the current and
L is the distance between the electrodes, the constant of proportionality, is the apparent resistivity.
Data from resistivity surveys are represented by apparent resistivity which takes into account, the arrangement and spacing of
electrodes.
From the relationship above, the potential at any point is given by:
V= …………………………………………………………………………………………………………………………..(3.4)
Where: V is the potential in volts is the resistivity of the medium and r is the distance from the electrode.
The resistivity of the ground is calculated from the potential difference between P1 and P2. The potential VP1 at the internal
electrode P1 is given by the algebraic sum of the potential contributions VC1 and VC2 from the current source at C1 and the sink at
C2.
VP2 ……………………………………………………………………………………………….(3.6.3)
Because various ecologic materials can be expected to have different resistance to current flow, it might seem fairly straight
forward to measure current and voltage to calculate resistance and determine the material in the subsurface. One immediate
complication is that resistance depends not only on the material but also on its dimensions.
𝐼𝜌 𝐼𝜌 𝐼𝜌 𝐼𝜌
( − )−( − ) … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … . . . (3.6.4)
2𝜋𝑟1 2𝜋𝑟2 2𝜋𝑟3 2𝜋𝑟4
The two pairs of electrodes P1 and P2 (Fig. 2.3) carry no current but are used to measure the potential difference between the
points P1 and P2.
………………………………………………………..……(4.3)
Where
k=2
is the geometric factor, and is only a function of the geometry of the electrode arrangement
Resistivity can be found from measuring values of V, I and K. So the apparent resistivity ( a) equation becomes; In the
resistivity method current is entered into the ground, potential difference ismeasured and resistivity is determined. For 𝜌, equation
7 becomes.
2𝜋∆𝑣 1 1 1 1 ……………………………………………………………………………………………………………………………………………………………
𝜌= ( − − + )-1 (5)
𝐼 𝑟
1 𝑟 𝑟
2 𝑟 3 4
Where R is the measured resistance and G is the Geometric Constant which is a function of the electrode configuration
employed during the survey.
CHAPTER THREE
MATERIAL AND METHOD
A. Introduction
The electrical resistivity method measures potential differences at points on the Earth’s surface that are produced by directing
current flow through the subsurface. This leads to the determination of resistivity distribution in the subsurface and to an
interpretation of earth materials. The movement of charges through the conducting wire is termed current.
B. Material
Ohmega 1000C Terrameter: - This equipment involves VES and electrical profiling method that is based on the four-electrode
principle used by geophysicist for evaluation of greater penetration depth of the injected current for successive readings. The
current injected into the ground could be DC or AC.
Electrode: - This enables the electrode convey current into the subsurface when driven into the earth with the hammer. Two
pairs of electrodes were used; a pair is current electrodes while the other pair is potential electrodes.
Reels and Cables: It is usually connected with a clip attached to one end of the electrode while the other end connects sat the
terrameter. The cable is made of an insulating material outside housing the flexible copper wire. Due to the flexibility of the
cable, it can be easily recoiled back into its reel.
Global Positioning System (GPS): - Garmin GPS72H GPS was used to take the coordinates as well as the elevation of the
surveyed areas. This is necessary so that the exact position of each VES point could be located in the future if the need arises.
Hammer: This was used to drive the electrodes into the ground.
C. Research Methodology
There are two basic procedures in resistivity work. The procedure to be used depends on whether we are interested in lateral
or vertical variations in resistivity. Resistivity surveys are made to satisfy the needs of two distinctly different kinds of interpretation
problems:
The variation of resistivity with depth, reflecting more or less horizontal stratification of earth materials; and
Lateral variations in resistivity that may indicate soil lenses, isolated ore bodies, faults, or cavities. For the first kind of problem,
measurements of apparent resistivity are made at a single location (or around a single center point) with systematically varying
electrode spacing.
D. Schlumberger Array
This type of configuration deals with vertical electrical sounding. Electrical sounding is the process by which the variation of
resistivity with depth below a given point on the ground surface is deduced and it can be correlated with the available geological
information in order to infer the depths (or thicknesses) and resistivity of the layers (formations) present.
The procedure is based on the fact that the current penetrates continuously deeper with the increasing separation of the current
electrodes.
The current electrodes are spaced much further apart than the potential electrodes. The potential electrodes remain fixed while
the current electrodes are expanded symmetrically about the center of the spread.
E. Field Procedures
The Schlumberger array is a commonly used electrode configuration in electrical resistivity surveys for both vertical electrical
sounding (VES) and profiling. The array consists of four collinear electrodes: two outer current electrodes (A and B) and two inner
potential electrodes (M and N). The current is transmitted through the outer electrodes, and the resulting voltage is measured between
the inner electrodes.
Select the survey area: Choose the location where you want to investigate the subsurface resistivity. The area should be clear
of any obstructions, such as trees or buildings, and the ground should be relatively flat
Set up the equipment: Connect the resistivity meter or multi-meter to the electrodes using appropriate cables. Ensure that the
meter is set to measure resistance or resistivity and that it is calibrated correctly.
Establish the initial electrode spacing: The outer electrodes (A and B) were positioned at a specified distance from each other
(AB), typically several times the expected investigation depth.The inner electrodes (M and N) were also positioned at a smaller
distance (MN) from each other, usually between 10% and 30% of AB.
Take measurements: Switch on the resistivity meter and allow it to stabilize. Record the measured resistance or resistivity value
and the corresponding electrode spacing (AB and MN). Ensure that the measured value is stable before recording it.
Increase electrode spacing: Move the outer electrodes (A and B) further apart, and keep the inner electrodes (M and N) at a
constant distance from each other. Take another measurement at this new electrode spacing and record the value.
Repeat the process: Continue increasing the distance between the outer electrodes and taking measurements until you reach the
desired depth of investigation. Keep the inner electrode spacing constant throughout the survey.
Analyze the data: Plot the measured resistivity values as a function of electrode spacing. The resulting curve can be interpreted
to determine the subsurface resistivity distribution, which provides information on the geological layers and the presence of
groundwater or other subsurface features.
Profiling (optional): If profiling is required, keep the outer electrodes (A and B) at a fixed distance from each other and move
the inner electrodes (M and N) along the survey line, taking measurements at regular intervals. This provides information on
lateral resistivity variations in the subsurface.
Pack up the equipment: Once the survey is complete, disconnect the cables, remove the electrodes from the ground, and pack
up the resistivity meter and other equipment.
The Schlumberger array is advantageous for VES and profiling due to its efficient use of personnel and time during field surveys.
Additionally, the variable electrode spacing allows for detailed investigations of specific depths, making it suitable for various
geological and hydrogeological application.
CHAPTER FOUR
RESULTS AND DISCUSSION
The VES data obtained were in resistance (Ω), using ohmega 1000c terrameter which were then multiplied by their
corresponding geometric factor (k) in order to obtain an apparent resistivity in ohm-meter. Computer software called IPI2WIN was
used to processes the data, the apparent resistivity values for each point was plotted against the half current electrode spacing AB/2
on log-graph to obtain sounding curves. The purpose of these curves is to determine the subsurface layers beneath each VES point
as well as their thickness and resistivity variation, the results of the analysis of data were correlated with a borehole log data obtained
near the study area. Boreholes log data are a necessary and reliable source of data, and electrical resistivity method using vertical
electrical sounding (VES) interpretations provide basic information of the area. The data was firstly processed and analyzed using
IPI2win version 3.0.1 Software which gave automatic geo-electric parameters of the study area. The output of VES data is resistivity
layer, log resistivity graph and resistivity depth table. Sounding curve obtained from the study area showed variation of four - five
geo-electrical layers.
ELECTRODE
ELECTRODE
SPACING
SPACING
RESISTIVITY
RESISTIVITY
RESISTANCE
RESISTANCE
GEOMETRIC
GEOMETRIC
FACTOR “K”
FACTOR “K”
APPARENT
APPARENT
(OHMS M)
(OHMS M)
(M)
(M)
(OHMS)
(OHMS)
S/NO.
S/NO.
AB/2 MN/2 AB/2 MN/2
Fig 4.2.2 displays KHK curve type with seven different layers which indicates different lithology by the resistivity values
(Ohm) in location Giwoligyo(m)and depth(m) of each of these layers are displayed as well as one dimensional pseudo section which
image the distribution of the VES in the subsurface.
Fig 4.2.3 displays an HKH curve type with five different layers which indicates different lithology by the resistivity
values(Ohm) in location Ofokofi, thickness(m)and depth(m) of each of these layers are displayed as well as one dimensional pseudo
section which image the distribution of the VES in the subsurface.
Fig 4.2.4 displays KHK curve type with seven different layers which indicates different lithology by the resistivity values
(Ohm) in location Gwalawo, thickness (m) and depth (m) of each of these layers are displayed as well as one dimensional pseudo
section which image the distribution of the VES in the subsurface. Thicker aquifers with moderate to low resistivity values are more
likely to yield significant water flows. Based on the analysis, Fig 4.2.2 emerges as one of the promising drilling target.
Start drilling at Fig 4.2.4 (GWALAWO) to a depth of 138 m, targeting the potential aquifer in Layer 6, Consider drilling at
VES 14 (EFAKU) to a depth of 120 m, targeting the water-saturated sand layer in Layer 4.Drill at VES 16 (OLAJI) to a depth of
45 m, targeting the potential aquifer in Layer 4.If the above locations do not yield satisfactory results, consider drilling at VES 19
(CHIKORO) to a depth of 212.3 m.
The VES curve obtained from plotting the apparent resistivity against the corresponding half of electrode spacing (AB/2) gives
curves types such a s HAK, HKQ, AAH, AAK, AKH,KHK and HAA. Most of the obtained sounding curves were of the HKH type
(D1>D2<D3>D4<D5) as shown and the HKQ type (D1>D2<D3>D4) . These curves types indicated five to eight lithologies. The
HAK curves rose steeply into positive slopes and such curves are a reflection of a highly resistive sedimentary rock at depth
which 20 serves as underlying beds for unconfined aquifers. Fig 4.3.4 shows the curve types, sequence and number of layers for the
four resistivity sounding curves, while others can be seen in the appendices.
Lithology Description
In fig 4.4 each VES location exhibits a unique layered sequence, comprising laterite, clay, sandyclay, sandstone, and water-
saturated sand layers. The resistivity values range from low (<100 Ωm) to extremely high (>100000 Ωm), indicating varying degrees
of compactness, fracturing, and water saturation. The reconstruction was done by interpreting the resistivity values from each layer
curves, the resistivity and thickness of the individual layers. It was noted that different layer exhibit same resistivity. The
interpretation was solved b y understanding the deeper knowledge of the underlying geology. The results of the quantitative
interpretation revealed that typical 5- and 8 geoelectrical subsurface layers characterize the area, generating eight types of resistivity
curves. The surface underlying laterite is a layer correlated to reddish brown, sandy clay and clayey sands, with thickness range
from 1.4 m to 8.4 m, averaged 4.68 m, and resistivity range with an average of 168.94 Ωm
Spatial Distribution
The spatial distribution of groundwater depths reveals an increase in depth towards the eastern part of the study area. This
suggests a potential groundwater flow direction, with recharge areas located in the western and central regions. Shallow groundwater
zones (<50m) may require shorter drilling depths, reducing costs. Moderate groundwater zones (50-100m) require standard drilling
depths. Deep groundwater zones (100-150m) require deeper drilling, increasing costs. Very deep groundwater zones (>150m) may
require specialized drilling equipment.
Spatial Distribution
The spatial distribution of overburden thickness reveals a general increase in thickness towards the eastern part of the study
area. This suggests a potential relationship between overburden thickness and groundwater flow direction.
Thin overburden zones may require shorter drilling depths, reducing costs. Moderate overburden zones require standard
drilling depths. Thick overburden zones require deeper drilling, increasing costs.
Thin Overburden Layers (<25m): VES 13 (19.1m), VES 14 (18.8m), VES 16 (15.5m), VES 20 (3.5m), and VES 1 (8.024m)
exhibit relatively thin overburden layers. These areas likely have shallow water tables or high permeability, facilitating rapid
groundwater recharge.
Moderate Overburden Layers (25-50m): VES 2 (22.3m), VES 3 (34.13m), VES 6 (21.13m), VES 7 (29m), VES 9 (33.7m),
VES 10 (11.3m), VES 12 (37.1m), and VES 18 (37m) display moderate overburden thicknesses. These locations are scattered
throughout the study area, indicating a relatively uniform distribution.
Thick Overburden Layers (50-75m): VES 5 (59.86m), VES 8 (33.7m), and VES 11 (54.2m) exhibit thicker overburden
layers. These areas likely have low permeability or high clay content, restricting groundwater flow.
Very Thick Overburden Layers (>75m): VES 4 (95.3m) and VES 19 (113m) shows exceptionally thick overburden layers.
These locations likely represent areas with complex geological structures or high water tables, influencing groundwater flow and
recharge.
Spatial Distribution
The spatial distribution of overburden thickness reveals a general increase in thickness towards the eastern part of the study
area. This suggests a potential similarities between overburden thickness and groundwater flow direction.
CHAPTER FIVE
SUMMARY, CONCLUSION AND RECOMMENDATION
A. Summary
Geophysical investigation using Vertical Electrical Sounding (VES) was adopted in Some part of Ofu Local government Kogi
State for the exploration of groundwater potential in the area. The interpreted results of the VES data were done using IPI2WIN and
suffer software packages and the results are presented as sounding curves, geo-electric sections and contour map . A total of 4-5
earth modeled layers were delineated which include, Topsoil, clay, sandy clay, sandstone and water-saturated sand layers were
identified. Resistivity values for topsoil ranges from 1103-294 Ωm with depth varying from of 0.146-1.61m; sandstone with
resistivity values 2964-9655 Ωm with depth varying from 0.364-59.86m; clay with 135-2006 Ωm and depth 0.788-11.0m; sandyclay
with 927-7111 Ωm with depth 0.364-59.86m; clayeysand of 224-2238 Ωm with depth 1.29-10.4m.
The resistivity parameters of the geo electric layers across the entire area were used to delineate the parameters of the identified
aquifers in the study area.
B. Conclusion
20 Schlumberger vertical electrical soundings (VES) were used in the geophysical research survey. The designation of an
aquifer unit in the research area has been made possible by the survey's findings. In the study region, three (4) to five (5) unique
geo-electric layers were identified, including the lateritic sand/topsoil, sandy clay, clay, Clayey sand ,sandstone, and Water saturated
sand tone indicated in fig 4.5. The productive aquifer unit in the area has been identified by survey results. According to computer
study, the specific layers, which is made up of sandstones, where a groundwater search is necessary in the research region, its
subsurface soil type has been identified, and geophysical data on the groundwater has been provided.
C. Recommendation
It is recommended that a thorough study of any project be done before any operation so as to know whether the aquifer is
naturally sealed or not. And the presence of clay above the aquifer must be sought for before choosing any site for drilling.
REFERENCES
[1]. Ehinola, O. A., Olasehinde, P. I., Agbaire, A. O., Osagiede, J. O., Akpata, J. O., &Ebunoluwa, A. O. (2014). Geology and
mineral resources of Edo State, Nigeria. Journal of Global Biosciences, 3(5), 711-718. Retrieved from
https://globalbiosciences.org/vol3no5/article-8.pdf
[2]. Etu-Efeotor, J. O., Igbinosa, E. I., Adeyemi, A. O., &Ifeduba, C. E. (2020).
[3]. Vegetation assessment of Sobe and Ologbo forest reserves in Edo State, Nigeria. International Journal of Forestry
Research, 2020, 1-11. https://doi.org/10.1155/2020/9541418 Nwajide, C. S., &Reijers, T. J. A. (1996).
[4]. The stratigraphic evolution and hydrocarbon potential of rift basins in Nigeria. Journal of African Earth Sciences, 22(3), 295-
307. https://doi.org/10.1016/0899-5362(95)00047-4 Ogezi, A. E. (1977).
[5]. The geology of part of the Basement Complex of Nigeria with special reference to the petrology and chemical evolution of
the Banded Gneiss Complex of Egbe-Isanlu area (Unpublished doctoral dissertation). University of Leicester, Leicester, UK.
[6]. Olabode, B., Omosanya, K. O., Oyawoye, M. O., &Fadipe, B.Edo State Government. (2018). Edo State Investors' Guide
2018. Edo State Investment Promotion Office. Retrieved from https://www.edostate.gov.ng/wp-
content/uploads/2018/11/Edo-State-Investors-Guide.pdf
[7]. Imoleayo, J. O., Nwankwo, N. E., James, E. O., &Usip, E. E. (2018). Wind speeds variability at different heights in Southern
Nigeria. International Journal of Climatology, 39(3), 1695-1703. https://doi.org/10.1002/joc.5663
[8]. Köppen-Geiger climate classification. (n.d.). Climate-Data.org. Retrieved from https://en.climate-
data.org/africa/nigeria/edo/benin-city/benin-city-298068/
[9]. Oguntunde, P. G., Fadairo, J. O., &Adetunji, A. R. (2015). Characterization of rainfall variability over Southern Nigeria in
the period 1961-2010. American Journal of Climate Change, 4(2), 92-102. https://doi.org/10.4236/ajcc.20American Journal
of Climate Change, 4(2), 92-102. https://doi.org/10.4236/ajcc.2015.42010
[10]. Onyekwelu, J. C. (2011). Floristic composition and vegetation structure of a disturbed tropical rainforest ecosystem in Edo
state, southern Nigeria. Journal of Forestry Research, 22(4), 671-679. https://doi.org/10.1007/s11676-011-0169-9
[11]. Ojo, E. O., Ibitoye, M. O., &Adebisi, M. A. (2017). Geophysical investigation of groundwater potential of the permanent
site of University of Abuja, FCT, Nigeria. Egyptian Journal of Remote Sensing and Space Science, 20(2), 187-196.
https://doi.org/10.1016/j.ejrs.2017.02.003
[12]. Saleh, S., Ejeh, I. U., &Chinedu, N. C. (2019). Geoelectric investigation for groundwater potential of Kaltungo and Environs,
Northeastern Nigeria. Journal of Geoscience and Environment Protection, 7(9), 66-80.
https://doi.org/10.4236/gep.2019.79006
[13]. Nwachukwu, S. N., Bello, R. A., Omosun, O. A., & Daniel, E. O. (2013). Evaluation of groundwater potentials of Orogun,
South-South, Nigeria using electrical resistivity method. Research Journal of Applied Sciences, Engineering and Technology,
5(12), 3075-3081. https://doi.org/10.19026/rjaset.5.3978
ELECTR APPENDIX
ELECTR
SPACIN
SPACIN
RESISTIVITY
RESISTIVITY
RESISTANCE
RESISTANCE
GEOMETRIC
GEOMETRIC
FACTOR “K”
FACTOR “K”
ODE
ODE
APPARENT
APPARENT
(M)
(M)
(OHMS M)
(OHMS M)
G
G
(OHMS)
(OHMS)
S/NO.
S/NO.
AB/2 MN/2 AB/2 MN/2
Y (OHMS M)
Y (OHMS M)
ELECTR
ELECTR
RESISTIVIT
RESISTIVIT
RESISTANC
RESISTANC
GEOMETRI
GEOMETRI
APPARENT
APPARENT
SPACIN
SPACIN
C FACTOR
C FACTOR
ODE
ODE
(OHMS)
(OHMS)
(M)
(M)
S/NO.
S/NO.
G
G
“K”
“K”
E
E
AB/2 MN/2 AB/2 MN/2
1 2 0.5 11.78 69.43 817.25 1 2 0.5 11.78 16.17 190.4826
2 3 0.5 27.5 17.90 492.25 2 3 0.5 27.5 5.606 154.165
3 6 0.5 112.3 6.695 751.85 3 6 0.5 112.3 2.614 293.5522
4 9 0.5 254 4.446 1129.28 4 9 0.5 254 2.055 521.97
5 9 2.0 60.5 8.465 512.13 5 9 2.0 60.5 17.19 1039.995
6 15 2.0 173.6 5.240 909.66 6 15 2.0 173.6 13.32 2312.352
7 25 2.0 488 4.039 1971.03 7 25 2.0 488 9.015 4399.32
8 40 2.0 1253 2.187 3415.68 8 40 2.0 1253 5.667 7100.751
9 50 2.0 1960 2.726 4286.52 9 50 2.0 1960 4.436 8694.56
10 75 2.0 4410 1.566 6906.10 10 75 2.0 4410 2.686 11845.26
11 75 10.0 868 5.708 4954.54 11 75 10.0 868 13.02 11301.36
12 100 10.0 1555 4.395 6834.23 12 100 10.0 1555 10.02 15581.1
13 150 10.0 3520 2.503 8810.56 13 150 10.0 3520 5.809 20447.68
14 200 10.0 6270 1.638 10270.26 14 200 10.0 6270 3.510 22007.7
15 300 10.0 14120 0.839 11846.68 15 300 10.0 14120 1.648 23269.76
16 300 20.0 7040 1.648 11601.92 16 300 20.0 7040 2.024 14248.96
17 400 20.0 12530 1.261 15800.33 17 400 20.0 12530 1.129 14146.37
LOCATION: EJIKOLU LOCATION: ALOKURA
TOWN: EJIKOLU VES 7 TOWN: ALOKURA VES 8
L.G.A: OFU L.G.A: OFU
STATE: KOGI STATE: KOGI
LONGITUDE: E007009'. 219'' LONGITUDE: E0070 04'.321''
LATITUDE: N06059'. 526'' LATITUDE: N060 59'. 583''
ELEVATION: 267 m ELEVATION: 278 m
ELECTRODE
ELECTRODE
SPACING
SPACING
RESISTIVITY
RESISTIVITY
RESISTANCE
RESISTANCE
GEOMETRIC
GEOMETRIC
FACTOR “K”
FACTOR “K”
APPARENT
APPARENT
(OHMS M)
(OHMS M)
(M)
(M)
(OHMS)
(OHMS)
S/NO.
S/NO.
AB/2 MN/2 AB/2 MN/2
1 2 0.5 11.78 221.6 2610.45 1 2 0.5 11.78 23.40 275.65
2 3 0.5 27.5 168.2 4625.5 2 3 0.5 27.5 6.379 175.42
3 6 0.5 112.3 61.27 6880.62 3 6 0.5 112.3 3.367 378.11
4 9 0.5 254 46.90 11912.6 4 9 0.5 254 1.973 501.14
5 9 2.0 60.5 52.40 3170.2 5 9 2.0 60.5 6.837 413.64
6 15 2.0 173.6 18.72 3249.79 6 15 2.0 173.6 4.059 704.64
7 25 2.0 488 8.526 4160.69 7 25 2.0 488 2.482 1211.22
8 40 2.0 1253 4.368 5473.10 8 40 2.0 1253 1.556 1949.67
9 50 2.0 1960 3.811 7469.56 9 50 2.0 1960 1.221 2393.16
10 75 2.0 4410 1.008 4445.28 10 75 2.0 4410 0.793 3497.13
11 75 10.0 868 2.299 1995.53 11 75 10.0 868 4.283 3717.64
12 100 10.0 1555 1.699 2641.95 12 100 10.0 1555 3.093 4809.62
13 150 10.0 3520 1.119 3938.88 13 150 10.0 3520 2.004 7054.10
14 200 10.0 6270 0.8761 5493.15 14 200 10.0 6270 1.221 7655.67
15 300 10.0 14120 0.5361 7569.73 15 300 10.0 14120 0.513 7243.56
16 300 20.0 7040 1.236 8701.44 16 300 20.0 7040 0.969 6821.76
17 400 20.0 12530 0.6343 7947.78 17 400 20.0 12530
LOCATION: ASADAM LOCATION: SANIOKO
TOWN: ASADAM VES 9 TOWN: SANIOKO VES 10
L.G.A: OFU L.G.A: OFU
STATE: KOGI STATE: KOGI
LONGITUDE: E007004'. 277'' LONGITUDE: E0070 05'. 494''
LATITUDE: N06056'. 385'' LATITUDE: N060 56'. 579''
ELEVATION: 202 m ELEVATION: 275 m
ELECTROD
ELECTROD
E SPACING
E SPACING
RESISTIVITY
RESISTIVITY
RESISTANCE
RESISTANCE
GEOMETRIC
GEOMETRIC
FACTOR “K”
FACTOR “K”
APPARENT
APPARENT
(OHMS M)
(OHMS M)
(M)
(M)
(OHMS)
(OHMS)
S/NO.
S/NO.
AB/2 MN/2 AB/2 MN/2
1 2 0.5 11.78 40.49 476.97 1 2 0.5 11.78 22.28 262.46
2 3 0.5 27.5 19.84 545.6 2 3 0.5 27.5 8.750 240.63
3 6 0.5 112.3 7.936 891.21 3 6 0.5 112.3 3.642 408.99
4 9 0.5 254 4.171 1059.43 4 9 0.5 254 2.289 581.41
5 9 2.0 60.5 19.12 1156.76 5 9 2.0 60.5 7.061 427.20
6 15 2.0 173.6 11.09 1925.22 6 15 2.0 173.6 3.947 685.20
7 25 2.0 488 5.809 2834.79 7 25 2.0 488 2.248 1097.02
8 40 2.0 1253 3.846 4819.04 8 40 2.0 1253 1.231 1542.44
9 50 2.0 1960 3.062 6001.52 9 50 2.0 1960 0.914 1791.44
10 75 2.0 4410 1.862 8211.42 10 75 2.0 4410 0.473 2085.93
11 75 10.0 868 8.343 7241.72 11 75 10.0 868 2.167 1880.96
12 100 10.0 1555 6.847 10647.09 12 100 10.0 1555 1.302 2024.61
13 150 10.0 3520 2.503 8810.56 13 150 10.0 3520 0.819 2882.88
14 200 10.0 6270 1.271 7969.17 14 200 10.0 6270 0.473 2965.71
15 300 10.0 14120 0.588 8231.96 15 300 10.0 14120 0.237 3346.44
16 300 20.0 7040 1.398 9841.92 16 300 20.0 7040
17 400 20.0 12530 17 400 20.0 12530
LOCATION: EJULE LOCATION: IBOKO
TOWN: EJULE VES 11 TOWN: IBOKO VES 12
L.G.A: OFU L.G.A: OFU
STATE: KOGI STATE: KOGI
LONGITUDE: E007007'. 495'' LONGITUDE: E0070 07'. 030''
LATITUDE: N07055'. 509'' LATITUDE: N070 55'. 958''
ELEVATION: 290m ELEVATION: 334 m
ELECTROD
ELECTROD
E SPACING
E SPACING
GEOMETRIC FACTOR “K”
APPARENT RESISTIVITY
(M)
(M)
RESISTANCE
RESISTANCE
(OHMS M)
(OHMS M)
(OHMS)
(OHMS)
S/NO.
S/NO.
AB/2 MN/2 AB/2 MN/2
ELECTRO
ELECTRO
SPACING
SPACING
RESISTIVITY
RESISTIVITY
RESISTANCE
RESISTANCE
GEOMETRIC
GEOMETRIC
FACTOR “K”
FACTOR “K”
APPARENT
APPARENT
(OHMS M)
(OHMS M)
(M)
(M)
(OHMS)
(OHMS)
DE
DE
S/NO.
S/NO.
AB/2 MN/2 AB/2 MN/2
1 2 0.5 11.78 25.13 296.03 1 2 0.5 11.78 27.61 325.25
2 3 0.5 27.5 15.46 425.15 2 3 0.5 27.5 15.11 415.53
3 6 0.5 112.3 6.990 784.98 3 6 0.5 112.3 5.596 628.43
4 9 0.5 254 4.232 1074.93 4 9 0.5 254 3.296 837.18
5 9 2.0 60.5 22.06 1334.63 5 9 2.0 60.5 14.75 892.38
6 15 2.0 173.6 10.68 1854.10 6 15 2.0 173.6 7.142 1239.85
7 25 2.0 488 5.935 2896.28 7 25 2.0 488 3.398 1658.22
8 40 2.0 1253 3.768 4721.30 8 40 2.0 1253 1.994 2498.48
9 50 2.0 1960 2.771 5431.16 9 50 2.0 1960 1.353 2651.88
10 75 2.0 4410 1.531 6751.71 10 75 2.0 4410 0.6573 2898.69
11 75 10.0 868 6.429 5580.37 11 75 10.0 868 3.296 2860.93
12 100 10.0 1555 4.232 6580.76 12 100 10.0 1555 1.750 2721.25
13 150 10.0 3520 1.058 3724.16 13 150 10.0 3520 0.761 2678.72
14 200 10.0 6270 0.511 3203.97 14 200 10.0 6270 0.4772 2992.04
15 300 10.0 14120 0.233 3289.96 15 300 10.0 14120 0.1658 2341.10
16 300 20.0 7040 0.411 2893.44 16 300 20.0 7040
17 400 20.0 12530 0.237 2969.61 17 400 20.0 12530
LOCATION:OCHANJA LOCATION: ALLOMO
TOWN: OCHANJA TOWN: ALLOMO
L.G.A: OFU L.G.A: ALLOMO VES 18
STATE: KOGI STATE: KOGI
LONGITUDE: E007007'. 047'' LONGITUDE: E0070 07'. 526''
LATITUDE: N07057'. 809'' LATITUDE: N070 57'. 440''
ELEVATION: 300 m ELEVATION: 311 m
GEOMETRIC FACTOR
GEOMETRIC FACTOR
RESISTIVITY (OHMS
RESISTIVITY (OHMS
ELECTRODE
ELECTRODE
SPACING
SPACING
RESISTANCE
RESISTANCE
APPARENT
APPARENT
(M)
(M)
(OHMS)
(OHMS)
S/NO.
S/NO.
“K”
“K”
M)
M)
AB/2 MN/2 AB/2 MN/2
ELECTRODE
ELECTRODE
SPACING
SPACING
RESISTIVITY
RESISTIVITY
RESISTANCE
RESISTANCE
GEOMETRIC
GEOMETRIC
FACTOR “K”
FACTOR “K”
APPARENT
APPARENT
(OHMS M)
(OHMS M)
(M)
(M)
(OHMS)
(OHMS)
S/NO.
S/NO.
Sounding Curves