Integrated Geophysical and Geological Mapping of the Lithological and Geological Structures at a Proposed Generator Site

The discovery of oil in commercial quantity in Ghana has led to increase in infrastructural development to provide the needed services to the oil industry. The siting of these industries requires the characterization of the subsurface to determine its suitability to host these facilities. It is in this line that seismic refraction and geo-electrical surveys were conducted on a site near Takoradi, Ghana that has been planned to host a generator and its appurtenances for the oil industry. The objective of the survey was to provide subsurface information critical for the design of earthing system for the proposed generator plant and its appurtenances. The ABEM SAS4000 Terrameter with the Lund imaging system and a 24-channel ABEM Terraloc MK6 were respectively used for the electrical and seismic refraction data collection. The seismic refraction and resistivity results identify three lithological units namely, dry gravel, sand and silty clay layer, water saturated sand and silty clay layer and granodiorites basement rocks. The first layer has a thickness of 10 m and the second layer with thickness in the range of 5 to 10 m. The lithological classi-fications as obtained from these results correlate with the drill logs. Electrical grounding will be suitable at a depth of 10 m beneath resistivity profile two at a distance between 32 to 52 m.


Introduction
The discovery of significant crude oil in Ghana in 2007 and the subsequent pro-duction in 2010, has brought a tremendous increase in infrastructure development especially in and around the Sekondi-Takoradi Municipality. These facilities are used as servicing centres for the oil industry as well as homes for the people working in the oil industry. In order to maintain the integrity of these structures, there is the need for detailed geological knowledge of the subsurface for geotechnical site characterization purposes, especially for the siting of a generator plant, which is the object of study. This knowledge can be acquired by the application of geophysical methods. The use of geophysical techniques continue to increase from the conventional usage of exploring for natural resources into civil engineering, rock mechanics, disaster prevention, environmental preservation, archaeology, impact cratering science among others because of the geological model that is generated from the geophysical data (Takahashi et al., 2006;Knödel et al., 2007).
Integrated geophysical and geological methods have played significant roles in exploring for natural resources, environmental issues and geotechnical investigations.
Electrical resistivity distribution within the subsurface is mapped by injecting current into the ground and measuring the resulting potential differences. Conventionally, this is done either by profiling where the lateral variation of electrical resistivity is measured or vertical electrical sounding (VES) which investigates vertical variation of resistivity. Recently, electrical resistivity tomography (ERT) technique which combines profiling and VES has become the preferred choice for near surface investigations.
For the ERT, a single injection of electrical current can result in multiple potential difference measurements thereby reducing drastically the time spent on the field. ERT has been widely used in engineering, environmental and shallow surface investigations (van Schoor, 2002;Aning et al., 2014;Andrews et al., 2013;Ugwu & Ezema, 2013;Obasi et al., 2015;Binley et al., 2015). The electrical resistivity method has also been used in mining and abandoned mining areas (Martinez-Pagan et al., 2013;Martinez-Lopez et al., 2013) and limestone or dolomite dominant zones which are susceptible to cavity formation due to natural processes (van Schoor, 2002;Metwaly & AlFouzan, 2013) to study the subsurface before projects were embarked on. Subsurface cavities have also been surveyed to assess the likelihood of landslides in an area with the ERT (Panek et al., 2010). Additionally, the ERT technique has found usage in road construction (Farooq et al., 2012), site characterization for borehole development (Ogungbe et al., 2010) and the detection of groundwater pollution (Uchegbulam & Ayolabi, 2014;Wemegah et al., 2017). Other studies also confirm its ability to map shallow anomalies (Kumar, 2012;Nero et al., 2016) and geological units (faults, fractures and quartzite veins) related to the formation of groundwater systems (Mendoza & Dahlin, 2008). In Ghana, Aning et al. (2013aAning et al. ( , 2013b  give a low enough contact resistance for safe grounding of equipment in power plants or substations. For soil corrosion potential, Rajan & Venugopalan (1977), classified soils according to resistivity values as in Table 1.
Low resistivity soils have very high corrosion effect and the same time provide effective grounding. Seismic refraction uses seismic signal that returns to the surface after refraction at the boundaries of subsurface layers with different seismic velocities. Seismic refraction uses this change of direction to obtain some subsurface information. Seismic wave velocity is the most fundamental parameter in seismic methods and it relates to the elastic properties and in addition, mineral content, porosity, lithology, saturation of pore fluid and the degree of compaction. Geophysicists and geotechnical engineers use the seismic refraction method to obtain information about the subsurface geology, particularly the depth to bedrock, rippability and unconsolidated overburden (Milsom & Eriksen, 2013). There are several types of seismic refraction methods depending on the objectives, but the commonest one is the use of first arrivals of p-waves. The refraction method is appropriate in situations where p-wave velocity increases with depth (Knödel et al., 2007;Gadallah & Fisher, 2008). Seismic refraction also provides valuable velocity information for the interpretation of reflection surveys such as saturated aquifer thickness and weathered fault zones. Seismic refraction surveys are useful in delineating the velocity distribution within the subsurface. In 2016, Addai et al. (2016) carried out a seismic refraction survey at KNUST, Kumasi, Ghana and reported p-wave velocity of 4.3 km/s for granite.
Interpretation of geophysical results is enhanced with borehole data and it helps in understanding the conditions of the site because boreholes produce direct subsurface structure information below the ground surface. Borehole measurements give information about the geology of the area for example: lithology, porosity and permeability, presence and orientation of fractures, determination of dip and strike of strata, water quality among others, thus, providing detailed information about the subsurface physical properties. The results from neighbouring boreholes can be correlated to generate cross-section of the site (Knödel et al., 2007;Milsom & Eriksen, 2013;Reynolds, 2011 The objective of this study was to determine the seismic and resistivity parameters, and consequently determine the depth to bedrock and other major structural features critical for siting of the generator plant. The study area is basically a virgin area as far as subsurface geophysical investigations are concerned. The results will therefore serve as a baseline for further studies in other areas of similar geology.

Project Site Description
The research area is located in the Sekondi-Takoradi metropolis of Ghana, which is situated between latitudes 40˚52'30''N and 50˚04'00''N and longitudes 1037˚00''W and 1052˚30''W. The area is within the southwestern part of equatorial zone of Ghana and is 242 km west of Accra, the capital town. The metropolis is bounded to the north by the Mpohor-Wassa District, to the south by the Gulf of Guinea, to the west by the Ahanta West District, and to the east by Shama District (Figure 1). It has fairly uniform temperature, ranging between 220˚C in August and 300˚C in March. The metropolis has a bimodal raining season. The major season runs from April to June and the minor season from September to the end of November with a mean annual rainfall of 2350 mm. Relative humidity is generally high throughout the year between 50% and 70% in the dry season and 75% and 85% in the wet season. The topography is undulating with the central parts of Takoradi lying at an altitude of six meters below the mean sea level. These lowlying areas are interspersed with ridges and hills that range from 30 to 60 m above sea level. Vegetation is highly woodland in the northern and central parts, whilst thicket is intermingled with tall grass species along the coast, especially in areas where there are no permanent crops (Fei-Baffoe et al., 2014).
Geologically, the area falls within the sedimentary basin granitoids which form part of the extensive Eburnean Granitoid plutons that intruded the Early Proterozoic Birimian Supergroup in Ghana (Hirdes et al., 1988Taylor et al., 1992). These granitoids suites in Ghana are characterized into four main groups namely the sedimentary basin granitoid; belt granitoids; and the late post Eburnean, K-rich granitoid and the Winneba granitoids (Hirdes et al., 1988;Leube et al., 1990). The sedimentary basin granitoids within which the study area falls, is a large granitoid batholiths emplaced within the Birimian metasediments. This granitoid coincides with the central axes of the sedimentary basins (Leube et al., 1990). It has wide range of rocks that include quartz, diorite, tonalite to trondhjemite compositional types as well as adamellite and granites. These rocks are characterized by weak foliations. The main ferromagnesian mineral is biotite, which is commonly accompanied by muscovite (Leube et al., 1990). Geochronological analysis done using U/Pb zircon and monazite dating by Hirdes et al. (1992) on this granitoids plutons gives it an age of 2116 ± 2 Ma.

Data Acquisition
Field measurements were conducted using the 24 channel ABEM Terraloc Mk6  and Abem LUND imaging system to take the seismic refraction and apparent resistivity measurements respectively. The electrical resistivity tomography (ERT) technique was employed for the resistivity data collection. Two (2) profiles were surveyed using the gradient array with a minimum electrode separation of 2 m. The Grad4L8 and Grad4S8 protocols of this type of array were made use of. The Grad4L8 is the long layout protocol whilst the Grad4S8 is the short layout protocol and is designed to supplement the former for the shortest electrode spacings only. In all, two 120 m long orthogonal traverse lines R 1 -R 1 and R 2 -R 2 were surveyed.
The seismic refraction survey was conducted using 5 m geophone spacing.
Investigations were carried out on two profiles namely S 1 -S 1 (profile 1) and S 2 -S 2 (profile 2). The profile S 1 -S 1 was 155 m and S 2 -S 2 was 240 m in length. Shot points along profile 1 were located at 0, 20, 47.5, 77.5, 135 and 155 m. Shot Nine boreholes were drilled and logged across the site with depth ranging between the minimum of 9.7 m for hole 1 to a maximum of 22.2 m for hole 6.

Data Processing
The resistivity data was processed using the Res2DInv software. Bad datum points were first eliminated from the data, after which it was inverted using the Gauss-Newton least squares inversion based on the robust inversion routine. 2D models were obtained in the process and displayed using the 'user defined logarithmic contour intervals'. This step results in even contour values and spacing on the resistivity models for all the profiles for easy comparison. Topographical correction was not done because the surface was scooped and levelled.
The processing of the seismic data was done using ReflexW software version 7.2. All picked travel times were put together and the wave front inversion was used for the layer reconstruction. The travel time curves of the single shots for each profile were put together and assigned layer velocities. No static correction was done because the seismic lines were on a flat surface and limited filtering was done. The p-wave velocities of the acoustic layers and depths were done using the intercept time method.
The drill logs were generated with Goldensoftware Strater. The depth and lithology logs were created one after the other using the interval lithology information that was obtained during the logging of the hole. Figure 2 shows electrical resistivity tomographies for profiles 1 and 2.  m distance points could be a saturated silty clay or burried river channel. Figure 3(a) and Figure 3(b) are the seismic refraction models of lines S 1 -S 1 and S 2 -S 2 respectively. It exhibits a three-layer model with undulating interfaces. The topmost layer has a thickness averaging about 10 m and p-wave velocity of < 2 kms −1 . This layer is made up of the dry silty clay, sandy clay and clayey gravels.

Seismic Refraction
The second layer has a velocity between 2.00 -4.70 kms −1 and average thickness of about 10 and 5 m for profile 1 (Figure 3(a)) and profile 2 (Figure 3(b)) respectively. The second layer of profile 1 (Figure 3(a)) thickens at the centre and shrinks towards the end. The layer consists of mainly of saturated sandy and silty clays. Beneath it is the last layer with velocity > 4.70 kms −1 . This layer comprises the basement granodiorites.

Drilling
The

Discussion
Electrical resistivity and seismic velocity depend on whether the subsurface materials are massive or disseminated and wet or dry. The electrical resistivity and seismic velocity within the same formation can differ depending on these factors. The electrical resistivity results compare very well with that of the refraction seismic for the near surface in that they all have a top layer with thickness of about 10 m. This topmost layer which has a moisture content less than the one below has a higher resistivity and a lower seismic velocity. Due to evapotranspiration and seepage, resistivity values near the surface are higher than those deep down within the same formation because the near surface is drier and for seismic waves, the velocity tends to increase with increasing moisture content. The seismic results correlate very well with the drill logs. Borehole 5 which is at a point distance of 65 m (Figure 1) on profile 2, the basement granodiorites are at a depth of about 15 m on the seismic model (Figure 3(b)) and at 12 m in borehole 5 (Figure 4). Borehole 6 which is about 2 m from profile 2 (Figure 1) and at a point distance of 130 m, the granodiorite bedrock is at a depth of about 16 m on both the seismic profile 2 (Figure 3(b)) and the drill log ( Figure 5). The summary of the various seismic ground properties obtained from the seismic survey are shown in Table 2.  Clayey gravels were found in the drill holes which are on the eastern side Corrosion will be severe below 10 m on profile 2 (Figure 2(b)) between 32 and 52 m distance points as a result of the low resistivity. Corrosion will be very mild on the entire profile 1 (Figure 2(a)) because the resistivity beneath this survey line is greater than 50 Ωm. This is based on the resistivity-corrosivity classification of soils as reported by Rajan & Venugopalan (1977). Grounding for the proposed generator plant will be ideal at a depth of 10 m between 32 and 52 m distance points on profile 2 (Figure 2(b)) due to the low resistivity of this portion of the subsurface (Schaefer, 1955).

Conclusion
In this study, we conducted electrical resistivity, seismic refraction, as well as drilling in a site earmarked for sitting a generator plant. The results from the Res2DInv resistivity and ReflexW seismic models correlated well with the drill logs. From the models, we were able to delineate low and high resistivity zones which are good sites for grounding and positioning the generator and its accessories respectively. We were able to map the basement granodiorites. It's evident that a single geophysical method is not enough to identify all possible subsurface features at a site and the integrated approach must be adopted. The analysis of