In attempting to detect and map out underground facilities, whether they be large-scale hardened deeply-buried targets (HDBT's) or small-scale tunnels for clandestine border or perimeter crossing, seismic imaging using reflections from the tunnel interface has been seen as one of the better ways to both detect and delineate tunnels from the surface. The large seismic impedance contrast at the tunnel/rock boundary should provide a strong, distinguishable seismic response, but in practice, such strong indicators are often lacking. One explanation for the lack of a good seismic reflection at such a strong contrast boundary is that the damage caused by the tunneling itself creates a zone of altered seismic properties that significantly changes the nature of this boundary. This report examines existing geomechanical data that define the extent of an excavation damage zone around underground tunnels, and the potential impact on rock properties such as P-wave and S-wave velocities. The data presented from this report are associated with sites used for the development of underground repositories for the disposal of radioactive waste; these sites have been excavated in volcanic tuff (Yucca Mountain) and granite (HRL in Sweden, URL in Canada). Using the data from Yucca Mountain, a numerical simulation effort was undertaken to evaluate the effects of the damage zone on seismic responses. Calculations were performed using the parallelized version of the time-domain finitedifference seismic wave propagation code developed in the Geophysics Department at Sandia National Laboratories. From these numerical simulations, the damage zone does not have a significant effect upon the tunnel response, either for a purely elastic case or an anelastic case. However, what was discovered is that the largest responses are not true reflections, but rather reradiated Stoneley waves generated as the air/earth interface of the tunnel. Because of this, data processed in the usual way may not correctly image the tunnel. This report represents a preliminary step in the development of a methodology to convert numerical predictions of rock properties to an estimation of the extent of rock damage around an underground facility and its corresponding seismic velocity, and the corresponding application to design a testing methodology for tunnel detection.
In many parts of the United States, as well as other regions of the world, competing demands for fresh water or water suitable for desalination are outstripping sustainable supplies. In these areas, new water supplies are necessary to sustain economic development and agricultural uses, as well as support expanding populations, particularly in the Southwestern United States. Increasing the supply of water will more than likely come through desalinization of water reservoirs that are not suitable for present use. Surface-deployed seismic and electromagnetic (EM) methods have the potential for addressing these critical issues within large volumes of an aquifer at a lower cost than drilling and sampling. However, for detailed analysis of the water quality, some sampling utilizing boreholes would be required with geophysical methods being employed to extrapolate these sampled results to non-sampled regions of the aquifer. The research in this report addresses using seismic and EM methods in two complimentary ways to aid in the identification of water reservoirs that are suitable for desalinization. The first method uses the seismic data to constrain the earth structure so that detailed EM modeling can estimate the pore water conductivity, and hence the salinity. The second method utilizes the coupling of seismic and EM waves through the seismo-electric (conversion of seismic energy to electrical energy) and the electro-seismic (conversion of electrical energy to seismic energy) to estimate the salinity of the target aquifer. Analytic 1D solutions to coupled pressure and electric wave propagation demonstrate the types of waves one expects when using a seismic or electric source. A 2D seismo-electric/electro-seismic is developed to demonstrate the coupled seismic and EM system. For finite-difference modeling, the seismic and EM wave propagation algorithms are on different spatial and temporal scales. We present a method to solve multiple, finite-difference physics problems that has application beyond the present use. A limited field experiment was conducted to assess the seismo-electric effect. Due to a variety of problems, the observation of the electric field due to a seismic source is not definitive.
The West Pearl Queen is a depleted oil reservoir that has produced approximately 250,000 bbl of oil since 1984. Production had slowed prior to CO{sub 2} injection, but no previous secondary or tertiary recovery methods had been applied. The initial project involved reservoir characterization and field response to injection of CO{sub 2}; the field experiment consisted of injection, soak, and venting. For fifty days (December 20, 2002, to February 11, 2003) 2090 tons of CO{sub 2} were injected into the Shattuck Sandstone Member of the Queen Formation at the West Pearl Queen site. This technical report highlights the test results of the numerous research participants and technical areas from 2006-2008. This work included determination of lateral extents of the permeability units using outcrop observations, core results, and well logs. Pre- and post-injection 3D seismic data were acquired. To aid in interpreting seismic data, we performed numerical simulations of the effects of CO{sub 2} replacement of brine where the reservoir model was based upon correlation lengths established by the permeability studies. These numerical simulations are not intended to replicate field data, but to provide insight of the effects of CO{sub 2}.
Full-waveform seismic reflection responses of an isolated porous sandstone layer are simulated with three-dimensional (3D) isotropic poroelastic and isotropic elastic finite-difference (FD) numerical algorithms. When the pore-filling fluid is brine water with realistic viscosity, there is about a ∼10% difference in synthetic seismograms observed in an AVO recording geometry. These preliminary results suggest that equivalent elastic medium modeling is adequate for general interpretive purposes, but more refined investigations (such as AVO waveform analysis) should account for poroelastic wave propagation effects.
We present a technique that allows for the simultaneous acquisition and interpretation of both shear-wave and compressive-wave 3-D velocities. The technique requires no special seismic sources or array geometries, and is suited to studies with small source-receiver offsets. The method also effectively deals with unwanted seismic arrivals by using the statistical properties of the data itself to discriminate against spurious picks. We demonstrate the technique with a field experiment at the Facility for Analysis, Calibration, and Testing at Sandia National Laboratories, Albuquerque, New Mexico. The resulting 3-D shear-velocity and compressive-velocity distributions are consistent with surface geologic mapping. The averaged velocities and V{sub p}/V{sub s} ratio in the upper 30 meters are also consistent with examples found in the scientific literature.
Three-dimensional seismic wave propagation within a heterogeneous isotropic poroelastic medium is simulated with an explicit, time-domain, finite-difference algorithm. A system of thirteen, coupled, first-order partial differential equations is solved for the velocity vector components, stress tensor components, and pressure associated with solid and fluid constituents of the composite medium. A massively parallel computational implementation, utilizing the spatial domain decomposition strategy, allows investigation of large-scale earth models and/or broadband wave propagation within reasonable execution times.
An improved capability for subsurface structure detection is needed to support military and nonproliferation requirements for inspection and for surveillance of activities of threatening nations. As part of the DOE/NN-20 program to apply geophysical methods to detect and characterize underground facilities, Sandia National Laboratories (SNL) initiated an electromagnetic induction (EMI) project to evaluate low frequency electromagnetic (EM) techniques for subsurface structure detection. Low frequency, in this case, extended from kilohertz to hundreds of kilohertz. An EMI survey procedure had already been developed for borehole imaging of coal seams and had successfully been applied in a surface mode to detect a drug smuggling tunnel. The SNL project has focused on building upon the success of that procedure and applying it to surface and low altitude airborne platforms. Part of SNL`s work has focused on improving that technology through improved hardware and data processing. The improved hardware development has been performed utilizing Laboratory Directed Research and Development (LDRD) funding. In addition, SNL`s effort focused on: (1) improvements in modeling of the basic geophysics of the illuminating electromagnetic field and its coupling to the underground target (partially funded using LDRD funds) and (2) development of techniques for phase-based and multi-frequency processing and spatial processing to support subsurface target detection and characterization. The products of this project are: (1) an evaluation of an improved EM gradiometer, (2) an improved gradiometer concept for possible future development, (3) an improved modeling capability, (4) demonstration of an EM wave migration method for target recognition, and a demonstration that the technology is capable of detecting targets to depths exceeding 25 meters.
Borehole-to-surface electromagnetic (EM) methods are an attractive alternative to Surface-based EM methods for a variety of environmental and engineering applications. They have improved sensitivity to the subsurface resistivity distribution because of the closer proximity to the area of interest offered by the borehole for the source or the receiver. For the borehole-to-surface measurements the source is in the borehole and the receivers are on the surface. On the other hand, for the surface-to-borehole methods, the source is on the surface and the receiver is in a borehole. The surface-to-borehole method has an added advantage since measurements are often more accurate due to the lower noise environment for the receiver. For these methods, the source can be a grounded electric dipole or a vertical magnetic dipole source. An added benefit of these techniques is field measurements are made using a variety of arrays where the system is tailored to the application and where one can take advantage of some new imaging methods. In this short paper the authors describe the application of the borehole-to-surface method, discuss benefits and shortcomings, and give two field examples where they have been used for underground imaging. The examples were the monitoring of a salt water flooding of an oil well and the characterization of a fuel oil spill.
The multifrequency, multisource holographic method used in the analysis of seismic data is to extended electromagnetic (EM) data within the audio frequency range. The method is applied to the secondary magnetic fields produced by a borehole, vertical electric source (VES). The holographic method is a numerical reconstruction procedure based on the double focusing principle for both the source array and the receiver array. The approach used here is to Fourier transform the constructed image from frequency space to time space and set time equal to zero. The image is formed when the in-phase part (real part) is a maximum or the out-of-phase (imaginary part) is a minimum; i.e., the EM wave is phase coherent at its origination. In the application here the secondary magnetic fields are treated as scattered fields. In the numerical reconstruction, the seismic analog of the wave vector is used; i.e., the imaginary part of the actual wave vector is ignored. The multifrequency, multisource holographic method is applied to calculated model data and to actual field data acquired to map a diesel fuel oil spill.