3D Electrical Resistivity Imaging (ERI) for subsurface evaluation in pre-engineering construction site investigation

3D Electrical Resistivity Imaging


The success of engineering structures which are directly established on the surface of the earth is among other factors dependent on the support offered by the foundation materials that bear the load of the structure (Terzaghi et al., 1996). The ability of foundation rocks to offer necessary bearing support for engineering structures depend on the bearing capacity of the rock(s). Geological and structural factors such as mineralogical composition, rock fabric, rock association, degree of weathering, fluid saturation, rock deformation and occurrence of joints and faults would not only define the bearing strength of the foundation rocks, it could also impact spatial variability in the strength of the foundation material arising from heterogeneity in rock’s physical properties (Fishman, 1979). Spatial variability in rock bearing capacity is known to impact stress on the poorly supported engineering structures with associated failure which could occur as total, partial and differential settlement or even total collapse of the structure (Nobahar and Popescu, 2001, El-Kateb et al., 2003, Cho, 2007).

The prevalent occurrence of failed and collapse buildings and other distressed engineering structures in Nigeria with associated loss of lives and properties has necessitated the need to ensure that engineering structures are properly constructed (Ede, 2010). While much attention has been focused on proper engineering design, structural balance, reinforcement, material quality among others, less effort has been channelled to ascertain the integrity of the foundation material that would support the structure. Hence common occurrence of many stress induced cracks and other foundation related defect features on different public and private engineering structures, some of which are simply disasters waiting to occur (Folagbade, 2001, Osinowo et al., 2011, Matawal, 2012).

Geophysical investigation offers a cheaper and faster means of getting detailed and credible information about the subsurface. Different geophysical techniques possess ability to image the ground in terms of rock association, distribution and structural deformation, which offer credible information regarding the support in term of strength that a foundation rock is likely to offer to engineering structures (Olorunfemi and Mesida, 1987, Palacky, 1987, Olayanju et al., 2017). Electrical Resistivity Imaging (ERI) technique has been very effective in illuminating the subsurface and apt at providing information about the rock physical properties for economical, environmental and engineering purposes (Keller & Frischknecht, 1966). Different electrical data acquisition technique as well as electrode and profile configurations have been described to present different desirable subsurface imaging abilities (Telford et al., 1990, Sharma, 1997). Three dimensional (3D) ERI geophysical investigation approach offers ability to characterise the subsurface as well as determine heterogeneity in measured rock properties along the vertical (z) and the two orthogonal horizontal (x & y) directions. The determination of variation in foundation rock’s properties along the three orthogonal directions affords the ability to evaluate the spatial variation in rock strength as imposed by the heterogeneity of rock properties as they vary from place to place (Badmus et al., 2012).

In this study, 3D ERI technique was applied as a pre-engineering construction investigation tool to study a proposed construction site within the University of Ibadan campus, Ibadan, southwestern Nigeria. The study evaluated the suitability of the site’s foundation rocks to offer necessary foundation support to the proposed engineering structure, which become particularly essential due to the prevalence of some distressed engineering structures (office complexes and students hostels), commonly characterised with numerous stress induced cracks, within the campus. The occurrence of stress induced cracks on these distressed buildings have been shown through investigations to be related to the incompetent foundation rocks which fail to offer required support for the entire or part of building due to variation in rock strength imposed by heterogeneity. The complex rock association of Quartz-schist, Amphibolitic-schist and Quartzite rocks underlain by Biotite and Granite Gneiss often predisposes buildings constructed, especially in the southern parts of the University of Ibadan campus to foundation support challenges as these different rock units present different support capacity to the structure founded on them.

This study infer the spatial variability in rock’s strength generated by heterogeneity in rock’s physical properties due to both lateral and vertical variations in rock physical properties measured through subsurface georesistivity investigation along carefully established parallel equally spaced resistivity profiles lines. Through this study, useful information as regards distribution of rock’s physical properties shown by the variation of the measured electrical resistivity are generated and used to guide the design of the structures’ foundation.

Location and the Geology of the Study Area

The study area is a proposed building construction site within the University of Ibadan Campus, in Ibadan, southwestern Nigeria. The site is located in the southwestern part of the University campus, and situated between latitude 7° 26.001N and 7° 27.71N and longitudes 3° 53.001E and 3° 54.121E in the north western part of Ibadan. The geology of Ibadan and environs has been well described (Jones & Hockey, 1964: Grant, 1971, Burke et al., 1976) to consist mainly of metasediment series with associated metaigneous rocks, classified to belong to the Migmatite – Gneiss Complex (MGC) of the Basement Complex of southwestern Nigeria (Burke & Dewey, 1972).

The Basement Complex of southwestern Nigeria consists predominantly of metasedimentary, minor metaigneous, metavolcanic rocks and older granite series which comprises of Schists, Quartzites, Banded and Granite Gneisses and interbedded amphiboles generated from metamorphism of shales, grewackes and sandstones with interbedded basalts (Oyawoye, 1972, Rahaman, 1976). Evidence of Pan African Orogeny and later tectonic imprints have been reported and they are impressed on rocks as folds, faults, rock foliation and schistocity around Ibadan and environs (Burke et al., 1976, Shakleton, 1976).

Field mapping was carried out around the University of Ibadan campus to identify different rocks types and complex rock associations as well as some structural elements around the study area. The dominant rock types include quartz-schist, amphibolitic-schist, quartzite, granite and banded gneisses which dominantly trend along approximately north–south direction. Minor rocks also identified include dolerite and amphibolites. The granite and banded gneisses are generally observed to be mostly low lying and underlay the schist rocks and quartzite which often occur as ridges. The quartzite rocks present clear evidence of schistocity. Other visible tectonic imprints include joints, faults, folds and foliation of gneisses rocks, evident as alternation of dark and light minerals (Burke & Dewey, 1972). The geology of Ibadan and environ showing the University of Ibadan campus as the study area is presented in Fig. 1.

Fig. 1. The geological map of Ibadan and environs (modified after Amanambu, 2015).

Materials and methods

A carefully designed data acquisition procedure was adopted in this study to generate grid of information which were systematically measured to present the variation of ground electrical resistivity along three orthogonal directions. Ground apparent resistivity measurement was carried out using multi-electrode Supersting R8/IP/SP Geo-electrical resistivity meter which has ability to make simultaneous measurements with the aid of built-in internal transmitter and AGI electrode Switch box. This automatically switches between electrodes aligned along the established profile for rapid measurement. Resistivity measurements were carried out along ten (10) established equally spaced resistivity profiles which were 10 m apart. Eighty four (84) geo-referenced electrode stations were occupied on each profile line at 3 m station interval using dipole-dipole electrode configuration (Fig. 2). Dipole-dipole electrode configuration is able to overcome resistivity measurement challenges arising from electrode polarization as a result of reasonable separation between the current and potential electrodes (Sasaki, 1989, White et al., 2003). The electrode configuration is also capable of high integrity data measurement at depth and thus able to generate high resolution electrical resistivity data. Fig. 2 presents the 3D profile layout and the pictures while occupying different profile lines during surveys are presented in Fig. 3.

Fig. 2. 3D resistivity survey profile layout showing parallel equal-spaced profiles.
Fig. 3. Profiles layout for ERI survey.

Ground resistivity measurement involves injecting current (I) into the ground across two current electrodes which are firmly driven into the ground. One of the electrodes act as a sink (−I) for current injection into the ground while the second acts as a return (+I) into the current injecting source and thus guarantee current flow within the subsurface. The resultant potential difference (V) generated due to the flow of injected current is measured via two potential electrodes which are separated by a distance that corresponds to a factor (n) of the spacing (a) between individual current or potential electrodes. A measure of ground resistance to current flow is computed as the ratio of the generated ground voltage resulting from flow of current within the subsurface. The computed ground resistance is converted to apparent ground resistivity (

) by incorporating the geometric factor as () of the electrode measurement configuration (Eq. (1)) (Telford et al., 1990, Sharma, 1997)

The integrity of the generated apparent ground resistivity data from the geophysical survey was ascertained and the data also evaluated for consistency. Subsequently the data were filtered to clean out spurious and spiky data having no geologic relevance. The cleaned data were thereafter inverted using RES2DINV inversion algorithm developed by Loke & Baker, (1996) to generate 2D resistivity distribution model of the subsurface which consist of noise free (reduced to barest minimum) and depth matched resistivity data (Loke, 2000 & Loke, 2001, Loke et al., 2013). This is achieved through the forward modeling subroutine of RES2DINV which calculates the theoretical apparent ground resistivity values based on a non-linear least squares optimisation (Sasaki, 1989, DeGroot-Hedlin and Constable, 1990, Auken and Christiansen, 2004). A high quality result with characteristic low RMS value (<5%) usually obtained after four (4) to six (6) iterations indicate low degree of misfit between the observed and the theoretical model and thus reflect the integrity of the generated model as a representation of the subsurface across the sampled profiles.

The resulting inverted data from the generated 2D geo-resistivity models were sorted based on corresponding profiles into depth locations using 45° depth point projection from the surface current and potential electrodes’ mid points. The sorted data were geo-referenced with the aid of the individual electrode’s Global Positioning System (GPS) readings of the established electrode points. All the geo-referenced inverted resistivity data from all the ten (10) profiles were combined and gridded using 3D kriging gridding algorithm extension of Geosoft to generate 3D inverse resistivity model of the subsurface (Olea, 1974, Swain, 1976, Webring, 1981). This presents variation in ground resistivity along and across the established profiles stations as well as show variation in ground resistivity with depth. The generate 3D inverse resistivity model was further smoothened to generate the distribution of ground resistivity across the proposed engineering construction site.

Results and Interpretations

The subsurface resistivity investigation of proposed engineering construction site generated 2D resistivity models, sections, resistivity distribution depth maps as well as 3D resistivity model of the subsurface. The 2D resistivity distribution models which were generated by inverting the resistivity data acquired from individual profiles present variation in ground resistivity along the occupied electrode stations within each profile. The models present variations in ground resistivity value from near the surface to a depth of about 30 m.

All the 2D resistivity models (Fig. 4) apart from profile 1 presented three layers, which appear relatively well developed (layers well defined) at the centre and the western end of the W–E running profiles with the near surface (0–1.5 m) constituting the first layer. The near surface resistivity distribution indicate highly heterogeneous unit with resistivity value ranging from low (<10 Ω m) through medium (20–450 Ω m) to high resistivity value (>1500 Ω m). The middle layer presents characteristically low resistivity values (8–157 Ω m), which appear very prominent at the centre and western end of the studied field. The thickness of the low resistivity zone varies from less than 3 m to about 17 m.

Fig. 4. 2D inverted resistivity models of occupied profiles at the study area.

The third (lower) layer presents relatively high resistivity value (>450 Ω m), higher than that of the second (middle) layer. The resistivity of the third layer increases with depth except in profile 6 where low resistivity unit is encountered at depth of about 28 m. Relatively high resistivity value (>1650 Ω m) characterises the eastern end (right end) of all the profiles which indicate relatively high resistive section of the field.

The 2D resistivity sections generated from combined inverted resistivity data obtained from the study area also reflect variation in ground resistivity with depth. Georesistivity sections generated along E–W and N–S directions, parallel and perpendicular to the occupied profile stations respectively are presented in Fig. 5 (a & b). The sections generally indicate higher ground resistivity in the eastern end of the field, while the western end presents relatively lower resistivity distribution. The low resistivity zones as shown by the E–W and N–S (Fig. 5 (a & b) appear to divide the field into two, where the low resistivity distribution zone appears restricted to the west. The prominent low resistivity zone observed on the N–S trending georesisitivity sections (Fig. 5b) which occurs at the centre of the field is more extensive along the north south directions but generally decreases westward.

Fig. 5. (a) E–W and (b) N–S trending 2D georesistivity sections across the study area.

A resistivity distribution mesh which combines both the N–S and E–W trending sections (Fig. 6) also confirmed higher resistivity distribution pattern in the east (>650 Ω m) relative to the west (15–300 Ω m). The resistivity mesh also identified some isolated regions of low resistivity zones (<110 Ω m) in the western part of the studied field.

Fig. 6. 2D georesistivity sections mesh across the study area.

Two-dimensional (2D) depth resistivity maps present the georesistivity distribution across the entire field at a specified depth. Here, all georesistivity data points from all the occupied profiles, corresponding to the same inverted depth point were sorted and extracted to generate each 2D resistivity map of the study area. Fig. 7 presents the stacked 2D georesistivity distribution maps generated for different depths, right from near the surface (<1 m) to about 30 m below the ground surface. The stacked resistivity maps indicate the ground resistivity distribution across the field. The individual 2D resistivity depth maps are presented in Fig. 8(a)–(f)

Fig. 7. Stacked 2D georesistivity map of the field at different depths.
Fig. 8. 2D georesistivity distribution maps of the field at different depths.

The 2D georesistivity maps of the proposed engineering construction site indicates the resistivity distribution pattern across the field at different depths and these enabled the evaluation of the resistivity distribution pattern of the entire field at different depths. Fig. 8a presents the field resistivity distribution at very close to the surface (<1.0 m), while Fig. 8(b)–(g) present resistivity distribution across the field at 5, 10, 20, 25 and 30 m depths, respectively. The 2D resistivity distribution maps corroborates the earlier mentioned higher resistivity distribution pattern in the eastern part of the field relative to the west as presented by the georesistivity sections (Fig. 5, Fig. 6). The maps also identified low resistivity zone towards the western part of the site. The various depth resistivity distribution maps indicate increase is ground resistivity value with depth at the eastern part of the field as the resistivity value increases from 450 Ω m at near the ground surface to greater than 4500 Ω m at 30 m depth. The western part of the field on the other hand displayed slightly varied resistivity distribution with some low resistivity zones which dominate the central and the extreme western end of the field. The resistivity distribution in the field ranges from 8 Ω m to 350 Ω m. Low resistivity distribution zone (8–35 Ω m) characterises the south-western part (left hand side) of the field. This zone appears to increase in prominence with depth up till about 10 m as the spread/distribution of the low resistivity zone increases northward as well as westward (Fig. 8).

At deeper depth beyond 10 m, the low resistivity zones become restricted to the extreme west of the field which presents very low resistivity value (<10 Ω m). The resistivity of the low resistivity zone is generally observed to decrease with depth from 25 Ω m at 15 m to less than 8 Ω m at 30 m (Fig. 8).

The summary of ground resistivity distribution of the proposed engineering construction site along three orthogonal directions (x, y & z) is presented in Fig. 9 as 3D subsurface georesistivity distribution model of the field. The 3D model presents the combined georesistivity distribution of all inverted resistivity data from all the occupied ten (10) geoelectric profiles. The figure indicate the variation in ground resistivity distribution along the E–W and N–S directions as well as show variation in subsurface resistivity with depth from the surface to up 30 m below the ground surface. The figure indicate relatively high ground resistivity (450–4500 Ω m) in the eastern part of the field while the western end displayed lower resistivity values (10–50 Ω m). The eastern part of the field displayed consistently relatively high resistivity distribution right from the surface to a depth of about 30 m with resistivity distribution ranging from 450 to 4500 Ω m. The resistivity distribution in the east of the field increases with depth with the deepest part presenting resistivity distribution in excess of 4300 Ω m, which could be attributed to fresh basement resistivity signature. However, the resistivity distribution in the west end of the field varies from medium to low resistivity (10–600 Ω m). The central part of the field presents significantly low resistivity value (10–60 Ω m) which decreases with depth. The low resistivity zone extends in spread north and westward.

Fig. 9. 3D Georesistivity distribution model of the proposed engineering construction site.

Discussion of results

The carefully planned 3 dimensional survey through sets of equally spaced parallel survey profile lines along which geophysical data were acquired enabled the generation of grid of data which shows distribution and variation of ground resistivity data along the two horizontal (x & y) and a vertical (z) orthogonal directions. The variation along the vertical is made possible by measurements along the electrical survey levels where based on continuous increase in electrode spacing from one level to another, the sampling depth increases with electrode spacing. The distribution of subsurface resistivity value has been presented as 2D resistivity sections and maps as well as 3D resistivity model across the studied field. The presented subsurface georesistivity models show how the ground resistivity value varies across the field. From the presented sections, maps and models, it is obvious that the subsurface, up to the imaged depth of 30 m, is not uniform in terms of resistivity distribution (Fig. 4, Fig. 5, Fig. 6, Fig. 7, Fig. 8, Fig. 9). This is apparent as the eastern part of the field presented significantly higher resistivity values than the western part of the field. Isolated low resistivity zones which characterises part of the central and western end of the studied field also alluded to this variation. Furthermore, variation in the resistivity values with depth which characterises both the resistive and less resistive zones in the eastern and western parts of the field respectively further confirm spatial variability of ground physical parameter along the three orthogonal directions.

Several factors have been reported to impact variation of electrical resistivity of the subsurface. These include variation in rock type, rock deformation, rock fabrics, different degree of weathering and water saturation among others (Ward and Fraser, 1967, Zonge, 1972). These factors which are known to impact spatial variation on ground electrical resistivity are also capable of impacting variation in other rock’s physical properties (Utset and Castellanos, 1999, Carroll and Oliver, 2005). For example variation in ground resistivity distribution which arises due to rock deformation or weathering could also impact variation in the rock density which may in turn affects the bearing strength of the rock unit to efficiently support the load of the building imposed upon it (Chandler, 1972, Tiwari and Marui, 2005).

Geotechnical investigation which has the capacity to directly measure ground’s strength to bear the burden of the proposed load is often expensive and labour intensive, thus few points are often selected under assumption that the subsurface is uniform. This assumption is far from reality as subsurface physical properties could vary along both vertical and horizontal directions within short distance (few meters). Adequate subsurface information is essential for site selection and also important to design foundation that would adequately support the intended load of the engineering structure. The best investigation that would aid good foundation design is such that would not only sample the subsurface but also account for spatial variability in the distribution of physical properties of the ground. The result obtained from carefully planned 3D geophysical investigation of proposed engineering construction site could aid point selection in order to capture different zones as presented by the 3D georesistivity model of the subsurface. This way, both the measured ground resistivity parameter and the vertical as well as lateral spatial variation would be accounted for.


The ability of geophysical investigation to measure both subsurface physical parameters such as ground resistivity, as well as determine the spatial variation of the measured parameter along the three orthogonal directions has been demonstrated. The generated sections, maps and 3D models of the subsurface have provided relevant information which is able to guide point selection for geotechnical investigations, especially where the budget is limited and every part of the field could not be covered at very close grid. Geotechnical sampling along carefully selected points and along different delineated varying georesistivity zones, would provide relevant information that would account for heterogeneity of the subsurface and thus guide design of appropriate foundation that will support every part of the building amidst possible variation in rock strength.

Source: 3D Electrical Resistivity Imaging (ERI) for subsurface evaluation in pre-engineering construction site investigation

Authors: Olawale Olakunle Osinowo, Michael Oluseyi Falufosi

1 Trackback / Pingback

  1. Design of rock support system under rockburst condition - Engineering geology

Leave a Reply

Your email address will not be published.