Monitoring CO2 Storage at Cranfield, Mississippi with Time-Lapse Offset VSP – Using Integration and Modeling to Reduce Uncertainty

AbstractA time-lapse Offset Vertical Seismic Profile (OVSP) data set was acquired as part of a subsurface monitoring program for geologic sequestration of CO2. The storage site at Cranfield, near Natchez, Mississippi, is part of a detailed area study (DAS) site for geologic carbon sequestration operated by the U.S. Dept. of Energy's Southeast Regional Carbon Sequestration Partnership (SECARB). The DAS site includes three boreholes, an injection well and two monitoring wells. The project team selected the DAS site to examine CO2 sequestration multiphase fluid flow and pressure at the interwell scale in a brine reservoir. The time-lapse (TL) OVSP was part of an integrated monitoring program that included well logs, crosswell seismic, electrical resistance tomography and 4D surface seismic. The goals of the OVSP were to detect the CO2 induced change in seismic response, give information about the spatial distribution of CO2 near the injection well and to help tie the high-resolution borehole monitoring to the 4D surface data.The VSP data were acquired in well CFU 31-F1, which is the ∼3200 m deep CO2 injection well at the DAS site. A preinjection survey was recorded in late 2009 with injection beginning in December 2009, and a post injection survey was conducted in Nov 2010 following injection of about 250 kT of CO2. The sensor array for both surveys was a 50-level, 3-component, Sercel MaxiWave system with 15 m (49ft) spacing between levels. The source for both surveys was an accelerated weight drop, with different source trucks used for the two surveys.Consistent time-lapse processing was applied to both data sets. Time-lapse processing generated difference corridor stacks to investigate CO2 induced reflection amplitude changes from each source point. Corridor stacks were used for amplitude analysis to maximize the signal-to-noise ratio (S/N) for each shot point. Spatial variation in reflectivity (used to ‘map’ the plume) was similar in magnitude to the corridor stacks but, due to relatively lower S/N, the results were less consistent and more sensitive to processing and therefore are not presented. We examined the overall time-lapse repeatability of the OVSP data using three methods, the NRMS and Predictability (Pred) measures of Kragh and Christie (2002) and the signal-to-distortion ratio (SDR) method of Cantillo (2011). Because time-lapse noise was comparable to the observed change, multiple methods were used to analyze data reliability.The reflections from the top and base reservoir were identified on the corridor stacks by correlation with a synthetic response generated from the well logs. A consistent change in the corridor stack amplitudes from pre- to post-CO2 injection was found for both the top and base reservoir reflections on all ten shot locations analyzed. In addition to the well-log synthetic response, a finite-difference elastic wave propagation model was built based on rock/fluid properties obtained from well logs, with CO2 induced changes guided by time-lapse crosswell seismic tomography (Ajo-Franklin et al., 2013) acquired at the DAS site. Time-lapse seismic tomography indicated that two reservoir zones were affected by the flood. The modeling established that interpretation of the VSP trough and peak event amplitudes as reflectivity from the top and bottom of reservoir is appropriate even with possible tuning effects. Importantly, this top/base change gives confidence in an interpretation that these changes arise from within the reservoir, not from bounding lithology. The modeled time-lapse change and the observed field data change from 10 shotpoints are in agreement for both magnitude and polarity of amplitude change for top and base of reservoir. Therefore, we conclude the stored CO2 has been successfully detected and, furthermore, the observed seismic reflection change can be applied to Cranfield's 4D surface seismic for spatially delineating the CO2/brine interface

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 non-smooth spatial process. Timelapse 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 timelapse feature can be directly constrained. We develop a differential traveltime tomography algorithm which selects for compact solutions i.e. models with a minimum area of support, through application of model-space iteratively reweighted least squares. Our technique is an adaptation of minimum support regularization methods previously explored within the potential theory community. We compare our inversion algorithm to the results obtained by traditional Tikhonov regularization for two simple synthetic models; one including several sharp localized anomalies and a second with smoother features. We use a more complicated synthetic test case based on multiphase flow results to illustrate the efficacy of compactness constraints for contaminant infiltration imaging. We conclude by applying the algorithm to a CO[subscript 2] sequestration monitoring dataset acquired at the Frio pilot site. We observe that in cases where the assumption of a localized anomaly is correct, the addition of compactness constraints improves image quality by reducing tomographic artifacts and spatial smearing of target features.Massachusetts Institute of Technology. Earth Resources Laborator

Preseismic Velocity Changes Observed from Active Source Monitoringat the Parkfield SAFOD Drill Site

Measuring stress changes within seismically active fault zones has been a long-sought goal of seismology. Here we show that such stress changes are measurable by exploiting the stress dependence of seismic wave speed from an active source cross-well experiment conducted at the SAFOD drill site. Over a two-month period we observed an excellent anti-correlation between changes in the time required for an S 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 traveltime data that are coincident with two earthquakes that are among those predicted to produce the largest coseismic stress changes at SAFOD. Interestingly, 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

Continuous active-source seismic monitoring of CO2 injection in abrine aquifer

Continuous crosswell seismic monitoring of a small-scale CO2injection was accomplished with the development of a noveltubing-deployed piezoelectric borehole source. This piezotube source wasdeployed on the CO2 injection tubing, near the top of the saline aquiferreservoir at 1657-m depth, and allowed acquisition of crosswellrecordings at 15-minute intervals during the multiday injection. Thechange in traveltime recorded at various depths in a nearby observationwell allowed hour-by-hour monitoring of the growing CO2 plume via theinduced seismic velocity change. Traveltime changes of 0.2 to 1.0 ms ( upto 8 percent ) were observed, with no change seen at control sensorsplaced above the reservoir. The traveltime measurements indicate that theCO2 plume reached the top of the reservoir sand before reaching theobservation well, where regular fluid sampling was occuring during theinjection, thus providing information about the in situ buoyancy ofCO2

Reconstructing sea-level change in the Falkland Islands (Islas Malvinas) using salt-marsh foraminifera, diatoms and testate amoebae

Proxy records of past sea-level change provide a means of extending sea-level histories from tide gauges into the pre-industrial period. This is especially valuable in the South Atlantic region where sea-level data are limited to only a few tide-gauge records. Multi-proxy approaches to sea-level reconstruction are relatively rare but have distinct benefits when groups of micro-organisms are sparse or under-represented in modern or fossil sediments. Here, we address this challenge by utilising surface foraminifera, testate amoebae and diatoms from a salt marsh at Swan Inlet, East Falkland. All three micro-organism groups occupied distinct vertical niches in the contemporary salt-marsh. We investigated the relative performance of each group of micro-organisms in providing a sea-level reconstruction using individual (group-specific) regression models and with a multi-proxy regression model that combined all three groups. Foraminifera alone were not a suitable proxy. Surveyed sample elevations were closely matched by estimated elevations using Weighted-Average (WA) and Weighted-Average Partial-Least-Squares (WA- PLS) regressions. Relative sea-level reconstructions were derived by applying each model to microfossil assemblages recovered from a core (SI-2) from the same site. The combined transfer function yielded reconstructive precision (± 0.08 m) comparable to our best single-proxy transfer function (± 0.06 m) but only 18% of palaeo-samples were identified as having “close” or “good” analogues in the combined training data set. We highlight the benefit of a pragmatic approach to sea-level reconstructions whereby additional proxies should be employed if the use of only one proxy performs poorly across the width of the elevation gradient

Renewal processes and fluctuation analysis of molecular motor stepping

We model the dynamics of a processive or rotary molecular motor using a renewal processes, in line with the work initiated by Svoboda, Mitra and Block. We apply a functional technique to compute different types of multiple-time correlation functions of the renewal process, which have applications to bead-assay experiments performed both with processive molecular motors, such as myosin V and kinesin, and rotary motors, such as F1-ATPase

Cavity QED with a Bose-Einstein condensate

Cavity quantum electrodynamics (cavity QED) describes the coherent interaction between matter and an electromagnetic field confined within a resonator structure, and is providing a useful platform for developing concepts in quantum information processing. By using high-quality resonators, a strong coupling regime can be reached experimentally in which atoms coherently exchange a photon with a single light-field mode many times before dissipation sets in. This has led to fundamental studies with both microwave and optical resonators. To meet the challenges posed by quantum state engineering and quantum information processing, recent experiments have focused on laser cooling and trapping of atoms inside an optical cavity. However, the tremendous degree of control over atomic gases achieved with Bose-Einstein condensation has so far not been used for cavity QED. Here we achieve the strong coupling of a Bose-Einstein condensate to the quantized field of an ultrahigh-finesse optical cavity and present a measurement of its eigenenergy spectrum. This is a conceptually new regime of cavity QED, in which all atoms occupy a single mode of a matter-wave field and couple identically to the light field, sharing a single excitation. This opens possibilities ranging from quantum communication to a wealth of new phenomena that can be expected in the many-body physics of quantum gases with cavity-mediated interactions.Comment: 6 pages, 4 figures; version accepted for publication in Nature; updated Fig. 4; changed atom numbers due to new calibratio

Late-Holocene climate dynamics recorded in the peat bogs of Tierra del Fuego, South America

The ombrotrophic peat bogs of Tierra del Fuego are located within the southern westerly wind belt (SWWB), which dominates climate variability in this region. We have reconstructed late-Holocene water-table depths from three peat bogs and aimed to relate these records to shifts in regional climate. Water-table depths were quantified by the analysis of testate amoeba assemblages, and a regional transfer function was used to infer past water-table depths. During the last 2000 years, testate amoeba assemblages have been relatively stable, with a dominance of Difflugia pulex and Difflugia pristis type, and an increase in Assulina muscorum and other Euglyphida at the top of each section. Multivariate analyses show that water-table depth remained the main environmental variable explaining assemblages along the TiA12 core, but reconstructions were not significant for the two other cores. In line with the low variability in assemblages, water tables were relatively stable during the last 2000 years. Slightly wetter conditions were found between ~1400 and 900 cal. BP and a pronounced recent dry shift was reconstructed in all of the three peat profiles. Considering the regional climatic context, this recent shift may have been forced by a decrease in precipitation and warmer conditions linked to an increase in the importance of the SWWB. Nevertheless, we cannot exclude the influence of higher UV-B radiation resulting from the local degradation of the ozone layer since the late 1970s, which may have had an additional effect on the relative presence of A. muscorum in the southern Patagonian region

An integrated approach to seismic stimulation of oil reservoirs: laboratory, field and theoretical results from DOE/industry collaborations.

It has been observed repeatedly that low-frequency (10-500 Hz) seismic stress waves can enhance oil production from depleted reservoirs . Until recently, the majority of these observations have been anecdotal or at the proof-of-concept level. The physics coupling stress waves to multiphase fluid flow behavior in porous media is still poorly understood, even though numerous underlying physical mechanisms have been proposed to explain the observations . Basic research on the phenomenon is being conducted through a U .S. Department of Energy funded collaboration between Lawrence Berkeley National Laboratory, the University of California at Berkeley, Los Alamos National Laboratory and the U .S . oil and gas industry . The project has focused on three main areas of research: (1) laboratory core flow experiments, (2) field seismic monitoring of downhole stimulation tests, and (3) theoretical modeling of the coupled stress/flow phenomenon . The major goal is to obtain a comprehensive scientific understanding of the seismic stimulation phenomenon so that field application technologies can be improved. Initial developments and experimental results in all three research focus areas confirm historic observations that the stimulated flow phenomenon is real and that a fundamental scientific understanding can be obtained through continued research . Examples of project results and developments are presented here