Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Project seismic”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

Analysis of Microseismicity and Reactivated Fault Size to Assess the Potential for Felt Events by CO 2 Injection in the Illinois Basin

The results of monitoring of carbon dioxide (CO 2 ) injection at the Illinois Basin—Decatur Project (IBDP) and the companion Illinois Industrial Carbon Capture and Sequestration Sources (IL-ICCS) project—have shown that reservoir response to fluid pressure changes can vary significantly at different injection locations within the same reservoir. Predrill reservoir characterization is important to identify potentially seismogenic faults. However, interpretations of newly reprocessed 3D seismic reflection data illustrate the challenges related to their identification in a region dominated by faulting with small vertical offsets. Faults interpreted in the 3D seismic volume range from ~300 to 1200 m wide and are in the same size range as faults that could have been the source of historical events up to Mw 2.7 in central Illinois. The array of monitoring sensors that was installed for the IBDP continues to collect data, as injection operates in IL-ICCS, the second injection well. CO 2 injection rates for the IL-ICCS well are on average 1.7 times the rates injected in the IBDP well, but a significantly reduced rate of induced seismicity is observed. This article presents results of passive seismic monitoring for the duration of the project to date, integrating active and passive seismic data to develop a new interpretation of the subsurface structure at the Decatur site that explicitly identifies pathways for fluid flow into the basement leading to induced seismicity, and provides a geological explanation for the sharp reduction of induced seismicity during injection at higher rates into the second well. The use of seismic moment to estimate the length of seismogenic slip planes in the local subsurface suggests that faults large enough to produce felt seismicity are unlikely to be present at or near the Decatur site.

Geochemistry & Geophysics↗

Utah FORGE 5-2419: Seismicity-Permeability Relationships Probed via Nonlinear Acoustic Imaging - Workshop Presentation

This is a presentation on the Seismicity-Permeability Relationships Probed via Nonlinear Acoustic Imaging project by Pennsylvania State University, presented by Derek Elsworth. The project's objectives were to explore controls and acoustic signatures of aseismic through seismic evolution of friction-stability-permeability relationships on fractures in shear-reactivation, link this reactivation to key features of the pre-existing stress state, and upscale these indexes to reservoir scale to drive successful reservoir stimulation. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 8, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗

Large-Volume Stimulation of Rock for Greatly Enhanced Fluids Recovery using Targeted Seismic-Assisted Hydraulic Fracturing (Final Technical Report)

This project has developed and demonstrated a new technology for large-volume and targeted comminution of rock in low permeability formations to enhance recovery from unconventional oil and gas (UOG) resources. The technology is based on a strategically designed interaction of multiple induced seismic pulses that assist the hydraulic fracturing process to enhance shear and multi-planar crack formation. This greatly increased rock stimulation, through bulk comminution, is expected to cause significant increase in permeability leading to enhancement of recovery factors for sub- surface fluids. The proposed technology is especially applicable for enhanced recovery in emerging UOG plays, such as ductile shales that are resistant to opening-mode fracturing by conventional hydraulic fracturing processes. The project combines an integrated experimental and computation approach to develop and demonstrate a modular technology that can be easily implement in the field to augment current practices. The effort integrates a fundamental scientific understanding of dynamic material response under constraint and damage-induced permeability and porosity enhancements at multiple length scales, along with models of comminution due to the local release of kinetic energy associated with high shear strain rate of dynamic deformation. The results will be validated through small-scale experiments and then implemented in a lab-scale field test to demonstrate the developed technology. The project involves collaboration between a solid mechanician, a materials scientist, and a petroleum engineering geo-mechanician. It builds on a successful history of success and collaboration that includes both fundamental scientific exploration and technology development.

04 OIL SHALES AND TAR SANDS↗

Identification of Faults Susceptible to Induced Seismicity Final Report

The primary objective of this project (referred to as the State of Stress project or SoS project) was to predict the presence of faults susceptible to movement in the presence of fluid injection as a consequence of vertical pressure migration from the storage unit to the crystalline basement. These faults included those difficult to resolve with traditional interpretations of surface seismic: faults with no or minimal offset, faults in the crystalline basement where no seismic reflectors are present, subseismic faults (due to size), and faults undetected (e.g., due to unfavorable orientation with surface seismic survey).

25 ENERGY STORAGE↗

Oregon Rivers Characterization Tool

DHS’s Regional Resiliency Assessment Program (RRAP) undertook the Oregon State Transportation Systems RRAP project between 2018-2021 to assess the seismic vulnerabilities of the state’s transportation system to a Cascadia Subduction Zone (CSZ) earthquake (CISA 2021). The project primarily focused on the effects of a CSZ earthquake on the mobility of emergency response efforts into the region. In the event that bridges located at crucial river crossings fail during a CSZ earthquake, rivers will become natural barriers affecting the ability to support response and recovery efforts. The characteristics of rivers (e.g., width, depth, and flow) are critical data to inform planning for the recovery of transportation systems across river barriers and determine reopening times following a CSZ earthquake. As part of the transportation system assessment, an Argonne National Laboratory (Argonne) project team developed a River Characterization Tool (RCT) to quantify the physical characteristics of the river system in Oregon, which can be used to inform a system-level analysis of transportation recovery and reopening times. This report mainly discusses approaches, methods, development, and applications of the RCT. The use of the river data generated from this tool for a system-level analysis of transportation recovery is discussed in the Oregon Transportation Systems RRAP Resiliency Assessment report (CISA 2021). The primary goal of the RCT is to generate statewide dataset of widths, depths, and flows for all rivers across Oregon to assess the effects of river barriers on emergency response and recovery efforts following a CSZ earthquake. The RCT leverages remote sensing data that multispectral sensors collect from satellites, employs a deep learning model to characterize river predictions, and computes river channel geometric dimensions with hydraulic relationships and river predictions. The entire process of river characterization is automated and implemented through the RCT.

58 GEOSCIENCES↗

2021 Summer SPE Project

SPE is a project in order to develop new, more physics-based, seismic models of explosions (see Nelson et al. 2013). One key component of this effort is numerical modeling enabled by modern state-of-art code/software. Accurate modeling of the shape and amplitude of seismic waves from their generation to their propagation to remote monitoring seismic stations is important to our ability to determine the origin and strength of the source from remote recording. In this project, the modeling is performed by coupling two codes HOSS and SPECFEM3D. HOSS models the dynamic nonlinear processes happening near the explosion. SPECFEM3D computes the propagation of seismic waves as they travel through 3D complex Earth models. SPECFEM3D was modified in order to be driven by a set of time-series calculated by HOSS in lieu of a seismic source. The goal of this summer project is the investigation of several questions pertinent to the establishment of a full end-to-end modeling capability from the high-rate strain regime area to remote distances where seismic station record seismic waves generated by explosions. The investigated questions are: (1) How to perform proper filtering? Direct-solution modeling only sustains a limited range of frequency depending on the grid size. As we go from one modeling domain to the other via coupling, the mesh size is getting coarser to allow modeling at large scale but also to account for the fact that high-frequency waves do not physically travel to large distances. So filtering of the time-series generated by the near-field hydrodynamic modeling is a current practice often employed but its effects on the modeled waveforms has to be investigated. (2) Quantitative assessment of the efficiency of attenuation to remove high-frequency content of the wavefield. This assessment will allow to create meshes with a grid size appropriate to the actual physics of wave propagation for a given explosion. (3) Checking that the coupling process and the two codes respect the cylindrical symmetry that is expected in the case of a pure explosion in a half-space. (4) The effect of the state-of-stress in the near-source area on the modeled seismic waveforms. These questions will be investigated through the modeling of SPE-4P, the fourth explosion of this series because it was designed to have little interaction with the geologic surrounding and the free surface of the Earth so that it is the most explosion-like experiment, with the most symmetries to be verified.

58 GEOSCIENCES↗

Imperial Valley Dark Fiber Project Continuous DAS Data

The Imperial Valley Dark Fiber Project acquired Distributed Acoustic Sensing (DAS) seismic data on a ~28 km segment of dark fiber between the cities of Calipatria and Imperial in the Imperial Valley, Southern California. Dark fiber refers to unused optical fiber cables in telecommunications networks and is repurposed in this project for DAS applications. The objective, which is further detailed in the attached journal article from Ajo-Franklin et al., is to demonstrate dark fiber DAS as a tool for basin-scale geothermal exploration and monitoring. The included DAS data were recorded during two days at the beginning the project. Data is stored in the .h5 (HDF5) file format, readable using various software tools, including the 'h5read' and 'h5info' functions in Matlab. Provided here are examples of methods to read and use the data with the 'h5py' package in Python.

15 GEOTHERMAL ENERGY↗

Seismic monitoring of CO 2 geosequestration using multi-well 4D DAS VSP: Stage 3 of the CO2CRC Otway project

An important part of any CO 2 geosequestration project is to ensure CO 2 containment and conformance in the subsurface. This is generally done by implementing a comprehensive, risk-based Measurement, Monitoring and Verification plan, a key element of which is active time-lapse seismic monitoring. However, high cost and environmental impact of the standard surface seismic monitoring dictate the need for a cost-effective and environmentally friendly alternative. An opportunity to develop such method emerges with advances in distributed acoustic sensing (DAS) technology, which turns an optical fibre into a seismic sensor with dense spatial sampling. DAS can be permanently deployed in multiple wells across the geosequestration site providing a robust and non-intrusive network of seismic receivers. This approach was developed and tested in the CO2CRC Otway project, where injection of 15 kt of CO 2 at 1.5 km depth was monitored with a 4D vertical seismic profiling (VSP) using five borehole DAS arrays and mobile vibroseis sources. The 4D DAS VSP in each of the five wells provides broadly consistent images of the CO 2 plume with some differences due to different illumination of the target horizon, lateral variation of velocities, and seismic anisotropy. When the newly injected CO 2 reaches a CO 2 plume created as a result of an earlier injection into the same formation ~600 m updip, 4D DAS VSP shows a change in reflectivity in that area and beyond. Furthermore, this shows a potential of 4D DAS VSP for monitoring gas injection into gas-saturated reservoirs.

58 GEOSCIENCES↗

Seismic response to paleo-sand dunes in the Nugget Sandstone Formation, southwestern Wyoming

We have analyzed a 3D seismic survey acquired for a carbon sequestration project on top of the Moxa Arch in southwestern Wyoming. We observed a zone of discontinuous reflectors on vertical slices of seismic amplitude volume, whereas, the northwest–southeast lineations were observed on the time slices. We performed a seismic to well tie that suggested that the lineations occur within the Nugget Sandstone. The Nugget Sandstone is an eolian sandstone deposit of Early Jurassic age, deposited as a subtropical dune field, and equivalent to the Navajo Sandstone of southwestern Utah. Petrophysical analysis indicates that the Nugget Sandstone is dominated by clean sandstone (70%–80%), whereas evaporites, including halite and anhydrite, are present in certain zones. Previous outcrop studies on the Navajo Sandstone indicate the wind direction to be predominantly northeast–southwest. Seismic attributes, including coherence and curvature, displayed on stratal slices within the Nugget Sandstone interval indicate the presence of lineations in the northwest–southeast direction with irregular spacing. These lineations are approximately perpendicular to the inferred dominant wind direction. We computed the dominant wind direction from the average azimuth of the lineations as seen on the curvature attribute in the Nugget Sandstone interval. Geological feature: Eolian sand dunes with interdunal evaporites Seismic appearance: Parallel lineations with irregular spacing on seismic attribute horizon slices Alternative interpretations: Canyons at continental slopes; slope failures Features with a similar appearance: Marine bars; contourites Formation: The Nugget Sandstone — equivalent to the Navajo Sandstone Age: Early Jurassic Location: Moxa Arch, Wyoming Seismic data: Obtained by the University of Wyoming with U.S. DOE funding Contributors: Dhruv Agrawal, Brady Lujan, Sumit Verma, Shuvajit Bhattacharya, and Subhashis Mallick Analysis tools: Coherence and curvature attributes; seismic inversion; petrophysical inversion

Geochemistry & Geophysics↗

Connect the Dots: In Situ 4-D Seismic Monitoring of CO 2 Storage With Spatio-Temporal CNNs

4-D seismic imaging has been widely used in CO 2 sequestration projects to monitor the fluid flow in the volumetric subsurface region that is not sampled by wells. Ideally, real-time monitoring and near-future forecasting would provide site operators with great insights to understand the dynamics of the subsurface reservoir and assess any potential risks. However, due to obstacles such as high deployment cost, availability of acquisition equipment, exclusion zones around surface structures, only very sparse seismic imaging data can be obtained during monitoring. That leads to an unavoidable and growing knowledge gap over time. The operator needs to understand the fluid flow throughout the project lifetime and the seismic data are only available at a limited number of times. This is insufficient for understanding reservoir behavior. To overcome those challenges, we have developed spatio-temporal neural-network-based models that can produce high-fidelity interpolated or extrapolated images effectively and efficiently. Specifically, our models are built on an autoencoder, and incorporate the long short-term memory (LSTM) structure with a new loss function regularized by optical flow. We validate the performance of our models using real 4-D post-stack seismic imaging data acquired at the Sleipner CO 2 sequestration field. We employ two different strategies in evaluating our models. Numerically, we compare our models with different baseline approaches using classic pixel-based metrics. We also conduct a blind survey and collect a total of 20 responses from domain experts to evaluate the quality of data generated by our models. Finally, via both numerical and expert evaluation, we conclude that our models can produce high-quality 2-D/3-D seismic imaging data at a reasonable cost, offering the possibility of real-time monitoring or even near-future forecasting of the CO 2 storage reservoir.

4-D seismic imaging↗

Multi-Scale Seismic Measurements for Site Characterization and CO2 Monitoring in an Enhanced Oil Recovery/Carbon Capture, Utilization, and Sequestration Project, Farnsworth Field, Texas

To address the challenges of climate change, significantly more geologic carbon sequestration projects are beginning. The characterization of the subsurface and the migration of the plume of supercritical carbon dioxide are two elements of carbon sequestration that can be addressed through the use of the available seismic methods in the oil and gas industry. In an enhanced oil recovery site in Farnsworth, TX, we employed three separate seismic techniques. The three-dimensional (3D) surface seismic survey required significant planning, design, and processing, but produces both a better understanding of the subsurface structure and a three-dimensional velocity model, which is essential for the second technique, a timelapse vertical seismic profile, and the third technique, cross-well seismic tomography. The timelapse 3D Vertical Seismic Profile (3D VSP) revealed both significant changes in the reservoir between the second and third surveys and geo-bodies that may represent the extent of the underground carbon dioxide. The asymmetry of the primary geo-body may indicate the preferential migration of the carbon dioxide. The third technique, cross-well seismic tomography, suggested a strong correlation between the well logs and the tomographic velocities, but did not observe changes in the injection interval.

Energy & Fuels↗

Wabash CarbonSAFE. Geologic Analysis of the Potosi Dolomite Reservoir Interval and Potential Confining Units, at the Wabash #1 Well, Vigo Co., Indiana, Subtask 7.2, Technical Report

The Wabash CarbonSAFE project drilled the Wabash #1 stratigraphic test well (ID# 168045) at the Wabash Valley Resources (WVR) IGCC facility in Vigo County, Indiana, to evaluate the feasibility of commercial-scale CO 2 storage near the site. The Mt. Simon Sandstone was the initial target for storage evaluation and was found to have generally poor reservoir qualities in the Wabash #1 well. Focus was then placed on the secondary target, the Potosi Dolomite. For storage and containment evaluation an extensive suite of geophysical logs was collected throughout the borehole, and whole cores were obtained from the Mt. Simon Sandstone, Eau Claire Formation, and Maquoketa Group. The Maquoketa Group is identified as a regionally extensive secondary sealing interval above the Mt. Simon Sandstone and as a primary seal for the Potosi Dolomite storage complex. Due to concerns over potential lost circulation zones in the Potosi Dolomite during drilling, no core samples or Formation Micro Imager (FMI) logs were acquired in the Potosi Dolomite. However, cased hole well testing was performed over a 20 feet (6 m) interval within the Potosi Dolomite during which a fluid swab sample was obtained. Approximately 35 miles of 2D seismic information was acquired in the project area near Wabash #1 to aid in reservoir and caprock characterization. The highly porous and permeable zones in the Potosi Dolomite and lower Oneota Dolomite observed within the Wabash #1 well are, for the purpose of this report, considered the same reservoir package and referred to as the Potosi reservoir interval. This report presents a regional depositional context for the Potosi reservoir interval, the overlying confining zones in the Knox Group, and Maquoketa Group primary seal. This context is developed through detailed information obtained from the Wabash #1 well such as lithologic data collected from cuttings and core, geophysical logging, geomechanical analysis of core samples, and well testing and fluid sampling within the Potosi reservoir interval. Information discussed in this report was compiled as a basis for the development of static and dynamic geologic models of the Potosi Dolomite storage complex. The static and dynamic modeling of CO 2 injection in the Potosi reservoir interval are discussed in a separate report (Dessenberger et al., 2022) under the Wabash CarbonSAFE project.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Gas hydrate quantification in Walker Ridge block 313, Gulf of Mexico, from full-waveform inversion of ocean-bottom seismic data

The Gulf of Mexico (GOM) Joint Industry Project Leg 2 logging-while-drilling data in Walker Ridge lease block 313 (WR313) in the GOM detected gas hydrate in coarse- and fine-grained sediments at sites WR313-G and WR313-H. The coarse-grained units are thin ([Formula: see text]) and highly saturated, whereas the fine-grained unit is thick (approximately 200 m) with low saturation and fracture-filling gas hydrate. Unlike its coarse-grained counterpart, the seismic character of the fine-grained unit does not clearly indicate the presence of gas hydrate, which would likely have remained undiscovered in the absence of drilling. In this paper, through frequency-domain acoustic full-waveform inversion (FWI) of ocean-bottom seismometer data along a 2D multichannel seismic transect near sites WR313-G and WR313-H, we detect and quantify gas hydrate in the fine-grained unit. Key results are as follows: First, the base of the gas hydrate stability zone, which is not obvious in the reflection profile, can be discerned in the FWI results. Second, the gas hydrate in the fine-grained unit is mainly confined to the area between two sets of opposite-dipping normal faults implying that the fault architecture may be partially responsible for this gas hydrate accumulation and distribution.

Geochemistry & Geophysics↗

Understanding the Seismic Ground Motion Spatial Variability Using Network Analysis Community Detection

This project is to explore ground motion spatial distribution using a new approach graph-based network analysis. In this study, we combine a large-N seismic array and graph analytics to explore spatial variability and correlation at a local scale using small local and regional earthquakes. In this method, each seismic station is modeled as a node and the similarities of the waveforms that represent ground motions between two stations are modeled as edges. By analyzing this graph network using the similarity matrices and community detection algorithm, we can group the stations spatially with similar patterns. A random forest algorithm is used to reveal the important features that affect the spatial grouping. The result suggests site conditions, and how they interact with the incident seismic wavefield, strongly condition the spatial correlation of ground motion. Future progress in characterizing ground motion spatial variability will require dense wavefield measurements, either through nodal deployments, or perhaps distributed acoustic sensing measurements of seismic wavefields.

58 GEOSCIENCES↗

Integration of seismic-pressure-petrophysics inversion of continuous active-seismic monitoring data for monitoring and quantifying CO 2 plume (Final Report)

The overall objective of this project is to develop and validate an integrated package of joint seismic-pressure-petrophysics inversion (jSPPI) of continuous active-source seismic monitoring dataset capable of providing real-time monitoring of CO 2 plume during geologic carbon sequestration (GCS). The three specific developments include: (a) the methodologies for fast seismic full waveform inversion of continuous active source seismic monitoring, (CASSM) datasets for simultaneously estimating velocity and attenuation, and with data assimilation; (b) joint Bayesian petrophysical inversion of seismic models and pressure data for providing and updating CO 2 saturation models; (c) the methods using multiple datasets including (Crainfield and Frio-II borehole) synthetic, laboratory, and field CASSM datasets. The outcomes of jSPPI include (a) a workflow for processing CASSM data, (b) Bayesian inversion algorithms using CASSM data and pressure response data, and (c) integration with data assimilation algorithms for continuously updating site-specific models used for prediction and reservoir management. The validation of joint FWI will be conducted using synthetic models based on the Cranfield and Frio experiments as well as field CASSM datasets collected as part of the Frio-II pilot injection. To quantify and map the mass and distribution of CO 2 (saturation), we will jointly invert velocity and attenuation measurements from the FWI with a Bayesian approach using a rock physics model for attenuation (e.g., White’s attenuation model with two selected patch sizes (White, 1976; Dutta and Seriff, 1979)). The Bayesian inversion will be applied to each time step in the CASSM survey in an updating scheme, which integrates with an ensemble of reservoir simulations at each step. A more complete experimental validation dataset will be collected as part of a mesoscale (2-3 m) gas-CO 2 injection experiment utilizing a higher frequency version of the CASSM system developed for laboratory studies; the integrated inversion will be demonstrated using this dataset which will provide both a dense geometry as well as more precise secondary confirmation measurements (e.g. saturation) typically not available in the field. The resulting real-time map of CO 2 saturation is able to provide a deeper scientific understanding of the complex, time-varying dynamics of subsurface fluid flow migration path as well as the rapid detection of CO 2 leakage hazards.

25 ENERGY STORAGE↗

Monitoring spatiotemporal evolution of fractures during hydraulic stimulations at the first EGS collab testbed using anisotropic elastic-waveform inversion

The EGS Collab project acquired continuous active-source seismic monitoring (CASSM) data before, during, and after hydraulic stimulations at the first testbed at the depth of 4850 ft (1478 m) at the Sanford Underground Research Facility in Lead, South Dakota, for monitoring fracture creation and evolution. CASSM acquisition was conducted using 24 hydrophones, 18 accelerometers, and 17 piezoelectric sources within four fracture-parallel wells and two orthogonal wells. 3D anisotropic traveltime tomography and anisotropic elastic-waveform inversion of the campaign cross-borehole seismic data show that the rock within the stimulation region is a heterogeneous horizontal transverse isotropic medium. Here we use these inversion results as the initial models and apply 3D anisotropic first-arrival traveltime tomography and 3D anisotropic elastic-waveform inversion to the CASSM data acquired after each stimulation in May, 2018 and December, 2018. We observe the spatiotemporal evolution of seismic velocities and anisotropic parameters caused by hydraulic fracture stimulations, showing the regions of rock alternation caused by hydraulic fracture stimulation.

15 GEOTHERMAL ENERGY↗

GraphAlign: Graph-Enabled Machine Learning for Seismic Event Filtering

This report summarizes results from a 2 year effort to improve the current automated seismic event processing system by leveraging machine learning models that can operated over the inherent graph data structure of a seismic sensor network. Specifically, the GraphAlign project seeks to utilize prior information on which stations are more likely to detect signals originating from particular geographic regions to inform event filtering. To date, the GraphAlign team has developed a Graphical Neural Network (GNN) model to filter out false events generated by the Global Associator (GA) algorithm. The algorithm operates directly on waveform data that has been associated to an event by building a variable sized graph of station waveforms nodes with edge relations to an event location node. This builds off of previous work where random forest models were used to do the same task using hand crafted features. The GNN model performance was analyzed using an 8 week IMS/IDC dataset, and it was demonstrated that the GNN outperforms the random forest baseline. We provide additional error analysis of which events the GNN model performs well and poorly against concluded by future directions for improvements.

58 GEOSCIENCES↗