Frictional properties affect the propagation of high-amplitude seismic waves across rock fractures. Laboratory evidence suggests that these properties can be measured in active seismic surveys, potentially offering a route to characterizing friction in situ. We present experimental results from a subresonance torsional modulus and attenuation apparatus that utilizes micron-scale sinusoidal oscillations to probe the nonlinear stress-strain relation at a range of strain amplitudes and rates. Nonlinear effects are further quantified using harmonic distortion; however, time series data best illuminate underlying physical processes. The low-frequency stress-strain hysteretic loops show stiffening at the sinusoid’s static ends, but stiffening is reduced above a threshold frequency. This shape is determined by harmonic generation in the strain; the stress signal has no harmonics, confirming that the fractured sample is the source of the nonlinearity. These qualitative observations suggest the presence of rate-dependent friction and are consistent between fractures in three different rock types. We propose that static friction at the low strain rate part of the cycle, when given sufficient “healing” time at low oscillation frequencies, causes this stiffening cusp shape in the hysteresis loop. While rate-and-state friction is commonly used to represent dynamic friction, it cannot capture static friction or negative slip velocities. So, we implement a bristle friction model, which describes this process and produces similar results.
The sparsity of permanent seismic instrumentation in marine environments often limits the availability of subsea information on geohazards, including active fault systems, in both time and space. One sensing resource that may provide observational access to the seafloor environment are existing networks of ocean bottom fiber optic cables; these cables, coupled to modern distributed acoustic sensing (DAS) systems, can provide dense arrays of broadband seismic observations capable of recording both seismic events and the ambient noise wavefield. Here, we report the detailed analysis of the ambient seismic noise acquired using DAS on a 20 km section of a fiber optic cable offshore of Moss Landing, CA, in Monterey Bay. Using this dataset, initially discussed in Lindsey et al. 2019, we extract Scholte waves using ambient noise interferometry techniques and invert the resulting multimodal dispersion curves to recover a high resolution 2D shear-wave velocity image of the near seafloor sediments. We show for the first time that the migration of coherently scattered Scholte waves observed on DAS records can provide an approach for resolving sharp lateral contrasts in subsurface properties, particularly shallow faults and depositional features near the seafloor. Our results provide improved constraints on shallow submarine features in Monterey Bay, including fault zones and paleo-channel deposits, thus highlighting one of many possible geophysical uses of the marine cable network.
Continuous crosswell seismic monitoring of a small scale CO2 injection was accomplished with the development of a novel tubing‐deployed piezoelectric borehole source. This 'piezotube' source was deployed on CO2 injection tubing, near the top of a saline aquifer reservoir at 1657 m depth, and allowed acquisition of crosswell recordings at 15‐minute intervals during the multi‐day injection. The change in travel time recorded at various depths in a nearby observation well allowed hour‐by‐hour monitoring of the growing CO2 plume via the induced seismic velocity change. Travel time changes of 0.2 to 1.0 ms (up to 8%) were observed, with no change seen at 'control' sensors placed above the reservoir. The travel time measurements indicate that the CO2 plume reached the top of the reservoir sand before reaching the observation well, where regular fluid sampling was occurring during the injection, thus providing information about the in‐situ buoyancy of CO2
The detailed mechanisms of the sealing of a single fracture, from hydration to almost complete closure by increase of confining pressure, as monitored from in situ synchrotron X-ray microtomography during the flow of carbonated water, is here shown for the first time. Different mechanisms play the key role at different stages in the evolution of the fracture. Hydration mechanically weakens the surfaces of the fracture and induces a first closure due to microcracking at the asperity contacts, increasing their size and creating choke points. Increase in confining stress promptly hydraulically seals the fracture by closing the main choke point, with a relative small deformation of the sample. Finally, the more pervasive mechanical deformation observed at higher stresses almost completely seals the whole fracture. The evolution of the sample has been described and quantified using 4D image processing, focusing on the evolution of aperture and digital volume correlation. Hydraulic properties of the sample at different stages have been modeled via Stokes flow simulation, and the results compared to the morphometric analysis, finding positive correlations with the average fracture aperture variation along the flowpath in function of time. Opalinus Clay is found to be a rock markedly prone to sealing in case of flow with carbonated water; this behavior is the result of the large fraction of clays and of its microstructure, lacking both cementing phases and large stiff particles. CO2 in this sample has no evident role in the evolution of the fracture; chemically-induced weathering on the surface has not been detected, in contrast with the behavior observed in samples with carbonates as cementing phase.
Summary We present results from an intermediate scale field experiment conducted in December 2014, evaluating the use of Distributed Acoustic Sensing (DAS) for recording ambient noise generated by infrastructure use; the future aim of our effort is to utilize ambient noise for real-time geotechnical evaluation of soils beneath active infrastructure. We present details of the installation, recording comparisons of 4 selected DAS cables, and a noise comparison between DAS and classical 3C geophones. We conclude by showing preliminary examples of ambient noise interferometry using passive data acquired with DAS vs. 3C geophones.
The Southeast Regional Carbon Sequestration Partnership (SECARB) early project in western Mississippi has been testing monitoring tools and approaches to document storage efficiency and storage permanence under conditions of CO2 EOR as well as downdip injection into brine. Denbury Onshore LLC is host for the study and has brought a depleted oil and gas reservoir, Cranfield Field, under CO2 flood. Injection was started in July 2008 and has now achieved injection rates greater than 1.2 million tons/year though 23 wells, with cumulative mass injected as of August, 2010 of 2.2 million metric tons. Injection is into coarse grained fluvial deposits of the Cretaceous lower Tuscaloosa Formation in a gentle anticline at depths of 3300 m. A team of researchers from 10 institutions has collected data from five study areas, each with a different goal and different spatial and temporal scale. The Phase 2 study began at the start of injection and has been using pressure and temperature as a tool for assessing permanence mostly in the oil productive interval. Real-time read-out shows high sensitivity to distant changes in injection rate and confirms the geologic model of reservoir compartmentalization. Above-zone pressure monitoring ∼120 m above the injection interval is used to test the sensitivity of this approach for documentation of integrity of the confining system in an area of numerous well completions as pressure increase is induced in the reservoir by more than 70 bar. Monitoring of the High Volume Injection Test (HiVIT) area includes repeat measurements of aqueous geochemistry in the injection zone. Rock-water- CO2 interactions in the reservoir as CO2 dissolves are minimized by mineral "armoring" by abundant chlorite cement in high permeability reservoir sandstone. Geochemical monitoring of confined freshwater aquifers at depths of 70–100 m is underway. Groundwater analysis focuses on assessment of the sensitivity of this method to detect leakage above background variability. A repeat seismic survey of the HiVIT is planned for late 2010 to assess saturation change especially in downdip brine-only areas. A study focused on feasibility of monitoring the shallow subsurface to separate leakage from normal complex surface fluxes is underway at an monitoring array installed in October 2009 to assess the interactions of recharge, soil gas, and shallow groundwater aquifers. Recent well re-entry and tracer injection will provide further information to interpret observed elevated deep-sourced methane. The Detailed Area Study (DAS) is collecting dense time-lapse data from closely-spaced three well array of an injector and two observation wells. The observation wells were completed with fiberglass casing to facilitate electrical resistance tomography (ERT) measurements, and a diverse array of instrumentation was both cemented behind casing and suspended on tubing. Injection started at the DAS December 1, 2009. We have measured pulsed neutron and resistivity via wireline, downhole and above-zone pressure, distributed temperature, and fluid chemistry including introduced pulses of perfluorocarbons, noble gases, and SF6 as tracers. Between wells, time-lapse crosswell seismic and electrical resistance tomography (ERT) are used to measure saturation change. The goals are to measure changes as fluids evolve from single phase (brine) to two phase (CO2–brine) in order to document linkages between pressure and sweep efficiency. A time-lapse VSP survey bridges the vertical resolution and areal coverage between cross-well and surface seismic. The repeat surveys for many tools are scheduled for September, 2010. Reservoir characterization based on cores, historic and new wireline log data, production history, hydrologic tests, fluid analysis, and a three-D seismic survey have been used in multiple numerical models to predict reservoir response in order to design effective monitoring strategies and optimize deployment. History matching of observed response to predicted response is used to interpret results and improve confidence in conceptual models and numerical approaches. Probabilistic methods have been used to assess the significant uncertainties resulting from reservoir heterogeneity. This study is funded by the US Department of Energy, National Energy Technology Laboratory as part of the Regional Carbon Sequestration Partnerships program. SECARB is led by Southern States Energy Board.
Hydrogeophysical methods are presented that support the siting and monitoring of aquifer storage and recovery (ASR) systems. These methods are presented as numerical simulations in the context of a proposed ASR experiment in Kuwait, although the techniques are applicable to numerous ASR projects. Bulk geophysical properties are calculated directly from ASR flow and solute transport simulations using standard petrophysical relationships and are used to simulate the dynamic geophysical response to ASR. This strategy provides a quantitative framework for determining site-specific geophysical methods and data acquisition geometries that can provide the most useful information about the ASR implementation. An axisymmetric, coupled fluid flow and solute transport model simulates injection, storage, and withdrawal of fresh water (salinity ∼500 ppm) into the Dammam aquifer, a tertiary carbonate formation with native salinity approximately 6000 ppm. Sensitivity of the flow simulations to the correlation length of aquifer heterogeneity, aquifer dispersivity, and hydraulic permeability of the confining layer are investigated. The geophysical response using electrical resistivity, time-domain electromagnetic (TEM), and seismic methods is computed at regular intervals during the ASR simulation to investigate the sensitivity of these different techniques to changes in subsurface properties. For the electrical and electromagnetic methods, fluid electric conductivity is derived from the modeled salinity and is combined with an assumed porosity model to compute a bulk electrical resistivity structure. The seismic response is computed from the porosity model and changes in effective stress due to fluid pressure variations during injection/recovery, while changes in fluid properties are introduced through Gassmann fluid substitution.
