The Use of Seismic Refraction and Self Potential Surveys
To Evaluate Existing Embankments

by
G.C. Davenport, L.M. Hadley and J.A. Randall

ABSTRACT

Although long recognized that earth embankments need periodic inspection and monitoring, it has only been recently recognized that geophysical surveys can supplement the results of standard inspection and monitoring techniques. Geophysical surveys have been performed on a number of Bureau of Reclamation dams. Many of these structures were built over 50 years ago, and some of them are in areas of high seismic risk. Common offset seismic refraction surveys have been performed to provide a cross-section of the dam embankment and foundation materials in terms of seismic characteristics. Self Potential (SP) surveys have been performed to investigate seepage conditions within dam embankment and foundation materials, and abutment materials.

The purpose of this paper is to demonstrate that these geophysical techniques are also applicable to tailing dams. A seismic refraction survey and self potential survey were conducted on a large tailings dam in the State of Guanajuato, Mexico as a demonstration. Both surveys indicate that the embankment materials are not uniform, and seepage conditions exist.

These geophysical surveys may find future use as monitoring tools applied to earthern dams.

 

INTRODUCTION

In 1978, approximately 5 billion dollars worth of crops were harvested from farm receiving water from Bureau of Reclamation Dams. Bureau operated dam also supplied over 500 billion gallons of water for municipal and industrial use and 11 million kilowatts of electrical power in that year. In the event of a dam failure, the economic loss as well as the potential hazard to life and property could be enormous.

Many of the Bureau’s dams are earth embankments built over 50 years ago. Some of these structures are located in areas of high seismic risk. Although long recognized that earth embankments need periodic inspection and monitoring, it has only been recently recognized that geophysical surveys can supplement the results of standard inspection and monitoring techniques.

Bureau of Reclamation personnel modified a seismic refraction technique, and have pioneered the use of that technique (the common offset seismic method) to evaluate the condition of earth embankment dams. This technique is currently being used as part of the Safety Evaluation of Existing Dams (SEED) Program in the western United States. Approximately 25 embankment evaluations utilizing this refraction technique have been made to date.

In addition to the seismic surveys, self potential (SP) surveys have been run on a number of dams to investigate seepage conditions within the embankment, foundation and abutment materials. The results of the self potential surveys, when correlated with the results of a seismic survey for the same dam, provide information on changes in porosity, saturation and permeability.

SEISMIC SURVEYS

The common offset seismic method emphasizes changes in lateral- seismic velocity and/or layer thickness. Seismic techniques were selected over electrical resistivity and other geophysical methods since seismic velocity is generally the most sensitive to slight changes in density and saturation in the types of materials commonly used in earth embankments.

The common offset seismic method is an adaptation of the seismic refraction method. A standard seismic refraction line is laid out using 24 geophones at 10-foot intervals. Energy is input into the ground at a point located some distance from the first geophone, resulting in a single-ended refraction line. A plot of time versus distance (in seconds) will yield the apparent compressional wave velocities of each layer detected by the refraction line. Figure 1 presents a schematic diagram of a single-ended refraction line and the resulting time-distance plot. After a single-ended refraction line is completed, the shotpoint (point of energy input) and geophones are advanced across the dam crest or toe at 50-foot intervals, where another single-ended refraction line is completed. This process is continued for the total length of the dam crest or toe. This sequence results in a compilation of seismic energy arrival times at a number of different geophones from a number of different shotpoints. However, the distance from a shotpoint on the first seismic line to, for example, geophone number 9 of that seismic line will be the same distance to geophone number 9 of any seismic line from the shotpoint for that line. Since this represents a common geophone to shotpoint offset (distance), the method of data compilation is termed the common offset method. Figure 2 presents a schematic diagram of this method. The resulting data for a particular geophone (1-24) for all single-ended refraction lines is displayed as a common offset sort, as shown in Figure 3.

For uniform embankment conditions, these arrivals would form a straight line of equal time. Since the field and processing parameters are held constant for each profile, any anomalies observed in the seismic data must be the result of changes in the embankment or foundation materials.

Velocities determined from the single-ended refraction time-distance graphs are compared with variations observed on the common offset sort. For direct arrivals, variations should be very similar since all observed anomalies in the common offset sort must be due to velocity variations. Variations observed on more distant channels (generally refracted arrivals) can be due to changes in velocity, changes in thickness or encountering a different refractor.

The compressional wave velocity is affected by many conditions. However, in sediments (or compacted soil) the primary factors affecting the compressional wave velocity are density, porosity and saturation. In lithified materials (rock) , factors such as cementation, fracturing, alteration and stress generally have a greater affect on the velocity.

The compressional wave velocities determined from the single-ended refraction lines are only apparent velocities in that they are affected by dipping layers. To overcome this shortcoming reversed seismic refraction profiles are performed at selected locations on the dam crest and toe. From these reversed profiles true compressional wave velocities of the embankment and foundation materials can be determined.

Figure 1

wpe1.jpg (15006 bytes)

v1= compressional (p) wave velocity first layer
v2= compressional (p) wave velocity second layer
- Upper part of figure represents a basic layout for a single-ended seismic refraction line.

-Lower part of figure represents a standard time-distance plot resulting from a single-ended seismic refraction line.

Figure 2

wpe4.jpg (20576 bytes)

a. Illustration of "roll along" field technique for gathering seismic refraction data for the common offset sort.

b. Illustration of how the first break data, gethered by shotpoint, will form a line of equal time if subsurface conditions are uniform. This plot also illustrates why low-numbered geophones must be used to evaluate near surface conditions and higher numbered geophones used to evaluate conditions on the refractor.

Figure 3

wpe5.jpg (42674 bytes)

Dam Axis/ Common Offset

EXAMPLES

The illustration shown on Figure 3 was taken from a dam which appears very uniform based on the seismic data. The following examples are presented to illustrate some of the anomalies observed at various dam sites.

EXAMPLE 1

The dam embankment at this site is roughly 200 feet high and 750 feet long. The dam is a zoned earthfill structure which was completed in 1966. The common offset sorts,Figure 4 show a bowed effect resulting from progressively earlier arrivals (increased velocity) towards the center of the dam. The sorts also show that the effect is more pronounced with depth (Channel 13). It is believed that this effect is due to increasing density towards the center of the dam due to differential compaction.

EXAMPLE 2

This rolled earthfill dam is roughly 115 feet high and 900 feet long. Figures 5a and 5b show the common off set sorts for a low- and a high-numbered channel respectively. Figure 5a indicates very uniform conditions in the near surface materials. Figure 5b indicates an abrupt change in arrival times on the right side of the figure. This would indicate a much faster velocity or a sudden decrease in depth. The time-distance analysis for this data indicates that all the arrivals are direct arrivals with the exception of the far end of the profile. It appears, according to construction drawings, that the anomaly observed on the common off set data is due to an abrupt thinning of the embankment materials over the abutment.

Surface wave velocity and character are also analyzed in the common offset method. Traces from various high numbered channels are added or subtracted to create a simple spatial-wavelength filter. The far channels are used because surface waves generally are not well developed within one or two wavelengths from the energy source. This data is presented as another plot, called the difference trace plot. Plots representing wavelengths from 10 to 140 feet are generated. These plots are compared with the common offset plots and visually inspected for surface wave alignment. The wavelength with the best alignment is considered to be an indication of the layer thickness in which that surface wave is generated (Figure 6). This is based on a simple rule-of-thumb relationship which applies only under certain conditions. Although these conditions may not always be met, the difference trace plot gives an indication of anomalous conditions and yields some information on the deeper materials. This presentation is helpful, as the refraction data rarely penetrates deeper than about 70 feet, due to the limited length of the refraction lines. The difference trace plot is not entirely understood in term of knowing what physical conditions cause variations in the surface wave. However, when the surface wave is well developed, it gives an indication of changing layer thicknesses, and seems to enhance anomalous areas observed on the common offset sorts.

Figure 4

wpe6.jpg (38153 bytes)

Figure 5a & b

Dam / Axis / Comon Offset / Channel 3wpe7.jpg (60125 bytes)
Dam Axis / Common Offset / Channel 13

Figure 6

wpe8.jpg (38947 bytes)

Dam Axis / Difference Traces

EXAMPLE 3

This dam is approximately 90 feet high and 14,000 feet long. The dam is a homogeneous earthfill structure. The common offset sort, Figure 7a indicates a relatively uniform section on the left third of the dam and a fairly uniform, but different section on the right third of the dam. The central portion shows some rather dramatic changes between the two. Other common offset sorts show the same pattern, with the anomalies becoming more pronounced with depth. The surface wave data, Figure 7b, shows these anomalous areas even more distinctly.

Additionally, under certain conditions, the surface wave velocity is close to the shear wave velocity of a material, and thus the surface wave velocity can be used to obtain an estimate of shear wave velocity for use in a preliminary dynamic analysis of the embankment or foundation materials.

The results of the common offset seismic surveys yield a visual cross-section of the nipper portions of a dam embankment and foundation materials in terms of seismic characteristics. Anomalous areas observed on the common offset sorts and/or the difference trace plots are identified and related to stationing along the crest or toe of a dam. The seismic data is then compared to drawings of the dam (when available) in order to eliminate anomalies associated with the dam structure itself (outlet works, abutment contacts, etc.). The remaining anomalies are then utilized to outline a drilling program, if warranted.

SELF POTENTIAL SURVEYS

Self potential is a natural electrical potential existing within the earth which arises from a number of causes. These causes can be broadly classified into two groups (excluding man-made causes):

(1) Background potentials, which are primarily a result of, (a) two electrolytes of different concentration being in contact with one another, (b) electrolytes flowing through a capillary system or porous media, (c) an electrolyte in contact with a solid, and (d) electromagnetically induced telluric currents.

(2) Mineralization potentials, which are primarily the results of chemical concentration cells formed when conductive mineral deposits, such as graphite or sulfide, are intersected by the water table.

Mineralization potentials are almost always negative, and may have values up to several hundred millivolts. Background potentials can be either positive or negative, and usually have values of only a few tens of millivolts.

The background potentials developed by electrolytes flowing through a capillary system, or porous media, called electro-filtration or streaming potentials, are used for the study of seepage. As water flows through a capillary system, it collects and transports positive ions from the surrounding materials. The positive ions accumulate at the exit point of the capillary, leaving a net positive charge. The untransported negative ions accumulate at the entry point of the capillary, thus leaving a net negative charge. If the streaming potentials developed by this process are of sufficient magnitude to measure, the entry point and the exit point of zones of concentrated seepage may be determined due to the negative and positive (respectively) self potential anomalies.

Figure 8 indicates the type of self potential anomaly, in millivolts, that can be produced by an artesian water system, while Figure 9 presents the type of self potential anomaly that can be produced by ground water draining or flowing into a fractured bedrock system. Self potential surveys have been successfully used to map reservoir leakage and water seepage paths through embankments and foundation materials.

Self potentials are measured by the use of nonpolarizable electrodes embedded in soil or water. These types of electrodes are used to prevent electrode potentials from developing and masking the self potentials. The measured value of the self potential is recorded on a sensitive, high impedance multimeter (or d.c. voltmeter). The field procedure consists of placing one electrode (base electrode) in a quiet area, such as the reservoir. Another electrode is moved at selected intervals along each of the profiles, and the potential between the base electrode and the moving electrode is measured, recorded and plotted at each location. Where changes in the measured self potential of over 10 millivolts are observed along a profile, detailed measurements are made at small intervals to better define the limits of these anomalous zones. Figure 11 presents the results of a self potential profile made over a known leak in the downstream toe of a dam.

The results of a self potential survey are usually presented as a set of profiles, plotting self potential values against the distance of the measuring electrode from the base electrode. The interpretation of this data is mainly qualitative. The magnitude and sign (±) of self potentials will be affected by seepage flow (streaming potentials) and both lithological and structural factors.

Potentials resulting from seepage flow can generally be classified as follows:

  1. Positive anomalies can represent subsurface water flow, and areas of water discharge.

  2. Negative anomalies can represent areas of water infiltration.

Potentials resulting from lithological and structural factors can generally be classified as follows:

  1. Positive anomalies can represent areas of a higher content of clayey material.

  2. Negative anomalies can represent zones of accumulations of coarse detritus material.

  3. A sudden change in the magnitude and/or sign of the anomaly can represent a geologic contact.

  4. Negative anomalies in a homogeneous rock formation can represent areas of fracturing.

It is evident that to separate streaming potentials from lithologic and structural potentials, care must be exercised in the field to note all soil and geological changes along each profile. In many cases, this information may be sufficient to eliminate certain factors from causing a self potential anomaly. For example, to determine if a negative anomaly is the result of water infiltration or fractured rock may be very difficult; however the geological position of the anomaly, its size and shape, observations of local outcrops, and correlation of the anomaly with known subsurface geological conditions may all be factors that can be used to eliminate water infiltration as the cause of the anomaly.

When both seismic refraction surveys and self potential surveys have been performed in a dam embankment, it may be possible to correlate the results of both surveys. The following table presents a summary of those factors that may cause changes in the compressional wave velocity and/or self potential values of a material. These factors are only those that are common to both types of geophysical surveying. There are other factors, such as density, that are not common to both seismic and self potential surveying; and these have not been evaluated. One major factor not presented in this table, but common to both types of surveying, is a change in the type of material. When considering the dam construction, it has been assumed that the materials are relatively uniform within each material zone.


Table 1
 
Increase in Compressional Wave Velocity Incease in Self Potential
Decrease in porosity
increase in saturation
decrease in permeability

Decrease in Compressional Wave Velocity

increase in porosity
increase in saturation
increase in permeability

Decrease in Self Potential

increase in porosity
decrease in saturation
decrease in permeability
decrease in porosity
decrease in saturation
decrease in permeability

Figure 7a & b

Dam Axis / common offset sort for geophone 8

wpe10.jpg (42596 bytes)

Dam Axis / difference trace 20-22 / wavelength = 40 ft


Figure 8

wpe11.jpg (6856 bytes)


Figure 9

wpe13.jpg (6769 bytes)

Figure 10

wpe14.jpg (11549 bytes)

Figure 11

wpe15.jpg (13303 bytes)

 

TAILINGS DAM NO. 7
GUANAJUATO, MEXICO

Construction was started on the No. 7 tailings dam in early 1982. The upper embankment, which consists of compacted conglomerate, is currently 35 meters high and has a crest length of 200 meters. Embankment material is placed and compacted in 5 - 7 meter lifts. Each succeedingly higher embankment is constructed upstream of the embankment below. The tailings are impounded bo a height of 10 meters below the dam crest. Figure 12 presents a cross section of this tailings dam. The present tailings pond is approximately 150 meters wide and 600 meters long. The dam straddles the Veta Madre, a fault with approximately 700 meters of vertical displacement. The right abutment, which is in the hanging wall of the fault, consists of upper Eocene conglomerate, while the left abutment consists of Jurassic phyllite.

As part of an annual field course in geophysical methods, a seismic refraction profile was performed along the crest of the dam, and self potential profiles were run along the crest of the dam, the upstream dam face and the downstream berm.

The only seismic equipment available for this survey was a single channel (one geophone) seismograph, therefore the common offset method was not used for the seismic work. The three seismic refraction profiles run on the crest of the dam were, however, all reversed to obtain true compressional wave velocities of the embankment materials. The maximum depth of penetration of the seismic survey was 25 meters. Reliable surface wave information was not obtained during the survey, thus penetration of the embankment deeper than 25 meters was not possible.

The resulting seismic section of the embankment is shown in the upper part of Figure 14. It is readily apparent that the embankment does not consist of homogeneous materials. The upper 10 meters of the embankment consist of materials with compressional velocities ranging from 200 to 320 meters per second (660 to 1,050 feet per second). This range of velocities is representative of loose, dry, angular materials. Between stations 0+50 and 0+70, the near surface materials have a compressional wave velocity of 750 meters per second (2,460 feet per second). This material extends downward, underlying, as a second layer, the right portion of the embankment from stations 0+70 bo 1+68 (the end of the seismic surveys). The compressional wave velocity of this material (750 meters per second) is representative of an unsaturated, relatively well compacted material.

Between stations 0+00 and 0+50, the upper layer is underlain by material having an expressional wave velocity of 280 meters per second (920 feet per second). This velocity is representative of loose, dry, granular material. This layer is underlain by material with a compressional wave velocity of 900 meters per second (2,950 feet per second) which may represent unsaturated, relatively well compacted material or badly fractured phyllite in the left abutment.

The results of the self potential survey conducted along the dam crest are shown in the lower part of Figure 14.

Based on the self potential survey data, the area between stations 0+10 and 0+60 appears to be an area of seepage as does the area between stations 0+75 and 1+15. The sharp spike at station 0+65, with a value of approximately +50 millivolts, is believed due to the boundary effects between the 280 and 750 meters per second materials. The sharp anomalies centered on stations 1+25 and 1+35 appear to be narrow areas of active seepage. The smaller anomaly centered on station 1+60 may be due to a small seepage area or a minor change in the embankment materials.

The self potential survey conducted along the upstream toe, in the tailings, shows a general agreement with the dam crest profile (Figure 15). The large amplitude, broad anomaly between stations 0+20 and 0+80 is believed due to the boundary effects of the different embankment materials. A general increase in self potential values is noted from station 1+20 to station 2+00. This agrees with the increase noted on the dam crest profile. The two sharp anomalies (stations 1+25 and 1+35) noted on the dam crest profile are not evident. This suggests that these two seepage areas are less than 20 meters wide (the profile was run with 20 meter electrode spacings) on the upstream toe.

Figure 15 also shows the results of the self potential profile conducted on the downstream berm (also using 20 meter electrode spacings). Only two areas are of note on this profile. The small (8 millivolt amplitude) anomaly at station 0+70 may represent a minor seepage path of water out of the embankment (hence the positive sign) , or the anomaly may be due to the boundary effects of the different embankment materials. (In the latter case, the small size of this anomaly may be due to diminishing effects of the potential due to the thickness of the embankment along the downstream berm.) The other area of note is between stations 1+40 and 1+60. There is a sharp rise in the values of self potential to 500 millivolts at station 1+60. A six inch diameter water pipe, approximately three inches deep was discovered to be broken and leaking water at station 1+60.

Figure 14

wpe16.jpg (22777 bytes)

Figure 15

wpe17.jpg (30956 bytes)

 

CONCLUSIONS

The common offset and self potential methods, as described in this paper, are used in order to accumulate data very rapidly. Currently these methods have been used to investigate a number of large dams. A number of anomalous areas have been observed, including those outlined in the examples above. Drilling is presently underway at one dam to investigate the sources of the observed anomalies. It is anticipated that once the drilling information is available, a correlation between seismic characteristics and physical properties can be made. If such a correlation can be demonstrated, the amount of drilling required in future investigations could be significantly reduced. Even without this correlation, the method is a quick and economical means of characterizing embankment and foundation materials.

The self potential surveys offer an inexpensive method of seepage detection, and, when correlated with seismic surveys, self potential surveys can provide added information on porosity, saturation and/or permeability changes in embankment materials.

These two methods may also find use in the future as monitoring tools by establishing a baseline data file immediately after dam construction. Subsequent surveys can be run after initial filling and after an event such as rapid drawdown or an earthquake to determine if any significant changes have occurred in the seismic or self potential data.