# Efficacy of Hummel (Modified Schlumberger) Arrays of Vertical Electrical Sounding in Groundwater Exploration: Case Study of Parts of Ibadan Metropolis, Southwestern Nigeria

Authors: Micheal Oladunjoye; Solomon Jekayinfa
Journal: International Journal of Geophysics
 (2015)

---

## Abstract

This research compared the interpretation results of the Vertical Electrical Sounding data acquired using the conventional Schlumberger and modified Schlumberger arrays with a view to assessing the effectiveness of the modified Schlumberger arrays of vertical electrical sounding as an alternative to the conventional Schlumberger array at sites with space constraint during groundwater exploration. A total of thirty-seven (37) sounding locations were occupied and one hundred (100) sounding data for both conventional Schlumberger and modified Schlumberger arrays were collected across different rock units within Ibadan metropolis, south-western Nigeria, with electrode spacing (A
B
/
2) ranging from 1 to 75 m. The field data were interpreted qualitatively by curve matching and computer iterative methods. Also, statistical analysis of subsurface units and the coefficient of correlation “R” of the statistical plots of the field data shows the relationship between the different arrays. The raw data plot of the different arrays shows significant similarities while statistical analysis of the geo-electric parameters obtained from the different arrays across varied lithologic units show that strong relationships exist between the different field methods. The coefficient of correlation R with values ranging from 0.7 to 0.99 implies that a good similarity exists between the different field methods employed in this study. Hence, modified Schlumberger arrays can be said to be a good alternative to the conventional Schlumberger array for groundwater exploration especially in urban settings where space constraint is a major challenge.


---

## Body

## 1. Introduction

The development of available groundwater resources is a significant aspect of the provision of potable water in urban and rural settings. For years, various geophysical techniques have been employed either singly or in combination in the search for groundwater. In recent times however, technology advancements have brought about the use of state of the art geophysical techniques in groundwater exploration in complex terrains such as the use of magnetic and DC resistivity surveys in semiarid or arid terrains [1]. Dowsing, an ancient traditional method of exploration, used jointly with the modern slingram electromagnetic survey method and vertical electrical sounding in finding groundwater, revealed more accurate results when the three methods are analysed jointly such that the vulnerability of errors linked to a single method is overcome; thus information from the three methods can be used to pinpoint prospective drilling sites for groundwater with a high degree of accuracy [2]. The use of ground penetrating radar in determining groundwater surfaces by studying the responses of the electromagnetic waves reflected from the water table has been used in recent times in prospecting for groundwater occurrences mostly in sedimentary terrains with low electrically conductive sediments [3–6] and sometimes in crystalline rock terrains [7]. Other modern geophysical techniques for groundwater exploration include the use of airborne electromagnetic inductions [8] and seismic refraction survey in conjunction with vertical electrical sounding (VES) resistivity method [9].The most popular configuration in vertical electrical sounding (VES) is the conventional Schlumberger array, which has a symmetrical layout with electrodes spread on either sides of the array spread [11–17]. However, geophysical investigation especially in highly developed or congested and thickly vegetated areas proved very challenging due to the inaccessibility of sites to conduct soundings of symmetrical configuration of electrodes to the required spacing. This is often responsible for incompleteness of field data. Consequently, the asymmetrical array of electrodes in vertical electrical sounding as described by Nick and Katerina [18] as well as Akintorinwa and Abiola [19] and referred to as modified Schlumberger array (Hummel method) could be employed as an alternative to alleviate problems of space constraint. In electrical resistivity method, there are many arrays which can be used. This includes Schlumberger, modified Schlumberger, Dipole-Dipole, Wenner, and Pole-Pole, among others. Most commonly used among the arrays mentioned above for vertical electrical sounding (VES) is the Schlumberger array, which requires large spacing at both ends for deeper information. The modified arrays enable vertical electrical sounding with the movement of current electrode (C
2) while the other current electrode (C
1) is fixed perpendicular (90°) to the line of spread, while the potential electrode remains temporarily fixed as in the conventional Schlumberger array. This research attempts a comparative study of the conventional Schlumberger and modified Schlumberger arrays in terms of deduced geoelectric parameters, in order to determine the efficacy of the modified Schlumberger arrays for groundwater exploration in developed areas with space constraint as an alternative to the conventional Schlumberger array.


## 2. Location and Physiography of the Study Area

Ibadan lies between longitudes 3°45′ and 4°00′E and latitudes 7°15′ and 7°30′N (Figure1(a)) within the basement complex terrain of south western Nigeria and covers about 300 km2. Olayinka and Adeyemi [20] described the scenery as being dominated by three major landforms units, namely, the hills, plains, and river valleys. The hills are the most striking features. Although they constitute less than 20% of the total surface area, two main types can be recognized, namely, quartzite ridges and gneissic inselbergs, of which the former are by the most extensive landform system in the area with an altitude generally ranging from 185 to 222 m above mean sea level [21]. They cover, essentially, the areas between the hill bases and the usually entrenched valley bottoms.(a) Location map of Ibadan showing the study area. (b) Geological map of Ibadan [10].
(a)
(b)


## 3. Geology of the Study Area

Ibadan, a region in the south-western part of Nigeria, is located within the Nigerian basement complex which lies between the West African craton which demarcates the pan-African belt [22]. The Precambrian basement complex of Nigeria may be distinguished into three major successions, a remobilized ancient complex of gneisses and migmatites, a series of approximately N-S trending schist belt, and an intrusive suite of gneisses of pan-African (600 Ma) granites and related rocks. The complex of gneisses and migmatites include minor quartzite and calc-silicate bearing units which bear imprints of Liberian (Ca 2500 Ma) and Eburnean (Ca 2000 Ma) events [23].Ibadan is dominated by rock types such as granite and schist of the metasedimentary scores, banded gneiss and granite gneiss, augen gneiss, and migmatite complex (Figure1(b)). Quartzite (quartz schist) outcrops occur as ridges with relatively high elevation and are commonly schistose in form. Their strike line runs in the north-south direction between 340° and 350° consistently dipping eastwards with characteristic cross-cutting features [24].The study area is underlain by banded and augen gneisses, schistose quartzite (quartz schist), and amphibolite [25]. The study of the weathered rock and soil cover indicates that they are being derived from the weathering of the basement rock which underlain them, the rock is highly weathered and fractured with well-developed foliation planes and strikes between 34° and 352° and dipping between 27° and 37°E [26].


## 4. Hydrogeology of the Study Area

Ibadan the study area is situated within the basement complex of southwestern Nigeria, consisting mainly of Precambrian crystalline rocks. Generally, the basement rocks in their unaltered form are characterized by low porosity and permeability. Therefore the occurrence of groundwater in this region depends on the existence of thick weathered and/or deep fractured zones. In Ibadan area, the weathered profile developed above crystalline basement rocks comprises, from top to bottom, the soil layer, the saprolite (product of in situ chemical weathering of the bedrock), the saprock (fractured bedrock), and the fresh bedrock [27]. The recharge of the aquifers is mainly by the annual rainfall infiltrations and discharges artificially through dug wells and boreholes [28].


## 5. Methodology

Vertical electrical sounding (VES) measures vertical variation of ground resistivity with depth in respect to a fixed station. This is carried out by gradually increasing the interelectrodes spacing at about a fixed centre of array [29]. Conventional Schlumberger array utilizes four electrodes system which are arranged linearly with different interelectrode spacing, the potential electrode remains partially fixed at the center of the spread, and the current electrode is expanded symmetrically about the centre of the spread (Figure 2) as described by Telford et al. [30] and Parasnis [31]. The geometric factor is
(1)
K
s
=
2
π
1
a
-
b
2
-
1
a
+
b
2
-
1
a
+
b
2
+
1
a
-
b
2
,
K
s
=
π
a
2
b
-
b
4
.
With the expression, the electric field is approximately determined; thus
(2)
K
s
=
π
a
2
b
,
ρ
a
=
R
·
K
,
ρ
a
=
π
a
2
b
·
R
.
In this array, the measuring electrodes are placed closely together such that the ratio of the measured voltage to the separation is approximately equal to the voltage gradient at the midpoint of the current spread. However, in the Hummel and “Half-Hummel” (modified Schlumberger) methods of vertical electrical sounding, the potential electrodes remain temporarily fixed while one of the current electrodes is kept fixed throughout the sounding operation at a distance equal to the mobile current electrode spacing “L” (known as the Half-Hummel method) or three times the mobile current electrode, that is, “3
L” (Hummel array) and perpendicular to the line of electrode spread as shown in Figure 3. The geometric factor of the modified Schlumberger arrays K
H was estimated to be twice that of the conventional Schlumberger array K
s [32] that is
(3)
K
H
=
2
K
s
,
ρ
a
=
R
·
K
H
,
where K
s = geometric factor for conventional Schlumberger array, K
H = geometric factor for Hummel and “Half-Hummel” (modified Schlumberger) arrays, ρ
a = apparent resistivity, and R = resistance of subsurface layer.Figure 2
Electrode configuration for conventionalSchlumberger array.Figure 3
Diagrammatic representation of the electrode configuration in Hummel method (modified Schlumberger array). Note that when the distanceO
C
1
=
O
C
2
=
L, we have the “Half-Hummel” method.This research was conducted on different geological units at different stations to ascertain the efficacy of the Hummel and “Half-Hummel” arrays. Conventional Schlumberger array was adopted in conjunction with the Hummel and “Half-Hummel” arrays with the current electrode (A
B
/
2) spacing ranging from 1.0 to 75 m. The fixed current electrode (C
1) was placed at a distance three times the current electrode spacing (i.e., 3
×
A
B
/
2) or equal to the current electrode spacing for the Hummel array and “Half-Hummel” array, respectively, and the results were compared with conventional Schlumberger array at every location of study. Figure 4 shows the application of the Hummel or “Half-Hummel” method in congested and/or developed areas where symmetrical spread of current electrodes is not feasible.Figure 4
Diagrammatic representation of electrode spread in congested area with limited access illustrating Hummel array.One hundred (100) vertical electrical soundings (VES) were carried out in the study area which was fairly distributed across different geological units using both the conventional Schlumberger and modified Schlumberger arrays at every station. The areas investigated include Kumapayi-Olodo community (banded gneiss terrain), Oki community (amphibolite terrain), University of Ibadan campus (augen gneiss and quartz schist terrain), Polytechnic Ibadan Campus (quartz schist terrain), and Ajibode community (augen gneiss terrain). The Geopulse Campus Tigre resistivity meter was used for data collection. Garmin’s global positioning system (GPS) was used to record the geographic coordinates of the VES stations.Field precautions observed during this research exercise include siting of sounding centres where potential electrodes will encounter little surface inhomogeneity, since surface inhomogeneities introduce misleading results in the potential reading by distorting the current flow pattern. Also, good contact was ensured at the current electrodes by driving the electrodes sufficiently deep into the surfaces and proper tightening of current and potential cables to their respective electrodes. Cables were checked properly to avoid current leakages during measurement for good quality data. Good conductors like wire fences, underground telephone lines, and pipes with tendencies to distort the current pattern were avoided by placing the potential electrodes as far as possible away from them; hence most of the readings were taken on open fields and along dirt roads where presence of these objects are minimal. Power lines, especially those accompanied by earth wires, create nuisance when measurements are taken close by as they introduce high voltages in the ground. Such were therefore avoided completely.The apparent resistivity measurements at each VES stations for both the conventional Schlumberger and modified Schlumberger (Hummel and “Half-Hummel”) arrays were plotted against electrode spacing on bilogarithmic graph sheet. The resulting curves were then inspected visually to determine the nature of the subsurface layering delineated by both arrays. Partial curve matching was carried out for the quantitative interpretation of the curves [32, 33]. The results of the curve matching (layer resistivity and thickness) were fed into the computer as starting model parameter in an iterative forward modelling technique using WINRESIST computer software. With the advent of computing facilities the inversion process has become widely used in geophysical data interpretation. Inversion is an iterative process that is controlled by a mathematical equation [34]. A model is proposed from the field data and fed into a computer and an inversion programme is run; the result of which is a new model whose parameters are compared to the field data. These parameters are altered and again run in the computer. This is repeated until a model is found whose calculated data best matches the field data. From the inversion results (layer resistivity and thickness), layer parameter charts were generated. Statistical calculation (test of significance) and the coefficient of correlation “R” of the statistical-plots of the layer parameters were carried out in order to estimate the relationship between the different arrays used.For the test of significance, both the raw data and layer parameters using the various means and variance were subjected to statistical treatment at 5% level of significance as follows.(i)
Null hypothesis (U
NH): there is a significant relationship between the parameters.
(ii)
Alternative hypothesis (U
T): U
T is 1.96 at 5% level of significance from standard normal variant table.
(iii)
Level of significance: 5%.
(iv)
Criterion: accept null hypothesis ifU
NH
<
U
T (significance relationship exists between parameters); reject if U
NH
>
U
T (i.e., no significance relationship exists between parameters).The standard normal variate is
(4)
U
NH
=
x
1
-
x
2
∂
/
n
1
+
∂
/
n
2
,
where x
1 = mean of first parameter, x
2 = mean of second parameter, ∂
1 = variance of first parameter, ∂
2 = variance of second parameter, n
1 = sample size of first parameter, n
2 = sample size of second parameter, U
NH
1 = null hypothesis test for conventional Schlumberger array versus Hummel array, and U
NH
2 = null hypothesis test for conventional Schlumberger array versus “Half-Hummel” array.


## 6. Results and Discussion

### 6.1. Presentation of Results

The results are presented as field curves showing layer parameters, tables, histograms showing test of hypothesis results, and statistical correlation plots. The summary of the VES interpretation results for the conventional Schlumberger and modified Schlumberger (Hummel and Half-Hummel) at VES stations occupied across different geological units within the study area are presented in Table1. Three to five layers were delineated from inverted curves generated from both arrays. The same curve types were identified by both arrays at almost all the sounding stations (Table 1).Table 1
Summary of geoelectric parameters obtained in the study area.
Location
Rock type
Layer resistivity for Schlumberger
Layer resistivity for Hummels
Layer resistivity for Half-Hummels
Thickness (m) for Schlumberger
Thickness (m) for Hummels
Thickness (m) for Half-Hummels
Depth (m) for Schlumberger
Depth (m) for Hummels
Depth (m) for Half-Hummels
Probable layer for Schlumberger
Probable layer for Hummels
Probable layer for Half-Hummels
Curve type for Schlumberger
Curve type for Hummel
Curve type for Half-Hummels
Aquifer type
Kumapayi
Banded gneiss
1
235
217
228
0.6
0.8
0.8
0.6
0.8
0.8
Top soil
Top soil
Top soil
H
H
H
WB
66
42
49
4.3
2.2
2.8
4.9
2.9
3.6
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
972
934
1183
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
Fresh bedrock
—
—
—
2
613
486
595
0.5
0.5
0.5
0.5
0.5
0.5
Top soil
Top soil
Top soil
H
H
H
WB
32
31
25
4
3.3
2.6
4.5
3.8
3.1
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
745
885
835
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
Fresh bedrock
—
—
—
3
63
36
250
0.9
1.6
0.5
0.9
1.6
0.5
Top soil
Top soil
Top soil
H
H
H
WB
32
21
16
5.9
2.7
3.5
6.8
4.3
3.9
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
1511
1391
1143
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
Fresh bedrock
—
—
—
4
73
67
53
0.7
0.6
0.7
0.7
0.6
0.7
Top soil
Top soil
Top soil
H
H
H
WB
25
21
20
2.1
2
1.9
2.8
2.7
2.8
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
982
1019
1338
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
Fresh bedrock
—
—
—
5
250
297
172
1
0.9
0.7
1
0.9
0.7
Top soil
Top soil
Top soil
H
H
H
WB
29
31
36
3.8
4
4
4.8
4.9
4.7
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
1045
1500
1807
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
Fresh bedrock
—
—
—
Alaagbon
Banded gneiss
6
61
53
1
1.3
—
1
1.3
—
Top soil
Top soil
—
H
H
—
WB
23
22
—
2.7
3.1
—
3.7
4.4
—
(WB) clayey
(WB) clayey
—
—
—
—
2600
2094
—
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
—
—
—
—
7
79
70
—
0.6
0.7
—
0.6
0.7
—
Top soil
Top soil
—
H
H
—
WB
23
20
—
1.8
1.6
—
2.4
2.3
—
(WB) clayey
(WB) clayey
—
—
—
—
3277
3210
—
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
—
—
—
—
8
187
156
—
1
1.1
—
1
1.1
—
Top soil
Top soil
—
H
H
—
WB
29
15
—
4.7
2.3
—
5.7
3.4
—
(WB) clayey
(WB) clayey
—
—
—
—
1236
1125
—
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
—
—
—
—
Ogoluwa
Banded gneiss
9
286
321
—
1.7
1.6
—
1.7
1.6
—
Top soil
Top soil
—
H
H
—
WB
31
31
—
3.9
3.2
—
5.6
4.8
—
(WB) clayey
(WB) clayey
—
—
—
—
474
524
—
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
—
—
—
—
10
149
200
—
0.9
0.5
—
0.9
0.5
—
Top soil
Top soil
—
H
H
—
WB
39
42
—
4.4
5.1
—
5.3
5.6
—
(WB) clayey
(WB) clayey
—
—
—
—
486
411
—
—
—
—
—
—
—
Fresh bedrock
Fractured bedrock
—
—
—
—
Oki
Amphibolite
11
23
19
18
1.5
1.5
1.3
1.5
1.5
1.3
Top soil
Top soil
Top soil
H
H
H
WB
14
15
13
3.7
4.1
3
5.1
5.6
4.3
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
265
262
595
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
Fresh bedrock
—
—
—
12
25
24
1.1
1.4
1.1
1.4
Top soil
Top soil
H
H
WFB
13
12
5.9
4.8
7
6.1
(WB) clayey
(WB) clayey
—
—
—
71
133
—
—
—
—
—
—
Fractured bedrock
Fractured bedrock
—
—
—
13
29
18
0.6
0.6
0.6
0.6
Top soil
Top soil
H
H
WB
18
12
2.4
1.8
3
2.4
(WB) clayey
(WB) clayey
—
—
—
359
483
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
—
—
—
Jagun Alekuso
Amphibolite
14
40
43
53
0.7
1.1
1
0.7
1.1
1
Top soil
Top soil
Top soil
H
H
A
WFB
35
40
65
2.1
1.7
5
2.8
2.7
6
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
368
309
583
—
—
—
—
—
—
Fractured bedrock
Fractured bedrock
Fractured bedrock
—
—
—
15
59
48
71
0.8
0.8
1
0.8
0.8
1
Top soil
Top soil
Top soil
H
H
H
WFB
18
21
17
6.1
5.9
6.8
6.9
6.7
7.9
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
110
62
200
—
—
—
—
—
—
Fractured bedrock
Fractured bedrock
Fractured bedrock
—
—
—
Oki
Amphibolite
—
—
—
16
125
125
138
1.1
1.2
1.4
1.1
1.2
1.4
Top soil
Top soil
Top soil
H
H
H
WFB
53
51
44
7.4
8.6
6.4
8.5
9.9
7.8
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
85
81
66
—
—
—
—
—
—
Fractured bedrock
Fractured bedrock
Fractured bedrock
—
—
—
17
52
43
—
0.7
0.7
0.7
0.7
Top soil
Top soil
—
H
H
—
WB
31
33
—
28.6
23
29.4
23.7
(WB) clayey
(WB) clayey
—
—
—
—
683
652
—
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
—
—
—
—
Gelato
Amphibolite
18
62
46
—
0.8
0.8
—
0.8
0.8
—
Top soil
Top soil
—
KH
KH
KH
WB
107
117
—
3.6
4.5
—
4.3
5.3
—
Lateritic clay
Lateritic clay
—
—
—
—
40
27
—
14.1
9.3
—
18.5
14.6
—
(WB) clayey
(WB) clayey
—
—
—
—
490
462
—
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
—
—
—
—
Oki
Amphibolite
19
20
17
17
0.8
1
0.6
0.8
1
0.6
Top soil
Top soil
Top soil
H
H
H
WB
13
12
15
2.6
2.2
3.1
3.4
3.3
3.8
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
458
264
298
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
Fresh bedrock
—
—
—
Monatan
Amphibolite
20
29
26
—
1
1.3
—
1
1.3
—
Top soil
Top soil
—
H
H
—
WB
27
22
—
1.4
0.4
—
2.4
1.7
—
(WB) clayey
(WB) clayey
—
—
—
—
318
367
—
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
—
—
—
—
Ajibode/UI
Augen gneiss
21
56
50
46
1.6
1.3
1.1
1.6
1.3
1.1
Top soil
Top soil
Top soil
AA
AA
AA
WB
227
437
428
17.9
18.5
16.3
19.5
19.8
17.4
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic clay
—
—
—
450
561
445
11
19.6
19.7
30.5
39.3
37.1
WB
WB
WB
—
—
—
4720
13333
12388
—
—
—
—
—
—
Fractured bedrock
Fresh bedrock
Fresh bedrock
—
—
—
22
504
311
303
0.6
1
1
0.6
1
1
Top soil
Top soil
Top soil
HA
HA
HA
WB
133
92
89
10
5.7
4.8
10.7
6.6
5.8
Lateritic clay
Lateritic clay
Lateritic clay
—
—
—
158
136
116
6.8
11.3
12.6
17.5
17.9
18.4
WB
WB
WB
—
—
—
1652
2803
2050
—
—
—
—
—
—
Fractured bedrock
Fresh bedrock
Fresh bedrock
—
—
—
23
297
264
277
0.9
0.8
0.7
0.9
0.8
0.7
Top soil
Top soil
Top soil
QH
QH
QH
WFB
175
194
201
3.5
4.4
4.5
4.4
5.2
5.2
Lateritic clay
Lateritic clay
Lateritic clay
—
—
—
33
21
16
11.8
12
12.7
16.2
17.2
18
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
46
46
180
—
—
—
—
—
—
Fractured bedrock
Fractured bedrock
Fractured bedrock
—
—
—
Abadina/UI
Augen gneiss
24
246
111
100
1.7
1.3
1
1.7
1.3
1
Top soil
Top soil
Top soil
K
K
K
748
814
700
9.6
8.6
9.8
10.2
9.9
10.8
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic clay
—
—
—
492
510
552
WB
WB
WB
—
—
—
25
221
137
139
1.6
1.1
1.1
1.6
1.1
1.1
Top soil
Top soil
Top soil
K
K
K
848
408
421
22.1
16.9
14.7
23.7
18
15.9
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic clay
—
—
—
233
222
112
—
—
—
—
—
—
WB
WB
WB
—
—
—
Botanical Garden/UI
Augen gneiss
26
420
364
324
0.9
0.8
0.9
0.9
0.8
0.9
Top soil
Top soil
Top soil
H
H
H
33
33
31
5.7
4.3
4.1
6.6
5.1
5
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
473
387
304
—
—
—
—
—
—
Fresh bedrock
Fractured bedrock
Fractured bedrock
—
—
—
Dairy Unit/UI
Augen gneiss
27
138
132
—
1.6
1.3
—
1.6
1.3
—
Top soil
Top soil
—
H
H
—
WB
24
20
—
15.9
18.1
—
17.5
19.4
—
(WB) clayey
(WB) clayey
—
—
—
—
1044
450
—
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
—
—
—
—
UI Stadium
Quartz schist
28
92
114
118
0.9
0.8
1
0.9
0.8
1
Top soil
Top soil
Top soil
AKH
AKH
AKH
110
120
185
1.3
1.1
1.8
2.2
2
2.7
Sandy clay
Sandy clay
Sandy clay
—
—
—
210
389
432
3.8
2.3
1.8
6
4.3
4.5
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic clay
—
—
—
29
30
28
22.9
28.6
27.7
28.9
32.9
32.2
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
331
377
3153
—
—
—
—
—
—
Fractured bedrock
Fresh bedrock
Fresh bedrock
—
—
—
UI Ijeoma Rd
Quartz schist
29
338
387
368
0.9
0.7
0.6
0.9
0.7
0.6
Top soil
Top soil
Top soil
KH
KH
KH
WFB
596
592
510
1.3
0.6
0.8
2.2
1.3
1.4
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic cay
—
—
—
250
108
102
7.9
7
5.8
10.1
8.4
7.2
WB
WB
WB
—
—
—
703
722
520
—
—
—
—
—
—
Fractured bedrock
Fractured bedrock
Fractured bedrock
—
—
—
Geology Dpt Poly
Quartz schist
30
161
185
179
1
1.1
1
1
1.1
1
Top soil
Top soil
Top soil
KH
KH
KH
WFB
235
207
222
1.6
1.4
1.4
2.6
2.5
2.4
Lateritic clay
Lateritic clay
Lateritic clay
—
—
—
59
76
73
16.1
16.2
15.6
18.7
18.6
18
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
191
170
152
—
—
—
—
—
—
Fractured bedrock
Fractured bedrock
Fractured bedrock
—
—
—
Central Admin Poly
Quartz schist
31
60
61
62
1.2
1.2
1.3
1.2
1.2
1.3
Top soil
Top soil
Top soil
A
A
H
74
62
50
11.8
8.6
6.1
13
9.8
7.3
Sandy clay
Sandy clay
Sandy clay
—
—
—
738
910
864
—
—
—
—
—
—
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic clay
—
—
—
Poly Main Gate
Quartz schist
32
169
183
181
1
1
1
1
1
1
Top soil
Top soil
Top soil
K
K
K
WFB
404
541
573
6.4
5.9
5.6
7.4
6.8
6.6
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic clay
—
—
—
86
87
78
—
—
—
—
—
—
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
Central Lib. Poly
Quartz schist
33
221
206
201
1.1
1.3
1.2
1.1
1.3
1.2
Top soil
Top soil
Top soil
K
K
K
WB
577
507
476
5.9
6.7
6.1
7
8
7.3
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic clay
—
—
—
152
139
96
—
—
—
—
—
—
WB
WB
WB
—
—
—
North Gate Poly
Quartz schist
34
70
57
58
1
1.2
1
1
1.2
1
Top soil
Top soil
Top soil
H
H
A
55
55
62
2.9
1.6
2.9
3.9
2.9
4
Sandy clay
Sandy clay
Sandy clay
—
—
—
400
219
621
—
—
—
—
—
—
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic clay
—
—
—
Olori Annex Poly
Quartz schist
35
206
292
306
0.8
0.7
0.7
0.8
0.7
0.7
Top soil
Top soil
Top soil
HK
HK
HK
WB
60
51
51
2.5
2.6
2.6
3.2
3.3
3.3
Sandy clay
Sandy clay
Sandy clay
—
—
—
197
210
202
9.8
9.8
8.8
13
13.1
12.1
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic clay
—
—
—
95
91
113
—
—
—
—
—
—
WB
WB
WB
—
—
—
Staff Qtrs Poly
Quartz schist
36
200
139
132
0.7
1
0.8
0.7
1
0.8
Top soil
Top soil
Top soil
H
H
H
80
86
93
13.5
11.4
16.4
14.1
12.4
17.2
Sandy clay
Sandy clay
Sandy clay
—
—
—
260
208
994
—
—
—
—
—
—
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic clay
—
—
—
Staff Qtrs Poly
Quartz schist
37
264
232
232
1.3
1.3
1.2
1.3
1.3
1.2
Top soil
Top soil
Top soil
QH
QH
KH
WB
209
169
272
0.8
2.5
1.1
2.1
3.8
2.2
Lateritic clay
Lateritic clay
Lateritic clay
—
—
—
49
27
31
18.7
19.2
23.8
20.8
23
26.1
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
1649
411
326
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
Fractured bedrock
—
—
—
Note:
WB—weathered basement.
WFB—weathered fractured basement.For locations with banded gneiss lithology, the maximum number of layers observed was three. The range of weathered layer resistivity is between 15 and 66 Ωm (clayey), with overburden thickness ranging between 2.3 and 6.8 m; the basement resistivity ranges between 411 and 3277 Ωm for the different array methods employed. For locations with amphibolite lithology, the maximum number of layers observed was four. The range of weathered layer resistivity is between 12 and 65 Ωm (clayey), with overburden thickness ranging between 1.7 and 29.4 m; the basement resistivity ranges between 62 and 683 Ωm for the different array methods employed. Also, locations with augen gneiss lithology have maximum number of four layers. The range of weathered layer resistivity is between 16 and 561 Ωm, with overburden thickness between 5.0 and 39.9 m; the basement resistivity ranges between 46 and 13333 Ωm for the different array methods employed. For locations with quartz schist lithology, the maximum number of layers observed was five. The range of weathered layer resistivity is between 27 and 250 Ωm, with overburden thickness between 7.2 and 32.9 m; the basement resistivity ranges between 152 and 3153 Ωm for the different array methods employed.Similarities were observed between the conventional Schlumberger array, Hummel array, and the Half-Hummel array methods employed as the interpretation obtained from the field data resulted in the same number of subsurface layers. Although, variation occurred in the nature of the bedrock (fractured or fresh) interpreted for conventional Schlumberger and modified Schlumberger array methods in locations 10, 21, 22, and 26 (Table1), this may be due to structural variations and differential weathering within the bedrock. However, overall similarities can be established between the results of data obtained using the conventional Schlumberger array, Hummel array, and the Half-Hummel array methods at all vertical electrical sounding stations within the area of study. The significant variation in overburden thickness (depth to basement), recorded for the conventional Schlumberger array and modified Schlumberger array methods in some VES locations such as 3, 8, 17, 18, 21, 28, 29, and 37 (Table 1), was largely due to differential weathering of the bedrock. Also, the significant variation in the values of basement resistivity recorded for the conventional Schlumberger array and modified Schlumberger array methods at some VES locations such as 21, 22, 27, and 37 (Table 1) may be due to differential weathering within the bedrock as well as varying degree of saturation of the bedrock unit from one point to another.


### 6.2. Comparison of Conventional and Modified Schlumberger Arrays Results

The results obtained from the different array methods at each VES stations within the study area were compared visually, through the use of statistical correlation plots and test of hypothesis.Figures5(a) and 5(b) show the similarities in the results of interpretation of the field data obtained from two different VES stations across quartz schist (locations 28 and 31). Figure 6 shows that the topsoil resistivity and thickness obtained from the conventional Schlumberger and modified Schlumberger arrays were virtually the same for VES stations on quartz schist while the layer parameters significantly correlate for the VES stations across the same lithology (Figure 7). Also, Figure 8 shows that the basement resistivity values from the different arrays are very similar for the quartz schist terrain. Figures 9 and 10 represent the statistical correlation plots of the apparent resistivity data obtained at the same stations as in Figures 5(a) and 5(b), respectively. The plots show a good linear relationship between apparent resistivity values determined from the different arrays with coefficient of correlation (R) for the Hummel/conventional Schlumberger and Half-Hummel/conventional Schlumberger arrays ranging between 0.98 and 0.99.Comparison of results obtained from the interpretation of field data collected using the different array methods at VES stations on quartz schist for (a) location 28, University of Ibadan, and (b) location 31, the Polytechnic Ibadan.
(a)
(b)Bar chart and statistical test of hypothesis table for the different array methods comparing (a) topsoil thickness and (b) topsoil resistivity, for locations on quartz schist.
(a)
(b)Figure 7
Bar chart and statistical test of hypothesis table for the different array methods comparing Overburden thickness for locations on quartz schist.Figure 8
Bar chart and statistical test of hypothesis table for the different array methods comparing basement resistivity for locations on quartz schist.Statistical plot of apparent resistivity data for location 28 (quartz schist) (a) Hummel against conventional Schlumberger (VES) and (b) “Half-Hummel” against conventional Schlumberger (VES).
(a)
(b)Statistical plot of apparent resistivity data for location 31 (quartz schist) (a) Hummel against conventional Schlumberger (VES) and (b) “Half-Hummel” against conventional Schlumberger (VES).
(a)
(b)Figures11(a) and 11(b) show the similarities in the results of interpretation of the field data obtained from two different VES stations across banded gneiss (locations 2 and 4). Figure 12 shows that the topsoil resistivity and thickness determined from conventional and modified Schlumberger arrays were virtually the same for the VES stations on banded gneiss. The layer parameters obtained show some degree of correlation for the VES stations across the same lithology (Figure 13). Also, Figure 14 shows that the basement resistivity values from the different arrays are very similar for the banded gneiss terrain. Figures 15 and 16 show the statistical correlation plots of the apparent resistivity data obtained at the same stations. The plots show a good linear relationship between apparent resistivity values determined from the different arrays with coefficient of correlation (R) for the Hummel/conventional Schlumberger and Half-Hummel/conventional Schlumberger arrays ranging between 0.96 and 0.99.Comparison of results obtained from the interpretation of field data collected using the different array methods at VES stations on banded gneiss for (a) location 2, Olodo Ibadan, and (b) location 4, Olodo Ibadan.
(a)
(b)Bar chart and statistical test of hypothesis table for the different array methods comparing (a) topsoil thickness and (b) topsoil resistivity, for stations on banded gneiss.
(a)
(b)Figure 13
Bar chart and statistical test of hypothesis table for the different array methods comparing overburden thickness for locations on banded gneiss.Figure 14
Bar chart and statistical test of hypothesis table for the different array methods comparing basement resistivity for locations on banded gneiss.Statistical plot of apparent resistivity data for location 2 (banded gneiss) (a) Hummel against conventional Schlumberger (VES) and (b) “Half-Hummel” against conventional Schlumberger (VES).
(a)
(b)Statistical plot of apparent resistivity data for location 4 (banded gneiss) (a) Hummel against conventional Schlumberger (VES) and (b) “Half-Hummel” against conventional Schlumberger (VES).
(a)
(b)The similarities in the results of interpretation of the field data obtained from two different VES stations (locations 22 and 24) across augen gneiss terrain using the different array methods are as shown in Figures17(a) and 17(b). Figure 18 shows that the topsoil resistivity and thickness determined from conventional Schlumberger and modified Schlumberger arrays were virtually the same for VES stations on augen gneiss. The layer parameters obtained show some degree of correlation for the VES stations across the same lithology (Figure 19). Also, Figure 20 shows that the basement resistivity values for the different arrays are very similar for the augen gneiss terrain. Figures 21 and 22 are the statistical correlation plots of the apparent resistivity data obtained at the same stations. The plots show a good linear relationship between apparent resistivity values obtained from the different arrays with coefficient of correlation (R) for the Hummel/conventional Schlumberger and Half-Hummel/conventional Schlumberger arrays ranging between 0.95 and 0.99.Comparison of results obtained from the interpretation of field data collected using the different array methods at VES stations on augen gneiss for (a) location 22, Ajibode Ibadan, and (b) location 24, Ajibode Ibadan.
(a)
(b)Bar chart and statistical test of hypothesis table for the different array methods comparing (a) topsoil thickness and (b) topsoil resistivity, for stations on augen gneiss.
(a)
(b)Figure 19
Bar chart and statistical test of hypothesis table for the different array methods comparing overburden thickness for locations on augen gneiss.Figure 20
Bar chart and statistical test of hypothesis table for the different array methods comparing basement resistivity for locations on augen gneiss.Statistical plot of apparent resistivity data for location 22 (augen gneiss) (a) Hummel against conventional Schlumberger (VES) and (b) “Half-Hummel” against conventional Schlumberger (VES).
(a)
(b)Statistical plot of apparent resistivity data for location 24 (augen gneiss) (a) Hummel against conventional Schlumberger (VES) and (b) “Half-Hummel” against conventional Schlumberger (VES).
(a)
(b)As shown in Figures23(a) and 23(b), there are similarities in the results of interpretation of the field data obtained from two different VES stations across amphibolite (locations 16 and 19). Figure 24 shows that the topsoil resistivity and thickness obtained from conventional Schlumberger and modified Schlumberger arrays were virtually the same for VES stations on amphibolite. The layer parameters obtained show some degree of correlation for the VES stations across the same lithology (Figure 25). Also, Figure 26 shows that the basement resistivity values from the different arrays are very similar for the amphibolite terrain. Figures 27 and 28 show the statistical correlation plots of the apparent resistivity data obtained at the same stations. The plots show a good linear relationship between apparent resistivity values obtained from the different arrays with coefficient of correlation (R) for the Hummel/conventional Schlumberger and Half-Hummel/conventional Schlumberger arrays ranging between 0.94 and 0.99.Comparison of results obtained from the interpretation of field data collected using the different array methods at VES stations on amphibolite for (a) location 16, Oki Ibadan, and (b) location 19, Oki Ibadan.
(a)
(b)Bar chart and statistical test of hypothesis table for the different array methods comparing (a) topsoil thickness and (b) topsoil resistivity, for stations on amphibolite.
(a)
(b)Figure 25
Bar chart and statistical test of hypothesis table for the different array methods comparing overburden thickness for locations on amphibolite.Figure 26
Bar chart and statistical test of hypothesis table for the different array methods comparing basement resistivity for locations on amphibolite.Statistical plot of apparent resistivity data for location 16 (amphibolite) (a) Hummel against conventional Schlumberger (VES) and (b) “Half-Hummel” against conventional Schlumberger (VES).
(a)
(b)Statistical plot of apparent resistivity data for location 19 (amphibolite) (a) Hummel against conventional Schlumberger (VES) and (b) “Half-Hummel” against conventional Schlumberger (VES).
(a)
(b)Based on the 1D resistivity and depth models for the different arrays used at each location, it was observed that there is good correlation in many of the locations in terms of respective resistivity, thickness of each layer, and depth to basement (Table1). This therefore suggests that for groundwater resources exploration, any of these arrays can be adopted to effectively characterize electrostratigraphic sections in the subsurface. The degree of correlation (statistical analysis) of the different arrays used was good in most of the locations of study which suggests that the data quality was good and subsurface models obtained from any of the arrays will be the actual representation of the subsurface configuration at such location.Although there is relatively high values of percentage difference obtained from the test of hypothesis for the conventional Schlumberger array and modified Schlumberger array methods across quartz schist and augen gneiss 12 (Figures8 and 20), this has no effect on the quality of correlation between the two methods as the high value was due to high basement resistivity values obtained across these lithologies.


## 6.1. Presentation of Results

The results are presented as field curves showing layer parameters, tables, histograms showing test of hypothesis results, and statistical correlation plots. The summary of the VES interpretation results for the conventional Schlumberger and modified Schlumberger (Hummel and Half-Hummel) at VES stations occupied across different geological units within the study area are presented in Table1. Three to five layers were delineated from inverted curves generated from both arrays. The same curve types were identified by both arrays at almost all the sounding stations (Table 1).Table 1
Summary of geoelectric parameters obtained in the study area.
Location
Rock type
Layer resistivity for Schlumberger
Layer resistivity for Hummels
Layer resistivity for Half-Hummels
Thickness (m) for Schlumberger
Thickness (m) for Hummels
Thickness (m) for Half-Hummels
Depth (m) for Schlumberger
Depth (m) for Hummels
Depth (m) for Half-Hummels
Probable layer for Schlumberger
Probable layer for Hummels
Probable layer for Half-Hummels
Curve type for Schlumberger
Curve type for Hummel
Curve type for Half-Hummels
Aquifer type
Kumapayi
Banded gneiss
1
235
217
228
0.6
0.8
0.8
0.6
0.8
0.8
Top soil
Top soil
Top soil
H
H
H
WB
66
42
49
4.3
2.2
2.8
4.9
2.9
3.6
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
972
934
1183
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
Fresh bedrock
—
—
—
2
613
486
595
0.5
0.5
0.5
0.5
0.5
0.5
Top soil
Top soil
Top soil
H
H
H
WB
32
31
25
4
3.3
2.6
4.5
3.8
3.1
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
745
885
835
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
Fresh bedrock
—
—
—
3
63
36
250
0.9
1.6
0.5
0.9
1.6
0.5
Top soil
Top soil
Top soil
H
H
H
WB
32
21
16
5.9
2.7
3.5
6.8
4.3
3.9
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
1511
1391
1143
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
Fresh bedrock
—
—
—
4
73
67
53
0.7
0.6
0.7
0.7
0.6
0.7
Top soil
Top soil
Top soil
H
H
H
WB
25
21
20
2.1
2
1.9
2.8
2.7
2.8
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
982
1019
1338
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
Fresh bedrock
—
—
—
5
250
297
172
1
0.9
0.7
1
0.9
0.7
Top soil
Top soil
Top soil
H
H
H
WB
29
31
36
3.8
4
4
4.8
4.9
4.7
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
1045
1500
1807
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
Fresh bedrock
—
—
—
Alaagbon
Banded gneiss
6
61
53
1
1.3
—
1
1.3
—
Top soil
Top soil
—
H
H
—
WB
23
22
—
2.7
3.1
—
3.7
4.4
—
(WB) clayey
(WB) clayey
—
—
—
—
2600
2094
—
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
—
—
—
—
7
79
70
—
0.6
0.7
—
0.6
0.7
—
Top soil
Top soil
—
H
H
—
WB
23
20
—
1.8
1.6
—
2.4
2.3
—
(WB) clayey
(WB) clayey
—
—
—
—
3277
3210
—
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
—
—
—
—
8
187
156
—
1
1.1
—
1
1.1
—
Top soil
Top soil
—
H
H
—
WB
29
15
—
4.7
2.3
—
5.7
3.4
—
(WB) clayey
(WB) clayey
—
—
—
—
1236
1125
—
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
—
—
—
—
Ogoluwa
Banded gneiss
9
286
321
—
1.7
1.6
—
1.7
1.6
—
Top soil
Top soil
—
H
H
—
WB
31
31
—
3.9
3.2
—
5.6
4.8
—
(WB) clayey
(WB) clayey
—
—
—
—
474
524
—
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
—
—
—
—
10
149
200
—
0.9
0.5
—
0.9
0.5
—
Top soil
Top soil
—
H
H
—
WB
39
42
—
4.4
5.1
—
5.3
5.6
—
(WB) clayey
(WB) clayey
—
—
—
—
486
411
—
—
—
—
—
—
—
Fresh bedrock
Fractured bedrock
—
—
—
—
Oki
Amphibolite
11
23
19
18
1.5
1.5
1.3
1.5
1.5
1.3
Top soil
Top soil
Top soil
H
H
H
WB
14
15
13
3.7
4.1
3
5.1
5.6
4.3
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
265
262
595
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
Fresh bedrock
—
—
—
12
25
24
1.1
1.4
1.1
1.4
Top soil
Top soil
H
H
WFB
13
12
5.9
4.8
7
6.1
(WB) clayey
(WB) clayey
—
—
—
71
133
—
—
—
—
—
—
Fractured bedrock
Fractured bedrock
—
—
—
13
29
18
0.6
0.6
0.6
0.6
Top soil
Top soil
H
H
WB
18
12
2.4
1.8
3
2.4
(WB) clayey
(WB) clayey
—
—
—
359
483
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
—
—
—
Jagun Alekuso
Amphibolite
14
40
43
53
0.7
1.1
1
0.7
1.1
1
Top soil
Top soil
Top soil
H
H
A
WFB
35
40
65
2.1
1.7
5
2.8
2.7
6
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
368
309
583
—
—
—
—
—
—
Fractured bedrock
Fractured bedrock
Fractured bedrock
—
—
—
15
59
48
71
0.8
0.8
1
0.8
0.8
1
Top soil
Top soil
Top soil
H
H
H
WFB
18
21
17
6.1
5.9
6.8
6.9
6.7
7.9
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
110
62
200
—
—
—
—
—
—
Fractured bedrock
Fractured bedrock
Fractured bedrock
—
—
—
Oki
Amphibolite
—
—
—
16
125
125
138
1.1
1.2
1.4
1.1
1.2
1.4
Top soil
Top soil
Top soil
H
H
H
WFB
53
51
44
7.4
8.6
6.4
8.5
9.9
7.8
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
85
81
66
—
—
—
—
—
—
Fractured bedrock
Fractured bedrock
Fractured bedrock
—
—
—
17
52
43
—
0.7
0.7
0.7
0.7
Top soil
Top soil
—
H
H
—
WB
31
33
—
28.6
23
29.4
23.7
(WB) clayey
(WB) clayey
—
—
—
—
683
652
—
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
—
—
—
—
Gelato
Amphibolite
18
62
46
—
0.8
0.8
—
0.8
0.8
—
Top soil
Top soil
—
KH
KH
KH
WB
107
117
—
3.6
4.5
—
4.3
5.3
—
Lateritic clay
Lateritic clay
—
—
—
—
40
27
—
14.1
9.3
—
18.5
14.6
—
(WB) clayey
(WB) clayey
—
—
—
—
490
462
—
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
—
—
—
—
Oki
Amphibolite
19
20
17
17
0.8
1
0.6
0.8
1
0.6
Top soil
Top soil
Top soil
H
H
H
WB
13
12
15
2.6
2.2
3.1
3.4
3.3
3.8
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
458
264
298
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
Fresh bedrock
—
—
—
Monatan
Amphibolite
20
29
26
—
1
1.3
—
1
1.3
—
Top soil
Top soil
—
H
H
—
WB
27
22
—
1.4
0.4
—
2.4
1.7
—
(WB) clayey
(WB) clayey
—
—
—
—
318
367
—
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
—
—
—
—
Ajibode/UI
Augen gneiss
21
56
50
46
1.6
1.3
1.1
1.6
1.3
1.1
Top soil
Top soil
Top soil
AA
AA
AA
WB
227
437
428
17.9
18.5
16.3
19.5
19.8
17.4
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic clay
—
—
—
450
561
445
11
19.6
19.7
30.5
39.3
37.1
WB
WB
WB
—
—
—
4720
13333
12388
—
—
—
—
—
—
Fractured bedrock
Fresh bedrock
Fresh bedrock
—
—
—
22
504
311
303
0.6
1
1
0.6
1
1
Top soil
Top soil
Top soil
HA
HA
HA
WB
133
92
89
10
5.7
4.8
10.7
6.6
5.8
Lateritic clay
Lateritic clay
Lateritic clay
—
—
—
158
136
116
6.8
11.3
12.6
17.5
17.9
18.4
WB
WB
WB
—
—
—
1652
2803
2050
—
—
—
—
—
—
Fractured bedrock
Fresh bedrock
Fresh bedrock
—
—
—
23
297
264
277
0.9
0.8
0.7
0.9
0.8
0.7
Top soil
Top soil
Top soil
QH
QH
QH
WFB
175
194
201
3.5
4.4
4.5
4.4
5.2
5.2
Lateritic clay
Lateritic clay
Lateritic clay
—
—
—
33
21
16
11.8
12
12.7
16.2
17.2
18
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
46
46
180
—
—
—
—
—
—
Fractured bedrock
Fractured bedrock
Fractured bedrock
—
—
—
Abadina/UI
Augen gneiss
24
246
111
100
1.7
1.3
1
1.7
1.3
1
Top soil
Top soil
Top soil
K
K
K
748
814
700
9.6
8.6
9.8
10.2
9.9
10.8
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic clay
—
—
—
492
510
552
WB
WB
WB
—
—
—
25
221
137
139
1.6
1.1
1.1
1.6
1.1
1.1
Top soil
Top soil
Top soil
K
K
K
848
408
421
22.1
16.9
14.7
23.7
18
15.9
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic clay
—
—
—
233
222
112
—
—
—
—
—
—
WB
WB
WB
—
—
—
Botanical Garden/UI
Augen gneiss
26
420
364
324
0.9
0.8
0.9
0.9
0.8
0.9
Top soil
Top soil
Top soil
H
H
H
33
33
31
5.7
4.3
4.1
6.6
5.1
5
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
473
387
304
—
—
—
—
—
—
Fresh bedrock
Fractured bedrock
Fractured bedrock
—
—
—
Dairy Unit/UI
Augen gneiss
27
138
132
—
1.6
1.3
—
1.6
1.3
—
Top soil
Top soil
—
H
H
—
WB
24
20
—
15.9
18.1
—
17.5
19.4
—
(WB) clayey
(WB) clayey
—
—
—
—
1044
450
—
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
—
—
—
—
UI Stadium
Quartz schist
28
92
114
118
0.9
0.8
1
0.9
0.8
1
Top soil
Top soil
Top soil
AKH
AKH
AKH
110
120
185
1.3
1.1
1.8
2.2
2
2.7
Sandy clay
Sandy clay
Sandy clay
—
—
—
210
389
432
3.8
2.3
1.8
6
4.3
4.5
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic clay
—
—
—
29
30
28
22.9
28.6
27.7
28.9
32.9
32.2
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
331
377
3153
—
—
—
—
—
—
Fractured bedrock
Fresh bedrock
Fresh bedrock
—
—
—
UI Ijeoma Rd
Quartz schist
29
338
387
368
0.9
0.7
0.6
0.9
0.7
0.6
Top soil
Top soil
Top soil
KH
KH
KH
WFB
596
592
510
1.3
0.6
0.8
2.2
1.3
1.4
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic cay
—
—
—
250
108
102
7.9
7
5.8
10.1
8.4
7.2
WB
WB
WB
—
—
—
703
722
520
—
—
—
—
—
—
Fractured bedrock
Fractured bedrock
Fractured bedrock
—
—
—
Geology Dpt Poly
Quartz schist
30
161
185
179
1
1.1
1
1
1.1
1
Top soil
Top soil
Top soil
KH
KH
KH
WFB
235
207
222
1.6
1.4
1.4
2.6
2.5
2.4
Lateritic clay
Lateritic clay
Lateritic clay
—
—
—
59
76
73
16.1
16.2
15.6
18.7
18.6
18
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
191
170
152
—
—
—
—
—
—
Fractured bedrock
Fractured bedrock
Fractured bedrock
—
—
—
Central Admin Poly
Quartz schist
31
60
61
62
1.2
1.2
1.3
1.2
1.2
1.3
Top soil
Top soil
Top soil
A
A
H
74
62
50
11.8
8.6
6.1
13
9.8
7.3
Sandy clay
Sandy clay
Sandy clay
—
—
—
738
910
864
—
—
—
—
—
—
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic clay
—
—
—
Poly Main Gate
Quartz schist
32
169
183
181
1
1
1
1
1
1
Top soil
Top soil
Top soil
K
K
K
WFB
404
541
573
6.4
5.9
5.6
7.4
6.8
6.6
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic clay
—
—
—
86
87
78
—
—
—
—
—
—
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
Central Lib. Poly
Quartz schist
33
221
206
201
1.1
1.3
1.2
1.1
1.3
1.2
Top soil
Top soil
Top soil
K
K
K
WB
577
507
476
5.9
6.7
6.1
7
8
7.3
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic clay
—
—
—
152
139
96
—
—
—
—
—
—
WB
WB
WB
—
—
—
North Gate Poly
Quartz schist
34
70
57
58
1
1.2
1
1
1.2
1
Top soil
Top soil
Top soil
H
H
A
55
55
62
2.9
1.6
2.9
3.9
2.9
4
Sandy clay
Sandy clay
Sandy clay
—
—
—
400
219
621
—
—
—
—
—
—
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic clay
—
—
—
Olori Annex Poly
Quartz schist
35
206
292
306
0.8
0.7
0.7
0.8
0.7
0.7
Top soil
Top soil
Top soil
HK
HK
HK
WB
60
51
51
2.5
2.6
2.6
3.2
3.3
3.3
Sandy clay
Sandy clay
Sandy clay
—
—
—
197
210
202
9.8
9.8
8.8
13
13.1
12.1
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic clay
—
—
—
95
91
113
—
—
—
—
—
—
WB
WB
WB
—
—
—
Staff Qtrs Poly
Quartz schist
36
200
139
132
0.7
1
0.8
0.7
1
0.8
Top soil
Top soil
Top soil
H
H
H
80
86
93
13.5
11.4
16.4
14.1
12.4
17.2
Sandy clay
Sandy clay
Sandy clay
—
—
—
260
208
994
—
—
—
—
—
—
Compacted lateritic clay
Compacted lateritic clay
Compacted lateritic clay
—
—
—
Staff Qtrs Poly
Quartz schist
37
264
232
232
1.3
1.3
1.2
1.3
1.3
1.2
Top soil
Top soil
Top soil
QH
QH
KH
WB
209
169
272
0.8
2.5
1.1
2.1
3.8
2.2
Lateritic clay
Lateritic clay
Lateritic clay
—
—
—
49
27
31
18.7
19.2
23.8
20.8
23
26.1
(WB) clayey
(WB) clayey
(WB) clayey
—
—
—
1649
411
326
—
—
—
—
—
—
Fresh bedrock
Fresh bedrock
Fractured bedrock
—
—
—
Note:
WB—weathered basement.
WFB—weathered fractured basement.For locations with banded gneiss lithology, the maximum number of layers observed was three. The range of weathered layer resistivity is between 15 and 66 Ωm (clayey), with overburden thickness ranging between 2.3 and 6.8 m; the basement resistivity ranges between 411 and 3277 Ωm for the different array methods employed. For locations with amphibolite lithology, the maximum number of layers observed was four. The range of weathered layer resistivity is between 12 and 65 Ωm (clayey), with overburden thickness ranging between 1.7 and 29.4 m; the basement resistivity ranges between 62 and 683 Ωm for the different array methods employed. Also, locations with augen gneiss lithology have maximum number of four layers. The range of weathered layer resistivity is between 16 and 561 Ωm, with overburden thickness between 5.0 and 39.9 m; the basement resistivity ranges between 46 and 13333 Ωm for the different array methods employed. For locations with quartz schist lithology, the maximum number of layers observed was five. The range of weathered layer resistivity is between 27 and 250 Ωm, with overburden thickness between 7.2 and 32.9 m; the basement resistivity ranges between 152 and 3153 Ωm for the different array methods employed.Similarities were observed between the conventional Schlumberger array, Hummel array, and the Half-Hummel array methods employed as the interpretation obtained from the field data resulted in the same number of subsurface layers. Although, variation occurred in the nature of the bedrock (fractured or fresh) interpreted for conventional Schlumberger and modified Schlumberger array methods in locations 10, 21, 22, and 26 (Table1), this may be due to structural variations and differential weathering within the bedrock. However, overall similarities can be established between the results of data obtained using the conventional Schlumberger array, Hummel array, and the Half-Hummel array methods at all vertical electrical sounding stations within the area of study. The significant variation in overburden thickness (depth to basement), recorded for the conventional Schlumberger array and modified Schlumberger array methods in some VES locations such as 3, 8, 17, 18, 21, 28, 29, and 37 (Table 1), was largely due to differential weathering of the bedrock. Also, the significant variation in the values of basement resistivity recorded for the conventional Schlumberger array and modified Schlumberger array methods at some VES locations such as 21, 22, 27, and 37 (Table 1) may be due to differential weathering within the bedrock as well as varying degree of saturation of the bedrock unit from one point to another.


## 6.2. Comparison of Conventional and Modified Schlumberger Arrays Results

The results obtained from the different array methods at each VES stations within the study area were compared visually, through the use of statistical correlation plots and test of hypothesis.Figures5(a) and 5(b) show the similarities in the results of interpretation of the field data obtained from two different VES stations across quartz schist (locations 28 and 31). Figure 6 shows that the topsoil resistivity and thickness obtained from the conventional Schlumberger and modified Schlumberger arrays were virtually the same for VES stations on quartz schist while the layer parameters significantly correlate for the VES stations across the same lithology (Figure 7). Also, Figure 8 shows that the basement resistivity values from the different arrays are very similar for the quartz schist terrain. Figures 9 and 10 represent the statistical correlation plots of the apparent resistivity data obtained at the same stations as in Figures 5(a) and 5(b), respectively. The plots show a good linear relationship between apparent resistivity values determined from the different arrays with coefficient of correlation (R) for the Hummel/conventional Schlumberger and Half-Hummel/conventional Schlumberger arrays ranging between 0.98 and 0.99.Comparison of results obtained from the interpretation of field data collected using the different array methods at VES stations on quartz schist for (a) location 28, University of Ibadan, and (b) location 31, the Polytechnic Ibadan.
(a)
(b)Bar chart and statistical test of hypothesis table for the different array methods comparing (a) topsoil thickness and (b) topsoil resistivity, for locations on quartz schist.
(a)
(b)Figure 7
Bar chart and statistical test of hypothesis table for the different array methods comparing Overburden thickness for locations on quartz schist.Figure 8
Bar chart and statistical test of hypothesis table for the different array methods comparing basement resistivity for locations on quartz schist.Statistical plot of apparent resistivity data for location 28 (quartz schist) (a) Hummel against conventional Schlumberger (VES) and (b) “Half-Hummel” against conventional Schlumberger (VES).
(a)
(b)Statistical plot of apparent resistivity data for location 31 (quartz schist) (a) Hummel against conventional Schlumberger (VES) and (b) “Half-Hummel” against conventional Schlumberger (VES).
(a)
(b)Figures11(a) and 11(b) show the similarities in the results of interpretation of the field data obtained from two different VES stations across banded gneiss (locations 2 and 4). Figure 12 shows that the topsoil resistivity and thickness determined from conventional and modified Schlumberger arrays were virtually the same for the VES stations on banded gneiss. The layer parameters obtained show some degree of correlation for the VES stations across the same lithology (Figure 13). Also, Figure 14 shows that the basement resistivity values from the different arrays are very similar for the banded gneiss terrain. Figures 15 and 16 show the statistical correlation plots of the apparent resistivity data obtained at the same stations. The plots show a good linear relationship between apparent resistivity values determined from the different arrays with coefficient of correlation (R) for the Hummel/conventional Schlumberger and Half-Hummel/conventional Schlumberger arrays ranging between 0.96 and 0.99.Comparison of results obtained from the interpretation of field data collected using the different array methods at VES stations on banded gneiss for (a) location 2, Olodo Ibadan, and (b) location 4, Olodo Ibadan.
(a)
(b)Bar chart and statistical test of hypothesis table for the different array methods comparing (a) topsoil thickness and (b) topsoil resistivity, for stations on banded gneiss.
(a)
(b)Figure 13
Bar chart and statistical test of hypothesis table for the different array methods comparing overburden thickness for locations on banded gneiss.Figure 14
Bar chart and statistical test of hypothesis table for the different array methods comparing basement resistivity for locations on banded gneiss.Statistical plot of apparent resistivity data for location 2 (banded gneiss) (a) Hummel against conventional Schlumberger (VES) and (b) “Half-Hummel” against conventional Schlumberger (VES).
(a)
(b)Statistical plot of apparent resistivity data for location 4 (banded gneiss) (a) Hummel against conventional Schlumberger (VES) and (b) “Half-Hummel” against conventional Schlumberger (VES).
(a)
(b)The similarities in the results of interpretation of the field data obtained from two different VES stations (locations 22 and 24) across augen gneiss terrain using the different array methods are as shown in Figures17(a) and 17(b). Figure 18 shows that the topsoil resistivity and thickness determined from conventional Schlumberger and modified Schlumberger arrays were virtually the same for VES stations on augen gneiss. The layer parameters obtained show some degree of correlation for the VES stations across the same lithology (Figure 19). Also, Figure 20 shows that the basement resistivity values for the different arrays are very similar for the augen gneiss terrain. Figures 21 and 22 are the statistical correlation plots of the apparent resistivity data obtained at the same stations. The plots show a good linear relationship between apparent resistivity values obtained from the different arrays with coefficient of correlation (R) for the Hummel/conventional Schlumberger and Half-Hummel/conventional Schlumberger arrays ranging between 0.95 and 0.99.Comparison of results obtained from the interpretation of field data collected using the different array methods at VES stations on augen gneiss for (a) location 22, Ajibode Ibadan, and (b) location 24, Ajibode Ibadan.
(a)
(b)Bar chart and statistical test of hypothesis table for the different array methods comparing (a) topsoil thickness and (b) topsoil resistivity, for stations on augen gneiss.
(a)
(b)Figure 19
Bar chart and statistical test of hypothesis table for the different array methods comparing overburden thickness for locations on augen gneiss.Figure 20
Bar chart and statistical test of hypothesis table for the different array methods comparing basement resistivity for locations on augen gneiss.Statistical plot of apparent resistivity data for location 22 (augen gneiss) (a) Hummel against conventional Schlumberger (VES) and (b) “Half-Hummel” against conventional Schlumberger (VES).
(a)
(b)Statistical plot of apparent resistivity data for location 24 (augen gneiss) (a) Hummel against conventional Schlumberger (VES) and (b) “Half-Hummel” against conventional Schlumberger (VES).
(a)
(b)As shown in Figures23(a) and 23(b), there are similarities in the results of interpretation of the field data obtained from two different VES stations across amphibolite (locations 16 and 19). Figure 24 shows that the topsoil resistivity and thickness obtained from conventional Schlumberger and modified Schlumberger arrays were virtually the same for VES stations on amphibolite. The layer parameters obtained show some degree of correlation for the VES stations across the same lithology (Figure 25). Also, Figure 26 shows that the basement resistivity values from the different arrays are very similar for the amphibolite terrain. Figures 27 and 28 show the statistical correlation plots of the apparent resistivity data obtained at the same stations. The plots show a good linear relationship between apparent resistivity values obtained from the different arrays with coefficient of correlation (R) for the Hummel/conventional Schlumberger and Half-Hummel/conventional Schlumberger arrays ranging between 0.94 and 0.99.Comparison of results obtained from the interpretation of field data collected using the different array methods at VES stations on amphibolite for (a) location 16, Oki Ibadan, and (b) location 19, Oki Ibadan.
(a)
(b)Bar chart and statistical test of hypothesis table for the different array methods comparing (a) topsoil thickness and (b) topsoil resistivity, for stations on amphibolite.
(a)
(b)Figure 25
Bar chart and statistical test of hypothesis table for the different array methods comparing overburden thickness for locations on amphibolite.Figure 26
Bar chart and statistical test of hypothesis table for the different array methods comparing basement resistivity for locations on amphibolite.Statistical plot of apparent resistivity data for location 16 (amphibolite) (a) Hummel against conventional Schlumberger (VES) and (b) “Half-Hummel” against conventional Schlumberger (VES).
(a)
(b)Statistical plot of apparent resistivity data for location 19 (amphibolite) (a) Hummel against conventional Schlumberger (VES) and (b) “Half-Hummel” against conventional Schlumberger (VES).
(a)
(b)Based on the 1D resistivity and depth models for the different arrays used at each location, it was observed that there is good correlation in many of the locations in terms of respective resistivity, thickness of each layer, and depth to basement (Table1). This therefore suggests that for groundwater resources exploration, any of these arrays can be adopted to effectively characterize electrostratigraphic sections in the subsurface. The degree of correlation (statistical analysis) of the different arrays used was good in most of the locations of study which suggests that the data quality was good and subsurface models obtained from any of the arrays will be the actual representation of the subsurface configuration at such location.Although there is relatively high values of percentage difference obtained from the test of hypothesis for the conventional Schlumberger array and modified Schlumberger array methods across quartz schist and augen gneiss 12 (Figures8 and 20), this has no effect on the quality of correlation between the two methods as the high value was due to high basement resistivity values obtained across these lithologies.


## 7. Conclusions

Geoelectrical survey at some locations within Ibadan metropolis underlain by Precambrian basement complex rocks has been undertaken using the conventional Schlumberger, Hummel, and “Half-Hummel” arrays at each point of investigation. The field data and curves obtained were interpreted quantitatively by curve matching technique and later subjected to computer assisted iterative interpretation. Various subsurface geologic units; topsoil, clay, weathered/fractured and fresh bedrock were identified. The results obtained from the different arrays were correlated with each other, that is, the field curves and statistical plots. The statistical plots of the apparent resistivities obtained at each VES station were generated and the coefficient of correlation “R” for the relationship between the different arrays was also determined.The correlation values obtained on the different geological units indicated that the results of VES interpretation by adopting conventional Schlumberger and modified Schlumberger arrays at the same location show some degree of good correlation. This therefore implies that the modified Schlumberger array can be employed instead of the conventional Schlumberger array method in locations where there is space constraint to accommodate the symmetrical electrode configuration of the conventional Schlumberger to reasonable values ofA
B
/
2 which in turn reduces the depth of investigation. In this scenario, the modified Schlumberger array utilises the little space available for its asymmetrical electrode configuration thus allowing study of subsurface to desired depth.



---

*Source: 612303-2015-02-28.xml*