Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic migration”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Potential Seismicity Along Basement Faults Induced by Geological Carbon Sequestration

Large-scale CO 2 sequestration into geological formations has been suggested to reduce CO 2 emissions from industrial activities. However, much like enhanced geothermal stimulation and wastewater injection, CO 2 sequestration has a potential to induce earthquake along weak faults, which can be considered a negative impact on safety and public opinion. This research shows the physical mechanisms of potential seismic hazards along basement faults driven by CO 2 sequestration under variation in geological and operational constraints. Specifically we compare the poroelastic behaviors between multiphase flow and single-phase flow cases, highlighting specific needs of evaluating induced seismicity associated with CO 2 sequestration. In contrast to single-phase injection scenario, slower migration of the CO 2 plume than pressure pulse may delay accumulation of pressure and stress along basement faults that may not be mitigated immediately by shut-in of injection. The impact of multiphase flow system, therefore, needs to be considered for proper monitoring and mitigation strategies.

58 GEOSCIENCES↗
DOE OSTI2021

Imaging Complex Subsurface Structures for Geothermal Exploration at Pirouette Mountain and Eleven-Mile Canyon in Nevada

Accurate imaging of subsurface complex structures with faults is crucial for geothermal exploration because faults are generally the primary conduit of hydrothermal flow. It is very challenging to image geothermal exploration areas because of complex geologic structures with various faults and noisy surface seismic data with strong and coherent ground-roll noise. In addition, fracture zones and most geologic formations behave as anisotropic media for seismic-wave propagation. Properly suppressing ground-roll noise and accounting for subsurface anisotropic properties are essential for high-resolution imaging of subsurface structures and faults for geothermal exploration. We develop a novel wavenumber-adaptive bandpass filter to suppress the ground-roll noise without affecting useful seismic signals. This filter adaptively exploits both characteristics of the lower frequency and the smaller velocity of the ground-roll noise than those of the signals. Consequently, this filter can effectively differentiate the ground-roll noise from the signal. We use our novel filter to attenuate the ground-roll noise in seismic data along five survey lines acquired by the U.S. Navy Geothermal Program Office at Pirouette Mountain and Eleven-Mile Canyon in Nevada, United States. We then apply our novel anisotropic least-squares reverse-time migration algorithm to the resulting data for imaging subsurface structures at the Pirouette Mountain and Eleven-Mile Canyon geothermal exploration areas. The migration method employs an efficient implicit wavefield-separation scheme to reduce image artifacts and improve the image quality. Our results demonstrate that our wavenumber-adaptive bandpass filtering method successfully suppresses the strong and coherent ground-roll noise in the land seismic data, and our anisotropic least-squares reverse-time migration produces high-resolution subsurface images of Pirouette Mountain and Eleven-Mile Canyon, facilitating accurate fault interpretation for geothermal exploration.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

A New Approach to Predict Hydrogeological Parameters Using Shear Waves from the Multichannel Analysis of Surface Waves Method

For near-surface contaminant characterization, the accurate prediction of hydrogeological parameters in anisotropic and heterogeneous environments has been a challenge since the last decades. However, recent advances in near-surface geophysics have facilitated the use of geophysical data for hydrogeological characterization in the last few years. A pseudo 3-D high resolution P-wave shallow seismic reflection survey was performed at the P Reactor Area, Savannah River Site, South Carolina in order to delineate and predict migration pathways of a large contaminant plume including trichloroethylene. This contaminant plume originates from the northwest section of the reactor facility that is located within the Upper Atlantic Coastal Plain. The data were collected with 40 Hz geophones, an accelerated weight-drop as seismic source and 1 m receiver spacing with near- and far-offsets of 0.5 and 119.5 m, respectively. In such areas with near-surface contaminants, a detailed subsurface characterization of the vadose zone hydraulic parameters is very important. Indeed, an inexpensive method of deriving such parameters by the use of seismic reflection surveys is beneficial, and our approach uses the relationship between seismic velocity and hydrogeological parameters together with empirical observations relating porosity to permeability and hydraulic conductivity. Shear wave velocity ( V s ) profiles were estimated from surface wave dispersion analysis of the seismic reflection data and were subsequently used to derive hydraulic parameters such as porosity, permeability, and hydraulic conductivity. Additional geophysical data including core samples, vertical seismic profiling, surface electrical resistivity tomography, natural gamma and electrical resistivity logs allowed for a robust assessment of the validity and geological significance of the estimated V s and hydrogeological models. The results demonstrate the usefulness of this approach for the upper 15 m of shallow unconsolidated sediments even though the survey design parameters were not optimal for surface wave analysis due to the higher than desired frequency geophones.

Engineering↗
DOE OSTI2024

Migration and accumulation of hydrous mantle incipient melt in the Earth's asthenosphere: Constraints from in-situ falling sphere viscometry measurements

The Earth's oceanic lithosphere-asthenosphere boundary (LAB) is marked by a notable decrease in seismic wave velocities and an increase in attenuation. This phenomenon is likely attributed to the accumulation of partial melt at the top of the asthenosphere. Nevertheless, the process involving the upward migration and aggregation of low-degree partial melts, highly likely to be mantle incipient melts in the asthenosphere, remains underexplored. Viscosity is a key factor controlling the flow of melts, thus in this study, we used In-situ X-ray falling sphere viscometry experiments to determine the viscosity of the mantle incipient melt containing 1.8–4.0 wt.% water at 1.5–6 GPa and 2100 K. We found that water and pressure can effectively decrease the viscosity of the mantle incipient melt. Therefore, the mantle incipient melts are highly mobile with a high segregation velocity in the deeper part of the asthenosphere. However, during ascent, their mobility diminishes due to an increase in viscosity and a decrease in water content. In conclusion, these mechanisms potentially contribute to the accumulation of partial melt at the top of the asthenosphere.

58 GEOSCIENCES↗
DOE OSTI2022

High-Precision Characterization of Seismicity from the 2022 Hunga Tonga-Hunga Ha'apai Volcanic Eruption

The earthquake swarm accompanying the January 2022 Hunga Tonga-Hunga Ha'apai (HTHH) volcanic eruption includes a large number of posteruptive moderate-magnitude seismic events and presents a unique opportunity to use remote monitoring methods to characterize and compare seismic activity with other historical caldera-forming eruptions. We compute improved epicentroid locations, magnitudes, and regional moment tensors of seismic events from this earthquake swarm using regional to teleseismic surface-wave cross correlation and waveform modeling. Precise relative locations of 91 seismic events derived from 59,047 intermediate-period Rayleigh- and Love-wave cross-correlation measurements collapse into a small area surrounding the volcano and exhibit a southeastern time-dependent migration. Regional moment tensors and observed waveforms indicate that these events have a similar mechanism and exhibit a strong positive compensated linear vector dipole component. Precise relative magnitudes agree with regional moment tensor moment magnitude ($M_w$) estimates while also showing that event sizes and frequency increase during the days after the eruption followed by a period of several weeks of less frequent seismicity of a similar size. The combined information from visual observation and early geologic models indicate that the observed seismicity may be the result of a complex series of events that occurred after the explosive eruption on 15 January, possibly involving rapid resupply of the magma chamber shortly after the eruption and additional faulting and instability in the following weeks. In addition, we identify and characterize an $M_w$ 4.5 event five days before the paroxysmal explosion on 15 January, indicating that additional seismic events preceding the main eruption could have been identified with improved local monitoring. As a result, our analysis of the HTHH eruption sequence demonstrates the value of potentially utilizing teleseismic surface-wave cross correlation and waveform modeling methods to assist in the detailed analysis of remote volcanic eruption sequences.

58 GEOSCIENCES↗
DOE OSTI2020

Imaging of a fluid injection process using geophysical data — A didactic example

In many subsurface industrial applications, fluids are injected into or withdrawn from a geologic formation. It is of practical interest to quantify precisely where, when, and by how much the injected fluid alters the state of the subsurface. Routine geophysical monitoring of such processes attempts to image the way that geophysical properties, such as seismic velocities or electrical conductivity, change through time and space and to then make qualitative inferences as to where the injected fluid has migrated. The more rigorous formulation of the time-lapse geophysical inverse problem forecasts how the subsurface evolves during the course of a fluid-injection application. Using time-lapse geophysical signals as the data to be matched, the model unknowns to be estimated are the multiphysics forward-modeling parameters controlling the fluid-injection process. Properly reproducing the geophysical signature of the flow process, subsequent simulations can predict the fluid migration and alteration in the subsurface. The dynamic nature of fluid-injection processes renders imaging problems more complex than conventional geophysical imaging for static targets. This work intents to clarify the related hydrogeophysical parameter estimation concepts.

58 GEOSCIENCES↗
DOE OSTI2026

Prairie State Generating Company 2D Seismic Interpretation

The objectives of the Illinois Storage Corridor (ISC) project are to accelerate commercial deployment of carbon capture utilization and storage at two individual sites and receive approvals for Underground Injection Control (UIC) Class VI permits for construction at each site (ISC Project Narrative, 2020). As part of this project, and as part of the subsurface geologic characterization, two-dimensional (2D) seismic data was acquired at both sites. This report summarizes the findings from two phases of 2D acquisition and seismic interpretation at the Prairie State Generating Company site near Marissa, Illinois. The seismic data indicates the presence of three faults that completely transect the storage and confining units. The initial 2021 2D seismic acquisition revealed the presence of a feature two miles east of the Lively Grove #1 characterization well. The feature is a polyphase fault with a component of strike-slip motion, forming a small positive flower structure. The second phase of 2D seismic acquisition in 2022 constrained the maximum extent of the fault, indicating it has a relatively limited length. Other 2022 seismic lines revealed the presence of two other faults that transect the storage and confining units in the southwest and northeast portions of the 2D seismic acquisition area. The 2021 and 2022 2D seismic surveys have identified three specific locations within the project area that may have an elevated risk of out of zone CO 2 migration due to faulting. Further technical work will be needed to quantify this risk and determine how and if this risk will impact the placement of CO 2 injection wells across the project area.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2026

Coupling geophysical, geological, geochemical and mineralogical assessments to examine preferential contaminant transport pathways in interbedded fractured bedrock

This study shows that a multi-faceted approach, combining borehole geophysical logging and surface seismic P-wave first-arrival tomography with confirmatory coring, well installation, and chemical and mineralogical analysis, is effective for identifying difficult-to-locate preferential contaminant transport pathways in deeper fractured bedrock. Seismic tomography detected porous 10–20 m wide elongated fractured conduits that allow acidic groundwater contaminated with uranium (U) and nitrate (NO 3 - ) to migrate within interbedded shale-limestone bedrock over 1000 m from a former disposal facility (S-3 Ponds site) located at the DOE Y-12 National Security Complex in Tennessee (USA). Conventional drilling techniques would easily miss these conduits because they are oriented parallel with fractured bedding planes. Synchrotron analysis of aquifer solids revealed that > 95 % of the U is hexavalent (U VI ). This uranyl (UO 2 2+ ) species is coordinated with carbonate, iron oxide, silicate and phosphate minerals within cemented white to yellow precipitates, which contain U concentrations as high as ∼21.6 % by weight fraction. Identifying the presence of these mineral phases, enables a further understanding of the potential effectiveness of remediation actions. The combination of methodologies presented here can also be applied to other explorations, such as the detection of water supply.

Fate and transport↗
DOE OSTI2024

Nowcasting Earthquakes With Stochastic Simulations: Information Entropy of Earthquake Catalogs

Earthquake nowcasting has been proposed as a means of tracking the change in large earthquake potential in a seismically active area. The method was developed using observable seismic data, in which probabilities of future large earthquakes can be computed using Receiver Operating Characteristic methods. Furthermore, analysis of the Shannon information content of the earthquake catalogs has been used to show that there is information contained in the catalogs, and that it can vary in time. So an important question remains, where does the information originate? In this paper, we examine this question using stochastic simulations of earthquake catalogs. Our catalog simulations are computed using an Earthquake Rescaled Aftershock Seismicity (“ERAS”) stochastic model. This model is similar in many ways to other stochastic seismicity simulations, but has the advantage that the model has only 2 free parameters to be set, one for the aftershock (Omori-Utsu) time decay, and one for the aftershock spatial migration away from the epicenter. Generating a simulation catalog and fitting the two parameters to the observed catalog such as California takes only a few minutes of wall clock time. While clustering can arise from random, Poisson statistics, we show that significant information in the simulation catalogs arises from the “non-Poisson” power-law aftershock clustering, implying that the practice of de-clustering observed catalogs may remove information that would otherwise be useful in forecasting and nowcasting. We also show that the nowcasting method provides similar results with the ERAS model as it does with observed seismicity.

58 GEOSCIENCES↗
DOE OSTI2024

Tracking Seismic Velocity Perturbations at Ridgecrest Using Ballistic Correlation Functions

Abstract We present results based on data of a dense nodal array composed of 147 stations, deployed in 2022 near the epicenter of the 2019 Mw 7.1 Ridgecrest earthquake to investigate characteristics of the seismic wavefields. Through array analyses, we identified two primary components. First, we observed far-field P waves dominating the 0.5–1.2 Hz frequency range, which are likely primarily generated by wind-driven oceanic swell activity. Second, we detected near-field body waves resulting from anthropogenic activities in the frequency range 2–8 Hz. We examined noise correlation functions derived from data of the dense deployment and regional stations to explore fault-zone seismic velocity changes using ballistic arrivals, with a focus on velocity perturbation shortly before and after the Ridgecrest earthquake sequence. Our findings exhibit distinct behavior compared to results obtained through standard coda-wave interferometry. Particularly, we observed a decrease in P-wave travel time on certain station pairs prior to the 2019 earthquake sequence. Supported by detailed investigation of the local seismic wavefields, we interpret the decreasing P-wave travel time as likely caused by a velocity increase away from the fault, possibly related to fluid migration. However, additional information is necessary to verify this hypothesis.

Geochemistry & Geophysics↗
DOE OSTI2022

Cretaceous sequence stratigraphy of the northern Baltimore Canyon Trough: Implications for basin evolution and carbon storage

We evaluate the Cretaceous stratigraphy and carbon sequestration potential of the northern Baltimore Canyon Trough (NBCT) using >10,000 km of multi-channel seismic profiles integrated with geophysical logs, biostratigraphy, and lithology from 29 offshore wells. We identify and map six sequences resolved primarily at the stage level. Accommodation was dominated by thermal and non-thermal subsidence, though sequence boundaries correlate with regional and global sea-level changes, and the record is modified by igneous intrusion, active faulting, and changes in sediment supply and sources. Our stratigraphic maps illustrate a primary southern (central Appalachian) Early Cretaceous source that migrated northward during the Aptian and Albian. During the Cenomanian, sedimentation rates in the NBCT increased and depocenters shifted northward and landward. We show that deposition occurred in three phases: (1) earlier Cretaceous paleoenvironments were primarily terrestrial indicated by variable amplitude, chaotic seismic facies, serrated gamma logs, and heterolithic sandstones and mudstones with terrestrial microfossils; (2) the Albian to Cenomanian was dominated by deltaic paleoenvironments indicated by blocky, funnel-shaped, gamma-ray logs and clinoforms characterized by continuous high-amplitude seismic reflections with well-defined terminations; and (3) the Cenomanian and younger was marine shelf, inferred from mudstoneprone lithologies, peak gamma-ray values in well logs, and foraminiferal evidence. Long-term transgression and maximum water depths at the Cenomanian/Turonian boundary correlative with Ocean Anoxic Event 2 were followed by a regression and relative sea-level fall. We show that porous and permeable sandstones of three Aptian to Cenomanian highstand systems tracts are high-volume reservoirs for supercritical CO 2 storage that are confined by overlying deep water mudstones.

58 GEOSCIENCES↗
DOE OSTI2024

Site Characterization of the Highest-Priority Geologic Formations for CO2 Storage in Wyoming

The project Site Characterization of the Highest-Priority Geologic Formations for CO2 Storage in Wyoming is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The Wyoming Carbon Underground Storage Project (WY-CUSP) consisted of CO2 storage site characterization and evaluation, focusing on Wyoming’s most promising CO2 storage reservoirs (the Pennsylvanian Weber/Tensleep Sandstone and Mississippian Madison Limestone) and premier CO2 storage site (Rock Springs Uplift). Results from the WY-CUSP project suggest the two reservoirs could store up to 17,000 million tons of CO2. The WY-CUSP team drilled a stratigraphic test well and acquired a 3-D seismic survey covering 25 square miles of the Rock Springs Uplift site. The team retrieved 916 feet of core from the 12,810-foot-deep well, along with a complete log suite, borehole images, fluid samples, and other data. Project partners (1) provided continuous visual documentation of the core, including grain size, mineralogy, facies distribution, and porosity; (2) performed continuous permeability and velocity scans of selected reservoir intervals; and (3) chemically analyzed the fluid samples. WY-CUSP scientists integrated seismic attributes with observations from log suites, a VSP survey, core, fluid samples, and laboratory analyses, including continuous permeability scans. From these integrations, researchers constructed 3-D spatial distribution volumes of reservoir and seal properties that represent geological heterogeneity at the targeted CO2 storage site. The WY-CUSP team used this data to perform new CO2 plume migration simulations. Baker Hughes, Inc., completed a series of small-scale, in-situ water injectivity measurements. A database was formed when observations, analyses, and experiments from the stratigraphic test well were integrated. Correlation of these data allowed petrophysical parameters to be extrapolated from the test well out into the storage domain (5x5 mile 3-D seismic survey volume). This resulted in an improved, realistic understanding of performance assessments for potential CO2 storage scenarios. The WY-CUSP team worked on (1) improving CO2 storage resource estimates, (2) establishing long-term integrity and permanence of confining layers, (3) designing a profitable strategy for pressure management, and (4) evaluating the utilization of stored CO2 at the Rock Spring Uplift. Finally, Baker Hughes developed a microseismic baseline for the test site using in-bore geophones to complete field operations.

3-D seismic↗
DOE OSTI2020

Integrated Modeling Report for CO 2 Storage with Enhanced Oil Recovery in Northern Michigan

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. MRCSP has been investigating various reservoir characterization and modeling technologies related to Carbon Capture, Utilization, and Storage (CCUS) in conjunction with carbon dioxide enhanced oil recovery (CO 2 -EOR) operations in multiple Silurian-age (Niagaran), oil-bearing, carbonate pinnacle reefs in northern Michigan, USA. This report provides a comprehensive discussion of reservoir modeling studies that were conducted for tracking oil production, forecasting CO 2 plume migration, and estimating associated storage in a number of reefs that were at different stages of their CO 2 -EOR life cycle. The modeling process for simulating oil production, CO 2 injection, and associated storage in these reefs entails two phases. The first phase, geologic framework modeling, integrates all pertinent geological and geophysical data (from logs, cores and seismic surveys) about reservoir structure, geometry, rock types, and property distributions (porosity, permeability, water saturation) into a 3-D distributed grid-based static earth model (SEM). The second phase, dynamic reservoir modeling, uses the SEM as a platform to simulate the movement of oil, gas, water, and CO 2 within the reservoir during primary hydrocarbon production, as well as during subsequent phases such as CO 2 -injection assisted EOR, plume migration, and associated storage. In addition, an assessment of coupled process effects is also carried out, where the impacts of geochemical and geo-mechanical processed induced by CO 2 injection are studied.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Preliminary Analysis of Source Physics Experiment Explosion-Triggered Microseismicity Using the Back-Projection Method

A series of four chemical explosions were detonated in a deep borehole within the Yucca Flat Dry Alluvium Geology (DAG) at the Nevada National Security Site between 2018 and 2019. The two larger chemical explosions of 50 tons (DAG-2) and 10 tons (DAG-4) TNT equivalent yield triggered energetic aftershock sequences numbering 1392 and 347 microearthquakes, respectively, within the first 10 days. No significant aftershock activity was observed for the two smaller 1-ton explosions (DAG-1 and DAG-3). Here, we used a back-projection method based on travel-time migration and stacking of signal-to-noise ratio traces to detect, associate and locate aftershocks from a subset of 22-geophones within a larger 2 urn:x-wiley:21699313:media:jgrb54910:jgrb54910-math-0001 2 km seismic array surrounding the borehole. The aftershocks located within 300 m of the borehole and the depths were above the working points of 300 and 50 m depths of DAG-2 and DAG-4, respectively, ruling out triggering slip on geologic faults or disturbances beneath neighboring collapse craters. DAG-2 and DAG-4 aftershocks decayed at similar rates, with power-law exponents of p = 1.48 and p = 1.49, respectively. These decay rates are comparable to aftershocks sequences triggered by earthquakes and historical nuclear explosions at Yucca Flat. A smooth power-law aftershock decay within the first 10 days suggests a triggering mechanism from explosion generated stress relaxation due to the diffusion of high gas pressures in the cavity and radial fractures. A more random and episodic aftershock rate would be expected due to cavity collapse or falling rubble in chimney formation.

58 GEOSCIENCES↗
DOE OSTI2021

Recommended Practices for Managing Induced Seismicity Risk Associated with Geologic Carbon Storage

The geologic storage of carbon dioxide (CO 2 ) is one method to help reduce or eliminate atmospheric CO 2 emissions. The sequestered CO 2 is originally captured from the atmosphere or from a stationary industrial source and subsequently injected into a deep subsurface porous rock formation. To facilitate the successful deployment of field scale carbon storage projects, the U.S. Department of Energy (DOE) is developing tools and protocols for defensible, science-based frameworks to quantify and mitigate risks associated with the long-term storage of CO 2 . This protocol specifically addresses the risk of induced seismicity due to injection in a geologic carbon storage (GCS) site. This integrated and risk-based protocol is a product of the U.S. DOE Fossil Energy’s National Risk Assessment Partnership (NRAP), a multi-year collaborative research effort of Los Alamos National Laboratory (LANL), Lawrence Berkeley National Laboratory (LBNL), Lawrence Livermore National Laboratory (LLNL), National Energy Technology Laboratory (NETL), and Pacific Northwest National Laboratory (PNNL). These recommended practices describe a set of 7 steps to evaluate, manage, communicate, and mitigate the risk of induced seismicity at GCS sites. The base methodology of the recommended practices follows a framework similar to the Protocol for Addressing Induced Seismicity Associated with Enhanced Geothermal Systems (Majer et al., 2012), developed for the Geothermal Technology Office of the U.S. DOE. These recommended practices present a framework to systematically assess the induced seismicity risk and quantify the associated uncertainties. These recommendations are based on current research and are sufficiently general to allow for modification and application to a variety of different types of sites. The substance of the recommended practices contained herein includes both technical and non-technical issues, and covers all operational stages of the GCS project lifecycle. They start at the preliminary risk assessment phase, continue through site assessment and characterization, include best practice communication and seismic monitoring plan methodologies, discuss the evaluation and mitigation of seismic hazard and risk, and closes with an exploration of operational management plans, which conclude when the induced seismicity risk abates back to background level. The focus of these recommendations is on actively managing the risks associated with induced seismicity by developing an actionable risk management plan that starts at the project proposal stage and continues through site closure through an iterative assessment and improvement process. The audience of this document is expected to include all interested stakeholders (e.g., operators, project developers, regulators, and the general public) and is expressly written to be accessible to this broad range of partners. This document is intended to disseminate knowledge gained through recent advances in the science of induced seismicity hazard and risk assessments, to provide updates based on recent experience gained by similar corollary injection-induced seismicity cases, and most importantly to establish a uniform framework to carry out a successful induced seismicity risk management plan for carbon storage projects in the future. These recommendations do not directly address any domestic or international regulations or standards. A complementary NRAP report makes recommendations for the assessment and management of environmental subsurface risks associated with unwanted fluid migration at GCS sites (Thomas et al., 2021) and should be referred to in order to address those additional GCS site risks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Recommended Practices for Managing Induced Seismicity Risk Associated with Geologic Carbon Storage

The geologic storage of carbon dioxide (CO 2 ) is one method to help reduce or eliminate atmospheric CO 2 emissions. The sequestered CO 2 is originally captured from the atmosphere or from a stationary industrial source and subsequently injected into a deep subsurface porous rock formation. To facilitate the successful deployment of field scale carbon storage projects, the U.S. Department of Energy (DOE) is developing tools and protocols for defensible, science-based frameworks to quantify and mitigate risks associated with the long-term storage of CO 2 . This protocol specifically addresses the risk of induced seismicity due to injection in a geologic carbon storage (GCS) site. This integrated and risk-based protocol is a product of the U.S. DOE Fossil Energy’s National Risk Assessment Partnership (NRAP), a multi-year collaborative research effort of Los Alamos National Laboratory (LANL), Lawrence Berkeley National Laboratory (LBNL), Lawrence Livermore National Laboratory (LLNL), National Energy Technology Laboratory (NETL), and Pacific Northwest National Laboratory (PNNL). These recommended practices describe a set of 7 steps to evaluate, manage, communicate, and mitigate the risk of induced seismicity at GCS sites. The base methodology of the recommended practices follows a framework similar to the $\textit{Protocol for Addressing Induced Seismicity Associated with Enhanced Geothermal Systems}$ (Majer et al., 2012), developed for the Geothermal Technology Office of the U.S. DOE. These recommended practices present a framework to systematically assess the induced seismicity risk and quantify the associated uncertainties. These recommendations are based on current research and are sufficiently general to allow for modification and application to a variety of different types of sites. The substance of the recommended practices contained herein includes both technical and non-technical issues, and covers all operational stages of the GCS project lifecycle. They start at the preliminary risk assessment phase, continue through site assessment and characterization, include best practice communication and seismic monitoring plan methodologies, discuss the evaluation and mitigation of seismic hazard and risk, and closes with an exploration of operational management plans, which conclude when the induced seismicity risk abates back to background level. The focus of these recommendations is on actively managing the risks associated with induced seismicity by developing an actionable risk management plan that starts at the project proposal stage and continues through site closure through an iterative assessment and improvement process. The audience of this document is expected to include all interested stakeholders (e.g., operators, project developers, regulators, and the general public) and is expressly written to be accessible to this broad range of partners. This document is intended to disseminate knowledge gained through recent advances in the science of induced seismicity hazard and risk assessments, to provide updates based on recent experience gained by similar corollary injection-induced seismicity cases, and most importantly to establish a uniform framework to carry out a successful induced seismicity risk management plan for carbon storage projects in the future. These recommendations do not directly address any domestic or international regulations or standards. A complementary NRAP report makes recommendations for the assessment and management of environmental subsurface risks associated with unwanted fluid migration at GCS sites (Thomas et al., 2021) and should be referred to in order to address those additional GCS site risks

58 GEOSCIENCES↗
DOE OSTI2020

Time-Lapse Vertical Seismic Profiling (VSP) for CO 2 Storage in a Depleted Oil Field in Northern Michigan

The goal of the Dover 33 VSP study was to test the effectiveness of time-lapse VSP for detecting and delineating a plume of more than 271,000 tonnes of CO 2 injected into the Brown Niagaran and A-1 Carbonate formations within the Dover 33 reef between March 2013 and September 2016. Five 2D walkaway VSP (WVSP) source lines were acquired by SIGMA3 in September 2016 to investigate the possible time-lapse response in both P-wave and PS-wave seismic data. The data was compared to the same survey geometry acquired in March 2013 by SR2020.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Influence of Fault Architecture on Induced Earthquake Sequence Evolution Revealed by High-Resolution Focal Mechanism Solutions

The increasing seismicity and improved seismic observation network in recent years provide an opportunity to explore factors that influence the triggering processes, spatiotemporal evolution, and maximum magnitude of induced sequences. We map the fault architecture and stress state of four induced sequences in Oklahoma to determine their influence on the seismicity. We systematically relocate the earthquakes and compute hundreds of focal mechanisms of small to medium events (1.0 < M < 5.1) using various techniques, including machine learning, for the Guthrie, Woodward, Cushing, and Fairview sequences. The detailed fault geometry and spatiotemporal evolution of seismicity and stress states reveal different dominant driving forces for each sequence. In Cushing and Fairview (largest event ≥M5.0), the main fault structures are near-vertical narrow strike-slip faults, with most of the small earthquake fault planes optimally oriented. The two sequences exhibit discontinuous temporal migration but strong earthquake self-driven rupture growth. In Guthrie and Woodward (largest event <M5.0), the two sequences show more complex diffuse fault structures with varying dipping angles along depth. The inverted focal mechanisms show a mix of strike-slip faulting and normal faulting in both sequences, and the normal faulting events are less optimally oriented than strike-slip events. The two sequences are dominated by continuous diffusive migration in time driven by pore pressure propagation. The above results suggest that fault architecture and stress state influence sequence evolution, major driving forces, and possibly maximum magnitude.

58 GEOSCIENCES↗