Application seismic refraction method to characterize the building site of Injibara University building site,injibara Ethiopia

doi:10.21203/rs.3.rs-1194882/v1

Download PDF

Research Article

Application seismic refraction method to characterize the building site of Injibara University building site,injibara Ethiopia

https://doi.org/10.21203/rs.3.rs-1194882/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

Geophysical investigation using seismic refraction method was conducted for engineering characterization of the foundation conditions of Injibara University buildings construction site located in Injibara town of Amhara Regional State, northwestern Ethiopia. The principal objective of the research was studying the suitability of the foundation earth materials underlying the site, where Injibara University is established. The seven refraction seismic spreads, seismic velocity models interpretation have provided valuable geotechnical information incorporated with available geologic information in the study area. Interpretation of geophysical data revealed that the subsurface geology of the area is composed of three layers. The topsoil consisted of clay, silt and sand mixtures having a 1-4 m thickness and 255-510 m/s p-wave velocity ranges are mapped over the whole area. The second layer attributed to the highly weathered and fractured vesicular basalt is characterized by 948-1802 m/s P-wave velocity range and revealed somewhat undulating morphology. The depth extent of this layer varies from about 10m on the North West end and southeastern parts and to about 27m around the central part. The third layer occurred in the depth range of 10-27m is characterized by greater than 2550m/s average high p-wave velocity and it is due to moderately weathered and fractured basaltic bedrock, which is deeper near to the center of the profiles and gets shallower towards North West end and southeastern portions. Besides, analyses of collected data have suggested the possible locations of minor structural discontinuities (maybe local fractures).The geophysical results show that the bedrock is found at shallow depth in the northwestern end and southeastern part of the study area, whereas in the central part of the survey area the bedrock is found relatively at high depth. Therefore, setting the building foundation is more recommended in the southeastern part of the construction site.

Geophysics

Geophysical investigation

construction site

Velocity Model

Seismic refraction

P-wave velocity

Construction of sustainable civil engineering structures require profound knowledge about the characteristics of subsurface earth materials, particularly physical peroprties of the underlying rocks/soils, distribution of tectonic elements, contents of moisture or fluid within them. Discontinuities in the form of bedding planes, joints, faults and folds highly determine the physical strength of rocks. Similarly, properties of materials filling voids, such as pure/mineralized water, air or both in unconsolidated soils or fractured rocks influence their physical characteristics. Therefore, the stability of civil engineering constructions depend on the correct assessment of the various physical and geotechnical properties of the underlying earth materials where the structures are intended to be constructed (Johnson and DeGraff, 1988).

On the other hand, constructions undertaken over formation lacking bearing capacities often result in failures, manifested by cracks, settlements, displacements or total collapses. Particularly, those structures erected over areas where expansive soils are widely distributed demand special attention as their shrinking and swelling characteristics can easily cause damages due to their property variations as a result of moisture / fluid content fluctuations associated with seasonal changes.

Therefore, geotechnical investigation of any construction sites is essential to obtain reliable inputs that enable to develop economically and technically feasible structural designs incorporating mitigation measures to anticipated geo-hazard events.

Like elsewhere in the world, in Ethiopia public officials require geotechnical investigation data acquired in accordance with the Ethiopian Building Construction Code with accompanying recommendations prior to issuing a building permit in order to protect the safety of the public the surrounding environment (Dagnachew, D., 2011).

Unlike drilling, pitting and trenching, geophysical methods are environmentally safe and also do not cause any substantial damages/ concerns to the communities. Geophysical measurement responds to change in the physical, chemical, mechanical, elastic, radioactivity or thermal properties of the underlying earth materials. Because of such diverse characteristics, usually one or more of the properties correspond to certain features of earth materials, i.e., contact, discontinuity (fracture/fault zones).

For engineering applications, seismic refraction method is widely used to map the subsurface structures. This method depends on the acoustic impedance contrast of the subsurface materials.

The unique tectonic setting of Ethiopia results in complex geological and geo-morphological setups where along with these and continuously deteriorating environmental conditions, the country is very vulnerable for such geo-hazard risks, as volcanic, seismic, landslide and alike. Every year Ethiopia allocates quite a substantial amount of budget to the expansion of infrastructures: roads, bridges, dams and building complexes. Particularly, to expand access to education the construction of universities are taking place in different parts of the country and among these is Injibara University.

To study the foundation conditions at site and evaluate its suitability for erecting a four story building to be used as dormitory for student, subsurface investigations were carried out employing seismic refraction method with an ultimate objective of generating inputs for civil applications.

The research area, Injibara University site, is situated, in Awi Zone, Amhara Regional State about 447km NW of Addis Ababa (Figure1.1). It is bounded by UTM coordinate 271602-272203m Easting and 1211020- 1211378m Northing, and characterized by flat to gently sloping topography bounded by mountains and small hills from the western and southern side. It has an average elevation of about 2552m amsl and located only about 1.5 km SW of the city center, just on the Injibara-Chagni asphalt road.

The geology of the study area and its surrounding is dominated by the following rock units, scoraceous basalt, vesicular basalt, scoria, trachyte and recent alluvial deposits.

Vesicular Basalt: This rock unit is covered large part of the study area and it is clearly out cropped at Ayo river bank and Mesni,Sutena and Zerket stream banks and beds. Ridges surrounding Injibara town are also covered by this unit. Weathering and fracturing are affecting this unit and it is not filled by secondary minerals. Previously drilled well data result indicates this unit gives potential amount of water to the wells.

Scoraceous Basalt: Dark to brown color, weathering, fracturing and vesiculation are the major characteristics of this unit. This unit is dominantly observed at Ayo riverbank and Sutena streambed.

Scoria: This unit is found in western part from the study area and it is clearly out cropped along Injibara-Chagni road cut. It is weathered and its thickness reaches up to 15 meter. This unit is important for ground water storage and conductance, however it is not suitable for engineering foundation purpose.

Trachyte: This unit is found in western ridges of Injibara town that cover small area and it is characterized as fractured, weathered and reddish in color.

Alluvial deposits: Alluvial deposits are covering low land plain part of the study area and it consists of clay, silt, sand, gravel, cobbles and boulders that are basaltic origin. The clay unit covers large part of the study area. Gravels, cobbles and boulders are dominant on the beds of streams and its thickness varies from place to place

Geological Structures

The dominant geological structures observed on the area are lineaments, local fault, and medium to large spacing fractures with three major sets trending N-S, SE-NW and NE-SW direction.

Data acquisition and Instrumentation

The objective of refraction seismic survey was to determine the velocity of elastic waves propagation along different paths within the subsurface and indirectly assess the density characteristics. Every wave reaching the geophone produces a momentary impulse on a record of ground vibration called seismogram.

The survey was made employing the 48-channel Dolang Seismograph(JEA247ESAC500), though in this study only 12 and 24 channels were used at geophone spacing of 4-10m. A 12kg sledge hammer was used as source generate elastic waves propagation along seven lines (spread), i.e, having a length of 55-110m.The sledge hammer was connected to the seismograph using the trigger cable to arm the seismograph as the impact is done. At each geophone and shot points location coordinates and elevation data were recorded using Garmin GPS map 62.

The geophones were laid on the ground at regular interval measured by a tape meter. Geophone spacing is determined based on accessibility conditions and the detail of the required information. The field procedure employed was an in-line spread in which the source and the geophones were placed in a straight line. The geophones were connected to the seismograph via seismic cables. The source was activated at three different points for the six spreads and six different points for one spread along the spread. The locations of the shots were 2m, 22m, 42m, 50m, 70m and 90m on spread 1 and 4m, 44m and 82m on spread 2.

The seismic refraction data were processed using the software programs PickWin95 and Plotrefa from the SeisImager software package. Raw field data in ‘seg2’ format were imported into PickWin95 and the first arrivals of the P-waves were chosen. A band pass filter frequency of 56.9Hz low cutoff and 409.6 Hz high cutoff were applied to take out the low and high frequency noise. These cutoff frequencies were chosen considering the range of seismic frequencies and observing the seismogram (waveform). This was performed for each of the shot points along the spreads. Then, the first arrival picked data were imported into Plotrefa, and a plot of time versus distance was generated (Fig.1.6). Plotrefa automatically checks reciprocal times for multiple shot locations. It is best if the root mean square (RMS) error is less than 5%. Most of the data points in the spreads have RMS values below 5%; however, there are some data points with higher RMS errors. Data processing requires much care and experience as noise can be deceptive in picking first breaks. The main noise sources in the project area were movement of trucks, excavator, people and animals in the proximities of the survey lines.

Layers were assigned by identifying crossover points, which occur where the slope (1/v) changes. The crossover point separating vesicular basalt and basalt is minor, but the change in slope between the top soil and bedrock is distinct. After the layer assignment, a time-term inversion model can be run. Velocity is calculated, and from the model depth is inferred. Velocity models for each spreads were generated using different first arrival picks in order to gain an understanding of model sensitivity.

This was done for the seven spreads and for each of the source. See a typical example of figure 1.5 and figure 1.6.

In this work, interpretation has been made based on the results of seismic refraction with the help of a borehole lithological log data. The borehole data helps to understand the vertical geological section of the study area and to correlate these different units with the seismic refraction velocity model sections. The depth of the boreholes used for lithological correlation is 296m; whereas the depths of the seismic velocity models are about maximum depth 20-35m, i.e. the depth of the geophysical sections is smaller as compared to the depth of the borehole depth. A borehole used in the interpretation of the geophysical data is found near the boundary of the study area drilled by Amhara design and supervision works enterprise collaboration with Amhara water, irrigation and energy bureau for water supply purpose to Injibara University for the coming year consumption.

Interpretation of seismic refraction data

Data were processed using the software programs PickWin95 and Plotrefa from the SeisImager software package. These programs allow cross-sectional areas of the subsurface beneath each spread to be plotted, thus modeling the bedrock interface using a time-term inversion method was employed. A three-layer model was employed to represent the basalt bedrock, vesicular basalt and a thin top soil layer. Raw field data were imported into PickWin95 and then before picking the first arrival times, band pass frequency filter using upper and lower cut off frequency of 409.6HZ and 56.95HZ respectively was applied to remove the noise and improve signal to noise ratio. The first arrival times were picked up using the auto-pick option of the system for all records and few adjustments were made manually where it seems necessary. This was performed for each of the shot points along the spreads. Then, the first arrival data were imported into Plotrefa, and a plot of time versus distance (1/v) was generated (Fig. 1.6). Plotrefa automatically checks reciprocal times for multiple shot locations. It is best if the root mean square (RMS) error is less than 5%. Most of the data points in the spreads have RMS values below 5%; however, there are some data points with higher RMS errors.

Layers are assigned by identifying crossover points, which occur where the slope of 1/v changes. The crossover point separating vesicular basalt and basalt is minor, but the change in slope between the top soil and bedrock is distinct. After the layer assignment, a time-term inversion model can be run. Velocity is calculated and the model depth is inferred. Velocity models for each spreads were generated using different first arrival picks in order to gain an understanding of model sensitivity.

The model produced using the above software packages were interpreted according to the area geology and the parameters determined from the model. Because of noise data of the seismic refraction spread three, first arrival picking of p-wave travel times was very difficult. Therefore, spread three was not considered in the processing and interpretation part.

Velocity model for spread-1

The seismic refraction velocity model for spread one is presented as in figure1.7. This is 92m long profile which runs NW-SE direction. The model is generated using time-term inversion method. The velocity model represents seismic velocities between 280m/s and 1758m/s. The top most layer shows low p-wave velocity varies 280-520m/s and is about 1.5-3m thick with slight difference along the spread. The velocity indicates that the top layer is composed of soil deposits. The second layer is located at a depth of about2-3m.The velocity of this layer ranges 1400-1355m/s and from the lithological log in the study area this layer is probably weathered and fractured vesicular basalt. The p-wave velocity in the third layer is relatively high and from local geology and lithological log this layer is possibly moderately weathered and fractured basalt. This moderately weathered and fractured basalt is regarded as the bed rock in the building site. The velocity of the model shows that the third layer is relatively strong rock type. Its velocity is 1758m/s. The depth of this layer is varies along the spread. It extends relatively high depth to the left ends of the model (figure1.7).

Velocity model for Spread-2

The seismic refraction of spread two lie in the same line as that of spread one laying in NW-SE direction and its velocity model is presented in figure1.8 with total spread length of 88m.From the velocity model generated from this profile, the time-term inversion model show three p-wave velocity layers. The model presents seismic velocities between 458m/s and 1500m/s. The top layer shows low p-wave velocity of 458m/s and the thickness of this layer is almost similar along the spread. The p-wave velocity indicates that the top layer consists soil deposits of clay, silt and sand. The p-wave velocity of the second layer is varying 900m/s-1250m/s and from the lithological log in the study area this layer more likely to be slightly weathered and highly fractured vesicular basalt. The p-wave velocity in the third layer is relatively high 1500m/s and from the lithological log this layer is probably moderately weathered and fractured basalt. This moderately weathered and fractured basalt is regarded as the bed rock of the site. The velocity of the model shows that the third layer is relatively competent rock formation. The depth of this layer is different at the left and right ends of the velocity model. It is relatively deeper in the left side of the velocity model than towards the right end. From the velocity model we can see that the velocity change among each layer is not gradational i.e. sharp or abrupt.

Velocity model for spread-4

In spread four, the first Layer in the velocity model essentially shows (fig.1.9) low-velocity material overlying on medium velocity layer. This 300m/s indicates that the top layer is dominated by soil deposits of clay, silt and sand with thickness ranging from 2-3m.

The velocity of the second Layer with p-wave velocity 1200m/s indicates moderately zone which is possibly weathered and highly fractured vesicular basalt from the lithological log near by the boundary of the study area. The depth for this layer extends up to about 12m.

The third layer which has a relatively high p-wave velocity material 3000m/s at the base of the velocity model is interpreted as to represent moderately weathered and fractured basalt and it is considered as bed rock in the area. This layer is suggested relatively good for setting civil structures.

Velocity model for spread-5

The seismic refraction of spread five laying in E-W direction which crosses profile four of the electrical resistivity sounding survey. As shown in velocity model (fig.1.10) the thickness of the first top layer is vary 2-3m with seismic velocity of 255m/s. Information from the lithologic log and the velocity value this layer is more likely made up of top soils of clay, silt, and sand. The second layer with p-wave velocity value 1092m/s is corresponding to weathered and highly fractured vesicular basalt. This layer from the geoelectric section has low resistivity due to its high moisture content. The third layer has relatively high p-wave velocity 2225m/s and it is found at a depth up to about 10-12m. From the calculated velocity and borehole information, this layer may be moderately weathered and fractured basaltic formation. This layer is relatively competent and therefore it is regarded as the bed rock in the study area.

Spread-6 Velocity model

Spread six is parallel to profile four of the resistivity sounding survey which lies along SW-NE direction. The seismic velocity model of this spread is shown in figure1.11.The top Layer of the velocity model for spread six is about 3-4 m thick with average P-wave velocity of 510m/s. Layer 2 and 3 have average velocities of 948m/s and 2555m/s respectively. The third Layer is buried about 11-13 m deep in the spread line. The second layer of spread six has relatively low velocity as compared to other spreads in the survey area. This may due to the presence of completely weathered and fractured vesicular basalt. Some intercalation of clay material in the vesicular basalt.

Velocity model for spread-7

Velocity model generated for spread seven located to the northwestern end of the survey area as shown in figure 4.5. As the p-wave velocity, model figure 1.12 shows the thickness of the top most layer varies from 2-3m with velocity of 273m/s. This very low p-wave velocity indicates the top layer is composed of clay, silt and sand deposits. From the lithological log in the study area the second layer is more likely to be weathered and highly fractured vesicular basalt with p-wave velocity varies 1400-1802m/s. This layer is mapped at depth ranges 3-15m deep with irregular morphology. Moderately weathered and fractured basalt is mapped below this layer with velocity of 2527m/s. The p-wave velocity of this layer suggests that it is competent enough to be the bed rock of the study area and it is located at a depth of about 20-25m.From the velocity model we can point out that for each layer velocity it is changed abruptly, it does not show gradational change. Relatively this bed rock is one the most competent one. The competent rock formation is relatively at shallow depth in the left side of the velocity model as compared to right end of the spread.

Geophysical studies involving seismic refraction method was carried out for engineering site characterization of the building construction site at Injibara University, southwest of Injibara town, Northwestern Ethiopia. Based on the results of collected data the following conclusions are made:

Based on geo-seismic property contrasts, 255-3000 m/s for p-wave velocity ranges, three different subsurface layers the study area are delineated that are associated with different geological units. The top soil, characterized by 255-510 m/s p-wave velocity ranges, is associated with the response of the upper layer composed of clay, silt and sand. The relative p-wave velocity variations are due to heterogeneous nature of these top soils in terms of their compositions, degree of compactions and moisture contents. Over the area, its thickness varies from about 1m at the SE part to 4m on NW flank of the study area. The second layer attributed to the highly weathered and fractured vesicular basalt is characterized by 948-1802 m/s P-wave velocity range and revealed somewhat undulating morphology. The third layer in the study area is described by relatively high p-wave velocity on average greater than 2250 m/s values, which are interpreted as responses of moderately weathered and fractured basaltic bedrock is assumed to acquire suitable geotechnical characteristics to bear loads from heavy civil engineering structures. The depth to the surface of this competent formation ranges from about 10m in the NW end and SE part to 27m near to the central part of the study area.

We confirm that the manuscript “Application of seismic refraction method to characterize Injibara University campus building site, Injibara, northwestern Ethiopia” has been read and approved by all the named authors, and they have contributed significantly to the paper. This is also to declare that the paper has not been published earlier in full or as part, or been sent to some other journal for consideration for publication. The paper is also free from any plagiarism/self-plagiarism. There is no conflicts of interest associated with this paper. The authors would be fully responsible if the paper is found to violate any copyright law in future.

Johnson, R.B., and DeGraff, J.V., 1998.Principles of engineering geology.
Dagnachew, D., (2011).Investigation on some of the engineering characteristics of soils in Adama town, Ethiopia, A.A.U

Download PDF

Version 1

posted

You are reading this latest preprint version

Application seismic refraction method to characterize the building site of Injibara University building site,injibara Ethiopia

Status:

Version 1

Abstract

Figures

1. Introduction

2. Description Of The Study Area

3. Local Geology Of The Study Area And Its Surroundings

Geological Structures

4. Methods And Material

Data acquisition and Instrumentation

5. Processing And Presentation Of Seismic Refraction Data

6. Result And Discussion

Interpretation of seismic refraction data

Velocity model for spread-1

Velocity model for Spread-2

Velocity model for spread-4

Velocity model for spread-5

Spread-6 Velocity model

Velocity model for spread-7

7. Conclusion

Declarations

References

Status:

Version 1