Frequency-Domain Modeling Techniques for the Scalar Wave Equation : An Introduction

Frequency-domain finite-difference (FDFD) modeling offers several advantages over traditional timedomain methods when simulating seismic wave propagation, including a convenient formulation within the context of wavefield inversion and a straight-forward extension for adding complex attenuation mechanisms. In this short paper we introduce the FDFD method, develop a simple solver for the scalar Helmholtz problem, and explore some possible approaches for solving large scale seismic modeling problems in the frequency domain.Massachusetts Institute of Technology. Earth Resources Laborator

The effect of CO2-induced dissolution on flow properties in Indiana Limestone: An in situ synchrotron X-ray micro-tomography study

The injection of CO -rich fluids in carbonate rocks results in an evolution of the pore space, with consequent changes in the hydraulic properties of the reservoir; how these properties evolve, particularly for parameters relevant to multiphase flow e.g. P (s), remains a topic of active research despite several decades of study. We have carried out an in situ synchrotron X-ray microtomography experiment to monitor pore structure evolution during dissolution of an Indiana Limestone core; the experiment involved flowing CO -saturated water through the core for 36 h and resulted in 10 volumes corresponding to different temporal stages of the dissolution process. The injection parameters corresponded to the flow velocities expected near the well-bore region of a shallow aqueous CO injection; fast flow rates with high reactant availability. Analysis of the tomographic data shows flow-enhanced dissolution i.e. channeling, and provides a time-resolved map of pore space alteration. Using the resulting 4D pore space volume, we modeled the evolution of capillary-pressure curves; this exercise demonstrates how pore structure evolution could impact the invasion and remobilization of non-wetting fluids, dramatically decreasing the entry pressure and the P in some parts of the sample. The modeling of permeability, using a Stokes solver approach, quantified the relationship of porosity vs. permeability; we found that a modest increase in porosity, especially when the channeling system is more developed, greatly affects permeability. These results demonstrate how movement of CO saturated brine near injected plumes might alter drainage dynamics near the plume boundary, thus leading to mobilization across subtle capillary barriers. 2 c 2 2 C

Extracellular electron uptake from a cathode by the lactic acid bacterium Lactiplantibacillus plantarum

A subset of microorganisms that perform respiration can endogenously utilize insoluble electron donors, such as Fe(II) or a cathode, in a process called extracellular electron transfer (EET). However, it is unknown whether similar endogenous EET can be performed by primarily fermentative species like lactic acid bacteria. We report for the first time electron uptake from a cathode by Lactiplantibacillus plantarum, a primarily fermentative bacteria found in the gut of mammals and in fermented foods. L. plantarum consumed electrons from a cathode and coupled this oxidation to the reduction of both an endogenous organic (pyruvate) and an exogenous inorganic electron acceptor (nitrate). This electron uptake from a cathode reroutes glucose fermentation toward lactate degradation and provides cells with a higher viability upon sugar exhaustion. Moreover, the associated genes and cofactors indicate that this activity is mechanistically different from that one employed by lactic acid bacteria to reduce an anode and to perform respiration. Our results expand our knowledge of the diversity of electroactive species and of the metabolic and bioenergetic strategies used by lactic acid bacteria

Non-Linear Constraints with Application to Self-Potential Source Inversion

We investigate the use of non-linear constraints for geophysical inverse problems, with specific examples applied to source inversion of self-potential data. Typical regularization methods often produce smooth solutions by introducing a quadratic term in the objective function that minimizes the L2 norm of a low-order differential operator applied to the model. In some cases, however, the properties of interest may not vary smoothly. Two alternative constraints are examined that provide inversion stability while allowing for solutions with non-smooth properties. One method, often referred to as ‘compactness’ or ‘minimum support’, seeks to minimize the area (in 2D) or volume (in 3D) occupied by non-zero model parameters. The second method, ‘total variation’, minimizes an approximation of the L1 norm of the gradient of the model. Both approaches involve a non-linear regularization functional, and must therefore be solved iteratively. We discuss the practical aspects of implementing these regularization methods and compare several examples using self-potential source inversion on a synthetic model. We also apply the compactness constraint for self-potential source inversion using a field data example.Kuwait-MIT Center for Natural Resources and the EnvironmentMassachusetts Institute of Technology. Earth Resources Laborator

Temporal Integration of Seismic Traveltime Tomography

Time-lapse geophysical measurements and seismic imaging methods in particular are powerful techniques for monitoring changes in reservoir properties. Traditional time-lapse processing methods treat each dataset as an independent unit and estimate changes in reservoir state through differencing these separate inversions. We present a general least-squares approach to jointly inverting time-varying property models through use of spatio-temporal coupling operators. Originally developed within the medical imaging community, this extension of traditional Tikhonov regularization allows us to constrain the way in which models vary in time, thereby reducing artifacts observed in traditional time-lapse imaging formulations. The same methodology can also accommodate changes in experiment geometry as a function of time thus allowing inversion of incremental or incomplete surveys. In this case, temporal resolution is traded for improved spatial coverage at individual timesteps. We use seismic traveltime tomography as a model problem although almost any geophysical inversion task can be posed within this formalism. We apply the developed time-lapse inversion algorithm to a synthetic crosswell dataset designed to replicate a CO2 sequestration monitoring experiment

A new mini-triaxial cell for combined high-pressure and high-temperature in situ synchrotron X-ray microtomography experiments up to 400°C and 24 MPa.

A new experimental triaxial cell for in situ synchrotron X-ray micro-computed tomography aimed at imaging small samples of (6 mm × 19 mm) at high temperatures (up to 400°C) and pressures (up to 24 MPa confining) is presented. The system has flow-through capabilities, independent axial and radial pressure control, and has been developed and tested at the 8.3.2. beamline at the Advanced Light Source. The characteristics of this new experimental rig are described, along with the challenges, mainly concerning the combination of X-ray transparency with vessel strength at high temperature, and solutions found during the development stage. An experiment involving oil shale pyrolysis under subsurface conditions, highlighting the importance of a device able to operate in this pressure and temperature range, is also introduced. The availability of this cell enables an unprecedented range of experiments in the Earth Sciences, with a special focus on subsurface geothermal processes

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

Applying Compactness Constraints to Seismic 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 seismic velocity anomalies is one such case; flow features are often characterized by spatial compactness or connectivity. We develop a traveltime tomography algorithm which selects for compact solutions through application of model-space iteratively reweighted least squares. Our technique is an adaptation of minimum support regularization methods previously developed within the potential theory community. We emphasize the application of compactness constraints to timelapse datasets differenced in the data domain, a process which allows recovery of compact perturbations in model properties. We test our inversion algorithm on a simple synthetic dataset generated using a velocity model with several localized velocity anomalies. We then demonstrate the efficacy of the algorithm on a CO2 sequestration monitoring dataset acquired at the Frio pilot site. In both cases, the addition of compactness constraints improves image quality by reducing spatial smearing due to limited angular aperture in the acquisition geometry.Toksoz, M. NafiMassachusetts Institute of Technology. Earth Resources Laborator

The Potential of DAS in Teleseismic Studies: Insights From the Goldstone Experiment

Distributed acoustic sensing (DAS) is a recently developed technique that has demonstrated its utility in the oil and gas industry. Here we demonstrate the potential of DAS in teleseismic studies using the Goldstone OpticaL Fiber Seismic experiment in Goldstone, California. By analyzing teleseismic waveforms from the 10 January 2018 M7.5 Honduras earthquake recorded on ~5,000 DAS channels and the nearby broadband station GSC, we first compute receiver functions for DAS channels using the vertical-component GSC velocity as an approximation for the incident source wavelet. The Moho P-to-s conversions are clearly visible on DAS receiver functions. We then derive meter-scale arrival time measurements along the entire 20-km-long array. We are also able to measure path-averaged Rayleigh wave group velocity and local Rayleigh wave phase velocity. The latter, however, has large uncertainties. Our study suggests that DAS will likely play an important role in many fields of passive seismology in the near future

The role of stress and fluid saturation on the acoustic response of fractured rock

Standard rock physics models are formulated to describe the behavior of porous sedimentary reservoirs, with clean sandstones being the archetypal system; however, many situations demand geophysical monitoring of rocks with significantly different structures, such as low porosity, fractured reservoirs. Conventional models also suggest that these “stiff” reservoirs can be challenging to monitor seismically due to small fluid substitution effects, but the presence of fractures leads to stress dependence which may be leveraged for remote monitoring purposes. Using samples from the Duperow Formation (dolostone) obtained from the Danielson test well in Kevin Dome, MT, we conducted ultrasonic and multi-scale structural (profilometry, synchrotron micro-tomography, pressure sensitive film) measurements on naturally fractured core in order to characterize the effects of fluid substitution and effective stress on the acoustic response of fractured reservoir rock with a focus in particular on the textural and seismic characteristics of natural fractures. We find that changes in effective stress can yield changes in velocity of up to 20% and changes in attenuation up to 200%. Measured fluid substitution effects are resolvable, but stress effects dominate. These measurements provide insight into the physical processes controlling acoustic response of fractured rocks in general and can also be used to inform monitoring efforts in fractured reservoirs