Cross Borehole Seismic Tomography Applications To Mineral Development

Introduction Cross borehole seismic tomography (CBST) is a geophysical technique that allows definition of the distribution of some physical rock property between two boreholes. The method is similar to the computer-assisted tomography (CAT) scans used by the medical industry. This paper describes the method, its applications and the developing trends in its use that interest the mining industry. The idea of using tomography in connection with geophysical methods for measuring rock properties between boreholes has been around for three or four decades. Until the advent of high-speed, large-capacity computers, the necessary calculations were too formidable for production use. Another factor in the development of the technique has been the refinement of a mathematical base to meet the needs of the medical industry. The geophysical applications have benefited from this refinement. An extensive discussion of the basic theory is given in Kak and Slaney (1988). CBST, although widely researched and written about, is still not extensively put to use in mining and engineering projects. Most real world examples in the literature have come from research carried out by government laboratories or by large companies. An exception is the regular use by Prospection Geophysique Francaise (Cottin et al., 1986) in damsite investigations. The Japanese appear to be doing extensive research and testing of CBST's uses in civil engineering projects (Saito et al., 1988). Canada is researching the use of the method in connection with nuclear waste site investigations (Wong et al., 1987). Finland has tested the method in connection with a foundation investigation for a nuclear power plant (Cosma et al., 1984). In the US, the US Bureau of Mines has taken an interest in using CBST as a means to enhance in situ mining operations (Tweeton, 1988). Major oil companies are examining the method as a tool in oil field development (Justice et al., 1989). Description of the method CBST involves placing a series of receivers down one borehole and recording transmissions from an energy source placed successively at a series of positions in another hole. The objective is to generate a dense network of seismic raypaths between the holes. The denser the network of raypaths, the better the features between the boreholes are resolved. Either seismic travel times or sine wave amplitudes are measured by the receivers. Analysis of travel times between boreholes develops a picture of rock velocities between the holes. Analysis of seismic wave amplitudes develops a picture of the distribution of energy transmission characteristics between the holes. These seismic properties can be used to define occurrences of fracture zones, lenticular sequences or other important geologic data not apparent from examination of the borehole logs. When longitudinal (compressional) and transverse (shear) waves are analyzed, the physical properties, such as elastic moduli, can also be mapped. When wave amplitudes are measured, rather extensive corrections to the data are necessary. There is some disagreement among researchers as to the feasibility of making some of the corrections. For this reason, this paper will concentrate primarily on travel time analysis. Wave amplitude analysis will be mentioned occasionally only in passing. This does not mean that amplitude analysis does not provide important data on rock parameters, only that for the present there are limitations on its uses. Vertical seismic profiling (VSP), extensively used in recent years in oil exploration, is similar to CBST. Instead of using two boreholes, however, the VSP method uses an array of seismic sources or detectors down one hole and an array of detectors or sources along the surface radiating outward from the hole. While adequate for determining vertical velocity distributions, the VSP method is not, by itself, suitable for most engineering or mining projects. Computer programs used to analyze the field data make use of the fact that the total travel time or change in signal amplitude between boreholes is described by Dines and Lytle (1979) as: [9k = f Rk f( x,Y) ds (1)] where gk = the total time or attenuation along the kth ray path,