Thomas M. Daley

Seismic low-frequency effects in monitoring fluid-saturated reservoirs

There is a complex relationship between seismic attributes, including the frequency dependence of reflections and fluid saturation in a reservoir. Observations in both laboratory and field data indicate that reflections from a fluid-saturated layer have an increased amplitude and delayed traveltime at low frequencies, when compared with reflections from a gas-saturated layer. Comparison of laboratory-modeling results with a diffusive-viscous-theory model show that low (<5) values of the quality factor Q can explain the observations of frequency dependence. At the field scale, conventional processing of time-lapse VSP data found minimal changes in seismic response of a gas-storage reservoir when the reservoir fluid changed from gas to water. Lowfrequency analysis found significant seismic-reflectionattribute variation in the range of 15‐50 Hz. The field observations agree with effects seen in laboratory data and predicted by the diffusive-viscous theory. One explanation is that very low values of Q are the result of internal diffusive losses caused by fluid flow. This explanation needs further theoretical investigation. The frequencydependent amplitude and phase-reflection properties presented in this paper can be used for detecting and monitoring fluid-saturated layers.

Preseismic velocity changes observed from active source monitoring at the Parkfield SAFOD drill site

Measuring stress changes within seismically active fault zones has been a long-sought goal of seismology. One approach is to exploit the stress dependence of seismic wave velocity, and we have investigated this in an active source cross-well experiment at the San Andreas Fault Observatory at Depth (SAFOD) drill site. Here we show that stress changes are indeed measurable using this technique. Over a two-month period, we observed an excellent anti-correlation between changes in the time required for a shear wave to travel through the rock along a fixed pathway (a few microseconds) and variations in barometric pressure. We also observed two large excursions in the travel-time data that are coincident with two earthquakes that are among those predicted to produce the largest coseismic stress changes at SAFOD. The two excursions started approximately 10 and 2 hours before the events, respectively, suggesting that they may be related to pre-rupture stress induced changes in crack properties, as observed in early laboratory studies.

Applying compactness constraints to differential traveltime tomography

Tomographic imaging problems are typically ill-posed and often require the use of regularization techniques to guarantee a stable solution. Minimization of a weighted norm of model length is one commonly used secondary constraint. Tikhonov methods exploit low-order differential operators to select for solutions that are small, flat, or smooth in one or more dimensions. This class of regularizing functionals may not always be appropriate, particularly in cases where the anomaly being imaged is generated by a nonsmooth spatial process. Time-lapse imaging of flow-induced velocity anomalies is one such case; flow features are often characterized by spatial compactness or connectivity. By performing inversions on differenced arrival time data, the properties of the time-lapse feature can be directly constrained. We develop a differential traveltime tomography algorithm whichselects for compact solutions, i.e., models with a minimum area of support, through application of model-space iteratively reweighted leas...

Fracture detection using crosswell and single well surveys

We recorded high‐resolution (1 to 10 kHz), crosswell and single well seismic data in a shallow (15 to 35 m), water‐saturated, fractured limestone sequence at Conocos borehole test facility near Newkirk, Oklahoma. Our objective was to develop seismic methodologies for imaging gas‐filled fractures in naturally fractured gas reservoirs. The crosswell (1/4 m receiver spacing, 50 to 100 m well separation) surveys used a piezoelectric source and hydrophones before, during, and after an air injection that we designed to displace water from a fracture zone. Our intent was to increase the visibility of the fracture zone to seismic imaging and to confirm previous hydrologic data that indicated a preferred pathway. For the single well seismic imaging (a piezoelectric source and an eight‐element hydrophone array at 1/4 m spacing), we also recorded data before and after the air injection. The crosswell results indicate that the air did follow a preferred pathway that was predicted by hydrologic modeling. In addition,...

Laboratory seismic monitoring of supercritical CO2 flooding in sandstone cores using the Split Hopkinson Resonant Bar technique with concurrent x-ray CT imaging

ABSTRACT Accurate estimation of CO2 saturation in a saline aquifer is essential for the monitoring of supercritical CO2 injected for geological sequestration. Because of strong contrasts in density and elastic properties between brine and CO2 at reservoir conditions, seismic methods are among the most commonly employed techniques for this purpose. However the relationship between seismic (P‐wave) velocity and CO2 saturation is not unique because the velocity depends on both wave frequency and the CO2 distribution in rock. In the laboratory, we conducted measurements of seismic properties of sandstones during supercritical CO2 injection. Seismic responses of small sandstone cores were measured at frequencies near 1 kHz, using a modified resonant bar technique (Split Hopkinson Resonant Bar method). Concurrently, saturation and distribution of supercritical CO2 in the rock cores were determined via x‐ray CT scans. Changes in the determined velocities generally agreed with the Gassmann model. However, both the velocity and attenuation of the extension wave (Youngs modulus or ‘bar’ wave) for the same CO2 saturation exhibited differences between the CO2 injection test and the subsequent brine re‐injection test, which was consistent with the differences in the CO2 distribution within the cores. Also, a comparison to ultrasonic velocity measurements on a bedded reservoir rock sample revealed that both compressional and shear velocities (and moduli) were strongly dispersive when the rock was saturated with brine. Further, large decreases in the velocities of saturated samples indicated strong sensitivity of the rocks frame stiffness to pore fluid.

Monitoring CO2 Intrusion and Associated Geochemical Transformations in a Shallow Groundwater System Using Complex Electrical Methods

The risk of CO(2) leakage from a properly permitted deep geologic storage facility is expected to be very low. However, if leakage occurs it could potentially impact potable groundwater quality. Dissolved CO(2) in groundwater decreases pH, which can mobilize naturally occurring trace metals commonly contained in aquifer sediments. Observing such processes requires adequate monitoring strategies. Here, we use laboratory and field experiments to explore the sensitivity of time-lapse complex resistivity responses for remotely monitoring dissolved CO(2) distribution and geochemical transformations that may impact groundwater quality. Results show that electrical resistivity and phase responses correlate well with dissolved CO(2) injection processes. Specifically, resistivity initially decreases due to increase of bicarbonate and dissolved species. As pH continues to decrease, the resistivity rebounds toward initial conditions due to the transition of bicarbonate into nondissociated carbonic acid, which reduces the total concentration of dissociated species and thus the water conductivity. An electrical phase decrease is also observed, which is interpreted to be driven by the decrease of surface charge density as well as potential mineral dissolution and ion exchange. Both laboratory and field experiments demonstrate the potential of field complex resistivity method for remotely monitoring changes in groundwater quality due to CO(2) leakage.

Crosswell seismic imaging in a contaminated basalt aquifer

Multiple seismic crosswell surveys have been acquired and analyzed in a fractured basalt aquifer at Idaho National Engineering and Environmental Laboratory. Most of these surveys used a high‐frequency (1000–10,000 Hz) piezoelectric seismic source to obtain P‐wave velocity tomograms. The P‐wave velocities range from less than 3200 m/s to more than 5000 m/s. Additionally, a new type of borehole seismic source was deployed as part of the subsurface characterization program at this contaminated groundwater site. This source, known as an orbital vibrator, allows simultaneous acquisition of P‐ and S‐waves at frequencies of 100 to 400 Hz, and acquisition over larger distances. The velocity tomograms show a relationship to contaminant transport in the groundwater; zones of high contaminant concentration are coincident with zones of low velocity and high attenuation and are interpreted to be fracture zones at the boundaries between basalt flows. The orbital vibrator data show high Vp/Vs values, from 1.8 to 2.8. In...

Orbital vibrator seismic source for simultaneous P‐ and S‐wave crosswell acquisition

A recently developed borehole seismic source, the orbital vibrator, was successfully deployed in a crosswell survey in a fractured basalt aquifer. This seismic source uses a rotating eccentric mass to generate seismic energy. Source sweeps with clockwise and counter‐clockwise rotations are recorded at each source location. Because this source generates circularly polarized waves, unique processing algorithms are used to decompose the recordings into two equivalent linearly oscillating, orthogonally oriented seismic sources. The orbital vibrator therefore generates P‐ and S‐waves simultaneously for all azimuths. A coordinate rotation based on P‐wave particle motion is used to align the source components from various depths. In a field experiment, both P‐ and S‐wave arrivals were recorded using fluid‐coupled hydrophone sensors. The processed field data show clear separation of P‐ and S‐wave arrivals for in‐line and crossline source components, respectively. A tensor convolutional description of the decompos...

Cost-effective imaging of CO2 injection with borehole seismic methods

Currently there is a critical need to increase oil and gas recovery from existing and new reservoirs. In addition, the ever-increasing need to sequester CO2 in the subsurface places further emphasis on accurate imaging methods to validate CO2 injection strategies. Two obstacles to increased efficiency are (1) a thorough understanding of the geologic complexity and fluid distribution and (2) the scaling relationships between fine scale/point measurements and larger scale/volumetric measurements. Although initially expensive, borehole methods may offer a cost-effective solution when integrated into a drilling and development program. New technology such as fiber-optic sensors emplaced during drilling and completion, microhole drilling, and other advances in sensors will make borehole technology much more cost effective when used over the long run. If deployed in a multicomponent and time-lapse fashion, seismic methods also offer the ability to define contrasts in properties, detecting subtle changes in prop...

Tube‐wave suppression in single‐well seismic acquisition

Single‐well seismic imaging is significantly hampered by the presence of borehole tube waves. A tube‐wave suppressor has been tested using single‐well seismic equipment at the Lost Hills (California) oil field. The suppressor uses a gas‐filled bladder kept slightly above borehole fluid pressure. Field tests show a measurable reduction in tube‐wave energy as compared to body waves propagating in the surrounding reservoir rock. When using a high‐frequency (500–4000 Hz) piezoelectric source, the P‐wave–tube‐wave amplitude ratio was increased by 33 dB. When using a lower frequency (50–350 Hz) orbital vibrator source, the S‐wave–tube‐wave amplitude ratio was increased by 21 dB while the P‐wave–tube‐wave amplitude ratio was increased by 23 dB. These reductions in tube‐wave amplitudes significantly improve single‐well data quality.