The Use of Seismic Refraction and Self
Potential Surveys by |
ABSTRACTAlthough 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. |
INTRODUCTIONIn 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 Bureaus 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 SURVEYSThe 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
|
Figure 2
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
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 1The 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 2This 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
|
Figure 5a & bDam / Axis / Comon Offset / Channel 3 |
Figure 6
Dam Axis / Difference Traces |
EXAMPLE 3This 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 SURVEYSSelf 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:
Potentials resulting from lithological and structural factors can generally be classified as follows:
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. |
|
|
| 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 10
|
Figure 11
|
TAILINGS DAM NO. 7
|
Figure 14
|
Figure 15
|
CONCLUSIONSThe 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. |