Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic activity”

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 253 records · Page 14

High–Pressure Deformation of Iron–Nickel–Silicon Alloys and Implications for Earth’s Inner Core

Earth’s inner core exhibits strong seismic anisotropy, often attributed to the alignment of hexagonal close-packed iron (hcp-Fe) alloy crystallites with the Earth’s poles. How this alignment developed depends on material properties of the alloy and is important to our understanding of the core’s crystallization history and active geodynamical forcing. Previous studies suggested that hcp-Fe is weak under deep Earth conditions but did not investigate the effects of the lighter elements known to be part of the inner core alloy. Here, we present results from radial X-ray diffraction experiments in a diamond anvil cell that constrain the strength and deformation properties of iron-nickel-silicon (Fe–Ni–Si) alloys up to 60 GPa. We also show the results of laser heating to 1650 K to evaluate the effect of temperature. Observed alloy textures suggest different relative activities of the various hcp deformation mechanisms compared to pure Fe, but these textures could still account for the theorized polar alignment. Fe–Ni–Si alloys are mechanically stronger than Fe and Fe–Ni; extrapolated to inner core conditions, Si-bearing alloys may be more than an order of magnitude stronger. Furthermore, this enhanced strength proportionally reduces the effectivity of dislocation creep as a deformation mechanism, which may suggest that texture developed during crystallization rather than as the result of postsolidification plastic flow.

36 MATERIALS SCIENCE↗

$S$ Hmax orientation in the Alpine region from observations of stress-induced anisotropy of nonlinear elasticity

The orientation of $S$ Hmax is commonly estimated from in situ borehole breakouts and earthquake focal mechanisms. Borehole measurements are expensive, and therefore sparse, and earthquake measurements can only be made in regions with many well-characterized earthquakes. Here, we derive the stress-field orientation using stress-induced anisotropy in nonlinear elasticity. In this method, we measure the strain derivative of velocity as a function of azimuth. We use a natural pump-probe (NPP) approach which consists of measuring elastic wave speed using empirical Green’s functions (probe) at different points of the earth tidal strain cycle (pump). The approach is validated using a larger data set in the Northern Alpine Foreland region where the orientation of maximum horizontal compressive stress is known from borehole breakouts and drilling-induced fractures. The technique resolves NNW-SSW to N-S directed $S$ Hmax which is in good agreement with conventional methods and the recent crustal stress model. We confirm that the NPP method can be applied to dense large-scale seismic arrays. The technique is then applied to the Southern Alps to understand the contemporary stress pattern associated with the ongoing deformation due to counterclockwise rotation of the Adriatic plate with respect to the European plate. Our results explain why the two major faults in Northeastern Italy, the Giudicarie Fault and the Periadriatic Line (Pustertal–Gailtal Fault) are currently inactive, while the currently acting stress field allows faults in Slovenia to deform actively. We have demonstrated that the pump-probe method has the potential to fill in the measurement gap left by conventional approaches, both in terms of regional coverage and in depth.

58 GEOSCIENCES↗

CarbonSAFE Illinois (Macon County Final Report)

This report summarizes an extensive work package to present key activities and findings towards evaluating sites for geologic storage of CO 2 . The project was undertaken as a series of tasks that addressed risk, stakeholder engagement, business and economic issues, permitting, drilling and well testing, geologic characterization, seismic surveys, modeling, CO 2 source suitability and infrastructure development. Detailed reports associated with these tasks are included as appendices to this report.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Microseismic Monitoring Study to Assess the Potential for Induced Seismicity in a Depleted Oil Field 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. Microseismic monitoring is the passive recording of very small-scale seismic energy events occurring underground. Microseismic monitoring has been proposed as a monitoring technology for CCUS sites to monitor fracturing before it becomes sufficiently extensive to cause leakage. This study aimed to determine if CO 2 injection into the Niagaran pinnacle reefs in northern Michigan is likely to generate microseismic events. The Niagaran reefs are relatively small, closed features making it possible to observe a large pressure increase by injecting a relatively small volume of CO 2 . In this study, two microseismic monitoring events were conducted 39 months apart during re-pressurization of the Dover 33 reef to evaluate the potential for CO 2 -injection induced seismicity in Silurian-age carbonate reef depleted oil reservoirs.

01 COAL, LIGNITE, AND PEAT↗

Summary of NDC Capacity Building Workshop and Regional Seismic Travel Time (RSTT) in combination with Data Sharing and Integration Training

The NNSA Seismic Cooperation Program (SCP) sponsored Stephen Myers (LLNL), Michael Begnaud (LANL), Brian Young (SNL) and Istvan Bondar (Research Center for Astronomy and Earth Sciences, Hungary) to serve as a presenters/trainers at the “NDC Capacity Building Workshop and Regional Seismic Travel Time (RSTT) in combination with Data Sharing and Integration Training” September 4-8 2022 in Muscat, Oman (See Appendix A for the agenda). The workshop and training (workshop from here forward) was organized by the Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO) Provisional Technical Secretariat (PTS). The first half of the week was devoted to NDC workshop activities, and the second half was devoted to RSTT training. Fifty-five participants from 27 countries and the CTBTO-PTS attended the 5-day workshop (See Appendix B for list of participants and countries of origin). Presentations from the PTS described the International Monitoring System (IMS), International Data Centre (IDC) products, and metrics of regional data utilization. Contributed presentations from each country’s scientists included descriptions of regional and national networks, methods of data analysis, and needs for material and technical assistance. Training included an overview of the RSTT method and instruction on how to locate seismic events with the iLoc program, which utilizes RSTT travel times to reduce bias in event location estimates. Methods of seismic tomography and the need for a high-quality tomographic set, including seismological “ground truth”, were emphasized. Seismological “ground truth” or “GT” is a term that has come to mean both events with known location and events with well-characterized locations that are estimated using seismological data, typically with epicenter accuracy of 5 km or better. Notably, the instructional platform has migrated from UNIX shell scripts to Jupyter Notebooks. Jupyter Notebooks have the advantage being more visually intuitive, including display of graphics within the notebook. Each notebook includes every processing step that participants need to reproduce the entire exercise.

58 GEOSCIENCES↗

Simulation of Physics-Based 0-10Hz Strong Motion Using High Performance Computing Supporting Refinements to Regional Ground Motion Models for the Central Eastern US

In collaboration with the U.S. Nuclear Regulatory Commission (NRC) the LLNL has developed a computationally efficient simulation platform designed to perform physics-based ground motion simulations for crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using high-performance computing. The main objective of the earthquake simulations was to use synthetic ground motion to provide constrains to refinements of existing ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-10Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model (1994). In simulations performed during the first stage of this project the Saikia’s velocity model demonstrated better performance in modelling high frequency regional wave propagation for the CEUS region recorded during the Mw5.0 November 7, 2016, Cushing Oklahoma (Taylor et al., 2017), and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes. The proposed regional 3D model includes random perturbations to the 1D background model using the stochastic scheme of Pitarka and Mellors (2021). In addition, validation analysis of the rupture generator and regional wave propagation models, using comparisons with different GMMs for Mw6.5 and Mw7.0 scenario earthquakes in the CEUS region resulted in a very good match between the simulated and empirical ground motion models. For the purposes of seismic hazard assessment at the existing and planned nuclear power plants, NRC is interested in studies aimed at improving the current ground motion models (GMM) for both Stable Continental Regions (SCR) in the Central and Eastern US and Active Crustal Regions (ACR) in the Western US. Due to lack of recorded data, these improvements require synthetic data for short fault distances and large magnitude earthquakes for which the existing recorded data is not enough to uniquely constrain the GMMs. The need for simulations and strong motion data is especially critical for the CEUS region where we do not have recorded data from potentially large damaging earthquakes with moment magnitudes 6.0 and higher. In this the project, we focused on 10Hz simulations of Mw7.0 scenario earthquakes with strike slip and thrust faulting mechanisms. We used more than 50 Mw7.0 earthquake rupture scenarios to investigate the ground motion uncertainty due to unknown earthquake rupture parameters, in particular, the slip distribution, rupture velocity, and faulting mechanism, and their implication on ground motion amplification due to forward rupture directivity effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Subtask 1.5 – CO2 Injection Monitoring with an Optimized Scalable, Automated, Semipermanent Seismic Array

The scalable, automated, semipermanent seismic array (SASSA) method is a flexible and relatively cost-effective surface geophysical method for regular time-lapse monitoring of the movement of injected carbon dioxide (CO2) in a reservoir for CO2 enhanced oil recovery (EOR) or geologic CO2 storage operations. It has the advantages of a low-environmental-footprint while monitoring regions of a reservoir from the surface without the need for a regular grid distribution of receivers. Automated data collection is possible. As only time-lapse amplitude changes at the reservoir level due to CO2 movement within the reservoir are monitored, the turnaround time to deliver results from the SASSA method can be short, without the need for long, time-consuming data-processing workflows. As data is collected and processed, incremental information can be provided to the field operator. The Energy & Environmental Research Center (EERC) conducted a SASSA field test from September 2018 to November 2020 in a portion of the Bell Creek Field in Montana, which implemented new CO2 EOR field activities during the study period. Lessons learned from a proof-of-concept study were incorporated to improve the data quality of the SASSA method and demonstrate the viability of the technology. The EERC implemented several enhancements to improve data quality, including 1) an iterative survey design, which allowed placing the receivers in strategic locations where the movement of the CO2 in the reservoir could be tracked with minimum interference by the cultural noise in the study area; 2) the use of powerful seismic sources in the form of surface orbital vibrators, and 3) data acquisition during optimal periods. History-matched reservoir simulation was performed to predict gas saturation and pressure response induced by CO2 injection in the study area. The results were compared with the SASSA-measured responses to CO2 injection as a partial validation technique. A match between the two methods was observed for most of the SASSA points predicted to have intersected a CO2 saturation change. The validated results provide confidence that the SASSA method can be used independently as a CO2 saturation monitoring technique. As data are collected and processed, incremental information can be provided to the field operator. The critical components of the SASSA workflow for a successful application of the method are the following: Iterative survey design with information about CO2 injection activities from the oilfield operator. A detailed CO2 injection plan is the key driver to select the strategic monitoring location of the SASSA sensors. After this information is incorporated in the initial distribution of sources and receivers in the study area, high-resolution satellite images are used to identify ground locations not affected by cultural noise sources, such as power lines, pipelines/flow lines, or roadways. In the next iteration of the survey design, a scouting trip to the study area is needed to understand more details of the noise sources identified in the previous step and the intensity of the field activities that can also generate noise during the monitoring. Integrating the information from the scouting trip into the survey design to select the optimum source and receiver locations is the final step. Noise attenuation. The variety of noise types during seismic monitoring of an oil field is enormous. Tailored noise characterization and processing at a node-by-node level can enhance the performance and sensitivity of the SASSA technique. Future advancements that could improve the efficiency and application of the SASSA technology include: Gaining a better understanding of the noise field produced by the seismic source to aid the choice of receiver location. Surface noise from the source can overwhelm the small signal changes due to CO2 that the SASSA method measures. Improved data-processing workflow to automatically analyze and adapt to dynamic noise conditions associated with industrial settings. This subtask was funded through the EERC–DOE Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE- FE0024233.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Cross-Well Seismic Monitoring of CO 2 Injected Into the A-1 Carbonate and Brown Niagaran Formations at the Chester 16 Reef

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. A cross-well seismic survey was acquired in the Chester 16 reef to attempt to locate 85,000 tonnes of carbon dioxide (CO 2 ) that were injected into the A-1 Carbonate and Brown Niagaran Formations earlier. The technique of conducting multiple cross-well seismic surveys over time, including a pre-CO 2 injection (baseline) survey, has been used elsewhere to monitor CO 2 injected into the subsurface. In this study, a baseline cross-well survey was not obtained; nevertheless, it was possible to generate an image that is a plausible representation of the CO 2 plume. This conclusion is supported by other monitoring and modeling results from the Chester 16 reef that provide an independent indication about the likely position of the injected CO 2 .

01 COAL, LIGNITE, AND PEAT↗

Building an EPA Class VI Permit Application

Summary To accelerate the commercialization of carbon capture and storage (CCS), the US Department of Energy (US DOE) is building on decades of characterization efforts and pilot-scale projects through their CarbonSAFE program. Administered through their National Energy Technology Laboratory, this program seeks to bring fully integrated projects to the sector that can store more than 50 million tonnes of CO2 over a 30-year period. The program, which was enacted before the enhancement of Internal Revenue Code Section 45Q, is in the capture assessment, characterization, and permitting phase. The objectives of this paper are to discuss (a) the injection permitting requirements of the CarbonSAFE projects; (b) information gathering in support of the permit; (c) the timelines of field development and permit-related activities; (d) the major technical components of the field development plan; and (e) early feedback from the regulators toward acceptance of the permit. In Mississippi, more than 30,000 acres have been characterized by six deep characterization wells, a deep groundwater well, and 92 line miles of 2D seismic as part of the CarbonSAFE Project ECO2S. During the acquisition of seismic data, all receiver lines were live, which resulted in the generation of a pseudo-3D seismic design. The incorporation of a 3D seismic survey was not included as part of this project due to logistical difficulties presented by the undulating, wooded surface terrain. A suite of openhole geophysical logs was taken from each well, allowing for a detailed interpretation of prospective storage reservoirs and confining intervals to complement the analysis carried out on the 290 ft of a whole core that was cut through the prospective confining zone and storage reservoir. The detailed geologic and reservoir data were assembled and entered into a 3D model to assess the injection capacity and the area of review (AoR). This information fed into the detailed corrective action, monitoring, testing, and postinjection site care (PISC) modeling. The results have been exceptional. The geologic assessment has revealed three primary storage targets, ranging in depth from 3,500 ft to 6,000 ft. These storage reservoirs net 1,300 ft of sandstone, with mean porosity and permeability of 29% and 3.6 darcies, respectively. Together, these reservoirs have storage capacities that may exceed 20 million tonnes per square mile, making this a gigatonne prospect. Forward modeling of the project resulted in an AoR of 16 sq miles, injecting about 8000 t/d, for 30 years, via two deep injection wells. The excellent confining characteristics of the caprock, relatively simple geologic structure, and lack of historical well drilling activity in this area provide excellent containment of the injected CO2. Based on this work, the project has proposed 20 years of PISC. To date, only two US CO2 injection permits have been granted. These projects relied on a singular capture point feeding a singular sequestration point (source to sink), and considerations have not been made to garner CO2 emissions from other industrial sources. The Kemper County Storage Complex is a first-of-its-kind storage hub concept that looks to develop an area capable of storing significant quantities of CO2 from the region. Also, this work will show how characterization efforts, geological and numerical modeling efforts, and plan development were constructed in support of permit and incentives acceptance.

Energy & Fuels↗

Time-Lapse Vertical Seismic Profiling (VSP) for CO 2 Storage in a Depleted Oil Field 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. This report describes the Dover 33 VSP study to test the effectiveness of time-lapse Vertical Seismic Profile (VSP) for detecting and delineating a plume of more than 271,000 tonnes of CO2 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. The seismic energy was recorded into an 80-level, three-component geophone array deployed into the Dover 1-33 well and placed just above the reef. P-wave and PS-wave reflection images were produced for each of the three source lines from both the 2013 baseline and 2016 monitor surveys. The images were then compared to look for changes in the reflectivity at and around the injection location that might indicate how the CO 2 has moved over this time period. This technique revealed several localized areas with sizable impedance differences inside the reef where CO 2 would be expected; however, a large number of similar impedance “hotspots” were also detected outside the reef in areas where injected CO 2 would not be expected. Therefore, a second analysis that involved calculating P-wave and S-wave travel time differences between the 2013 and 2016 VSPs was conducted to look for a change that could be caused by the CO 2 plume.

01 COAL, LIGNITE, AND PEAT↗

Distributed Acoustic Sensing (DAS) Seismic Monitoring of CO 2 Injected for Enhanced Oil Recovery in Northern Michigan (Volume II.C)

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. This report describes the monitoring study conducted to assess the effectiveness of DAS (Distributed Acoustic Sensing) -based VSP (Vertical Seismic Profiling) technology for delineating CO 2 injected into the Silurian-age pinnacle reefs in northern Michigan, the host rocks for the MRCSP Phase III demonstration project. The DAS VSP study was conducted in the Chester 16 reef, one of several reefs in Otsego County Michigan that is operated by Core Energy, LLC of Traverse City, Michigan. Time-lapse DAS VSP was implemented at the Chester 16 reef to attempt to detect approximately 85,000 tonnes of CO 2 injected into the A-1 Carbonate and Brown Niagaran Formations. A baseline survey was conducted in February 2017 prior to injecting CO 2 and a repeat survey was conducted in August 2018. During the interim period between the baseline and repeat surveys, CO 2 was injected into the Chester 16 reef via the 6-16 injection well without production (withdrawal) of fluids from the reef. A grid of 181 source positions consisting of 44 vibrator positions, plus 137 dynamite shot locations, was used to give approximately continuous spatial coverage of the injection zone (A-1 Carbonate and upper Brown Niagaran) in the area between the two wells.

01 COAL, LIGNITE, AND PEAT↗

Seismoacoustic Monitoring of Pyroprocessing Equipment

With the projected increase in world nuclear capacity comes the hurdle of spent nuclear fuel. Pyroprocessing is one method to process irradiated fuel by using high temperatures and electrochemical steps to separate radioactive components. High temperatures, high radiation levels, and equipment confined to a heavily shielded hot cell are a few of the challenges introduced regarding safeguards for pyroprocessing. The first stage of the process utilizes an element chopper to cut irradiated fuel into smaller pieces, offering the potential to identify equipment operation through seismoacoustic monitoring. Eleven seismic and infrasound sensors were deployed at various locations and distances near a pyroprocessing hot cell to evaluate signals emitted from the chopper. Signals were processed using short-time Fourier transforms and manually scanned to locate chopping events. Operator logs were then collected for a ground truth comparison. Results indicate that both infrasound and seismic signals offer the ability to accurately measure this stage of pyroprocessing, though seismic signals are more prominent. Sensor locations are also evaluated to determine where the signal is detectable. This type of monitoring offers a method to aid in safeguards and proliferation detection by validating a pyroprocessing facility’s schedule and identifying activity that does not align with records.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Controlling Factors Analyses of Sediment Dispersal From Source-to-Sink Systems of Lacustrine Rift Basin: A Case Study From Paleogene Shahejie-Dongying Formations of Nanpu Sag, Eastern China

Tectonic activity not only shapes the basic stratigraphic framework of rift basins, but also profoundly affects the sediment dispersal in rift basins. In this study, analyses of heavy mineral assemblages in different periods demonstrate that there are three obvious tectono-sedimentary evolutionary stages (E s 3 –E s 2 , E s 1 –E d 2 , and E d 1 , respectively) in the Paleogene provenance area of Nanpu Sag, and the volume of sand bodies increases from the bottom of the Paleogene Shahejie (E s ) Formation to the top of the Dongying (E d ) Formation in Nanpu Sag. Besides, this study comprehensively utilize the analyses of seismic interpretation, palynology, heavy mineral, and borehole core samples to investigate the controlling factors of sediment dispersal in the rift basin. The assemblages of heavy minerals in different periods reflect the rock composition and catchment area of different provenance areas, and their vertical differences reflect the evolution process of the provenance area and reflect the uplift-denudation process of the provenance area. The results reveal that the synergy of the evolution of tectonic activity and the adjustment of topographical evolution are the main controlling factors of sediment dispersal in Nanpu Sag, while climate change is not the main controlling factor. We conclude that an increased sediment supply rate in the long term reflects the control of tectonic activity on basin topography, rather than climate fluctuations. The differences in morphological modification result in differential sediment dispersal, which is principally related to the differential extrusion of the fault system. The catchment area and provenance distance adjustment is evidenced by the vertical changes of heavy mineral characteristics of single-well and interaction and linkage of boundary faults, and the adjustment of topography evolution. A consideration is that the interaction and linkage of boundary faults and complex subsidence history are multi-directional, and differential evolution of provenance area is universal in lacustrine rift basins, all of this highlights the adjustment of sediment pathways generated by this characteristic of rift basins and emphasizes the importance of controlling factors analyses in understanding differential sediment dispersal that presents in the rift basins. Besides, four sets of sediment dispersal patterns were delineated based on different developmental regions in the rift basin, which are fault segmental point and multi-stage fault terrace, single-stage fault terrace and axial fault valley, axial fault terrace, and paleo-terrace and axial fault valley, respectively. This study has a certain guiding significance for the prediction of the spatial distribution of sand bodies in the rift basin and the exploration of potential oil and gas targets in the rift basin.

Yu, Jifeng↗

An autoencoder compression approach for accelerating large-scale inverse problems

Abstract Partial differential equation (PDE)-constrained inverse problems are some of the most challenging and computationally demanding problems in computational science today. Fine meshes required to accurately compute the PDE solution introduce an enormous number of parameters and require large-scale computing resources such as more processors and more memory to solve such systems in a reasonable time. For inverse problems constrained by time-dependent PDEs, the adjoint method often employed to compute gradients and higher order derivatives efficiently requires solving a time-reversed, so-called adjoint PDE that depends on the forward PDE solution at each timestep. This necessitates the storage of a high-dimensional forward solution vector at every timestep. Such a procedure quickly exhausts the available memory resources. Several approaches that trade additional computation for reduced memory footprint have been proposed to mitigate the memory bottleneck, including checkpointing and compression strategies. In this work, we propose a close-to-ideal scalable compression approach using autoencoders to eliminate the need for checkpointing and substantial memory storage, thereby reducing the time-to-solution and memory requirements. We compare our approach with checkpointing and an off-the-shelf compression approach on an earth-scale ill-posed seismic inverse problem. The results verify the expected close-to-ideal speedup for the gradient and Hessian-vector product using the proposed autoencoder compression approach. To highlight the usefulness of the proposed approach, we combine the autoencoder compression with the data-informed active subspace (DIAS) prior showing how the DIAS method can be affordably extended to large-scale problems without the need for checkpointing and large memory.

Mathematics↗

North Dakota Integrated Carbon Storage Complex Feasibility Study. Final report

In spring 2017, the Energy & Environmental Research Center (EERC) initiated an effort to determine the feasibility of developing a commercial-scale CO 2 geologic storage complex able to store 50+ million tonnes (Mt) of CO 2 in central North Dakota safely, permanently, and economically. The objective was to fulfill the goals of the U.S. Department of Energy (DOE) Carbon Storage Assurance Facility Enterprise (CarbonSAFE) Initiative and address technical and nontechnical challenges specific to commercial-scale deployment of a CO2 storage project. The findings clearly show that the concept of capturing CO 2 from a lignite-fired electrical generation facility in central North Dakota and safely and permanently storing the CO 2 in the deep subsurface is indeed technically, economically, and socially feasible. This project evaluated two study areas and their respective geologic storage complexes located adjacent to separate coal-fired facilities in North Dakota: the Basin Electric Power Cooperative (BEPC)-owned Great Plains Synfuels Plant (GPSP) and the Minnkota Power Cooperative (Minnkota)-owned Milton R. Young Station (MRYS). These locations, one with CO 2 capture in place and an existing CO 2 pipeline, are bolstered by progressive North Dakota pore space ownership and long-term liability laws. These elements and a motivated team created an ideal synergistic scenario for ensuring success of the CarbonSAFE Initiative and promoting North Dakota’s statewide vision for carbon management. The project included drilling two new geologic characterization wells, integrating an existing 3-D seismic survey, creating a geologic model subsequently used for injection simulation, a risk assessment, public outreach, and generating a site development plan based on results. In addition, the performance of select National Risk Assessment Partnership tools was evaluated. The geologic characterization wells were drilled ~5600 feet deep to the Broom Creek Formation; ~350 feet of core was retrieved from each well. The core included the Broom Creek (targeted injection zone) and a portion of the overlying Opeche Shale (seal). The Flemmer-1 well, west of Beulah, North Dakota, yielded 169 feet of sandstone. The BNI-1 well located south of Center, North Dakota, yielded 124 feet of sandstone. In each case, laboratory analysis of the sandstone showed permeability in the 300–1000-mD range, with porosity of 20%–30%. The Flemmer-1 well was sited within the boundaries of an existing 3-D seismic survey. Colocating the well with the seismic survey maximized the relationship between new and legacy data and developed a first-of-its-kind interpretation of the geologic fabric of the Broom Creek. Geologic characterization data were integrated into a 5544-mi 2 geocellular model that encompassed both new wells and stratigraphy from the surface to the Amsden Formation (underlying the Broom Creek). The model was later expanded vertically to include the deeper Black Island–Deadwood interval to examine its potential viability as a storage target. The geocellular model provided the foundation for dynamic simulation of CO 2 into the Broom Creek. Results of the simulation suggest that the Broom Creek could accept the DOE target rate of 2 Mt/yr of CO 2 into as few as two wells. To bracket the expected capture from MRYS, simulations were also investigated for a 4-Mt/yr rate near MRYS. Although more wells are needed (two additional), the Broom Creek still has the storage resource to accept the CO 2 at the increased rate. A risk assessment exercise was conducted to identify and assess technical and nontechnical risks that could prevent potential candidate storage complexes within the study area from serving as commercial storage sites. The assessment identified and evaluated six technical risk categories: 1) CO 2 injectivity, 2) storage capacity, 3) lateral migration of CO 2 , 4) lateral pressure propagation, 5) vertical migration of CO 2 or formation brine, and 6) induced seismicity. Following two rounds of analysis and scoring, no risks were determined to preclude continued efforts to develop carbon capture, utilization, and storage (CCUS) in central North Dakota. The risk assessment results will be used to guide future site characterization, modeling and simulation, and monitoring activities. A specific nontechnical strategic risk based on challenges that may be realized in amalgamation of pore space resulted in vertically expanding the geologic model to incorporate the potential for stacked storage in multiple saline reservoirs. By including the Black Island–Deadwood interval (the basal sedimentary reservoir in this region), the amalgamated areal extent could be reduced by as much as 45%. Working with a smaller geographic area reduces risks and costs associated with monitoring and pore space leasing. An economic evaluation incorporating capture; transport (<5 mi); Class VI wells; permitting; and monitoring, verification, and accounting suggests implementing commercial-scale CCUS is economically attractive if the federal tax benefits of 45Q are included. This is validated by Minnkota’s continued pursuit of CO 2 capture and geologic storage at MRYS through its Project Tundra initiative, indicating that there is a business case for CCUS in central North Dakota. Currently, North Dakota is the only state with underground injection control (UIC) Class VI primacy. Built into the North Dakota Century Code is a series of regulatory requirements that guide the process to obtain a Class VI CO 2 storage facility permit. As part of this project, a site development plan was compiled to assure compliance with North Dakota’s requirements to permit a commercial-scale CO 2 storage operation and includes a prospective time line encompassing a general breakdown of activities. In total, an estimated 30 months would be needed to execute the necessary steps to attain a North Dakota CO 2 storage facility permit. Outreach was an integral part of the project and encompassed any project-related activity that had contact or exposure beyond the project team. The goals of outreach were to foster an environment from which stakeholders could make informed decisions about the project and gauge community receptiveness to a CCUS project. A consistent set of messages and outreach products were developed in conjunction with an outreach advisory board that integrated project partners and team members. To gauge public acceptability of geologically storing CO 2 , 5611 households in the project area were invited to participate in an online survey. The survey results indicate that the public attitude regarding CCUS is neutral to positive, with strong sentiment that CO 2 capture and storage may be an approach to maintain the economic vitality of the region. To achieve project objectives, critical support in the form of financial backing, engineering evaluations, site access, outreach collaboration, operations data, risk assessment/evaluation, and software access was provided by BEPC, the North Dakota Industrial Commission Lignite Research Council, ALLETE Clean Energy, BNI Energy, North American Coal Corporation, Minnkota, Prairie Public Broadcasting, Computer Modelling Group Ltd., and Schlumberger.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Phase-weighted slant stacking for surface wave dispersion measurement

SUMMARY Surface wave retrieval from ambient noise records using seismic interferometry techniques has been widely used for multiscale shear wave velocity (Vs) imaging. One key step during Vs imaging is the generation of dispersion spectra and the extraction of a reliable dispersion curve from the retrieved surface waves. However, the sparse array geometry usually affects the ability for high-frequency (>1 Hz) seismic signals’ acquisition. Dispersion measurements are degraded by array response due to sparse sampling and often present smeared dispersion spectra with sidelobe artefacts. Previous studies usually focus on interferograms’ domain (e.g. cross-correlation function) and attempt to enhance coherent signals before dispersion measurement. We propose an alternative technique to explicitly deblur dispersion spectra through use of a phase-weighted slant-stacking algorithm. Numerical examples demonstrate the strength of the proposed technique to attenuate array responses as well as incoherent noise. Three different field examples prove the flexibility and superiority of the proposed technique: the first data set consists of ambient noise records acquired using a nodal seismometer array; the second data set utilizes distributed acoustic sensing (DAS) and a marine fibre-optic cable to acquire a similar ambient noise data set; the last data set is a vibrator-based active-source surface wave data. The enhanced dispersion measurements provide cleaner and higher-resolution spectra without distortions which will assist both human interpreters as well as ML algorithms in efficiently picking curves for subsequent Vs inversion.

58 GEOSCIENCES↗

Seismotectonic evidence for subduction beneath the Eastern Greater Caucasus

SUMMARY We generated high-resolution 3-D seismic velocity models as well as a relocated earthquake catalogue across the eastern Greater Caucasus and Kura basins. This work was done using data from the recently upgraded Republic Seismological Survey Center's (RSSC) seismic network. We generated our tomographic images of crustal velocity structure in Azerbaijan using double-difference inversions (i.e. tomoDD and hypoDD). Earthquake catalogues from the RSSC between 2011 and 2016 were used; these catalogues include absolute arrival times of 103 288 P- and 120 952 S-wave traveltime picks for 7574 events recorded at 35 stations in Azerbaijan. Beginning with a layered, 1-D velocity model that was estimated using VELEST, we inverted simultaneously for relative location, Vp and Vs on a 3-D grid with dimensions 670 × 445 × 45 km, with a uniform grid spacing of 55 × 55 × 5 km for all of Azerbaijan. We observe that the relocated hypocentres cluster into two depth ranges, at the surface and at depth, that appear to correspond to major fault zones and the top of a subducting plate. Additionally, we note intermediate depth seismicity (∼50–60 km) beneath the Kura Basin, and a northward deepening of earthquake depths. Seismic velocities vary significantly throughout the study region; we observe very slow velocities throughout the Kura Basin between 5 and 15 km, and elevated velocities at 20–35 km. The wholesale velocity structure and seismic structure of Kura Basin strongly mirrors that of the Caspian Sea, which suggests that the geodynamics of the Caspian continue westwards into Azerbaijan. The key results of this study suggest that the northward subduction observed in the Caspian Sea continues beneath the Eastern Greater Caucasus, as well as provides evidence for active faulting along the southern margin of the mountain range.

Gunnels, Michael↗

Local Explosion Detection and Infrasound Localization by Reverse Time Migration Using 3-D Finite-Difference Wave Propagation

Infrasound data are routinely used to detect and locate volcanic and other explosions, using both arrays and single sensor networks. However, at local distances (< km) topography often complicates acoustic propagation, resulting in inaccurate acoustic travel times leading to biased source locations when assuming straight-line propagation. Here we present a new method, termed Reverse Time Migration-Finite-Difference Time Domain (RTM-FDTD), that integrates numerical modeling into the standard RTM back-projection process. Travel time information is computed across the entire potential source grid via FDTD modeling to incorporate the effects of topography. The waveforms are then back-projected and stacked at each grid point, with the stack maximum corresponding to the likely source. We apply our method to three volcanoes with different network configurations, source-receiver distances, and topography. At Yasur Volcano, Vanuatu, RTM-FDTD locates explosions within ~20 m of the source and differentiates between multiple vents. RTM-FDTD produces a more accurate location for the two Yasur subcraters than standard RTM and doubles the number of detected events. At Sakurajima Volcano, Japan, RTM-FDTD locates the source within 50 m of the active vent despite notable topographic blocking. The RTM-FDTD location is similar to that from the Time Reversal Mirror method, but is more computationally efficient. Lastly, at Shishaldin Volcano, Alaska, RTM and RTM-FDTD both produce realistic source locations (<50 m) for ground-coupled airwaves recorded on a four-station seismic network. We show that RTM is an effective method to detect and locate infrasonic sources across a variety of scenarios, and by integrating numerical modeling, RTM-FDTD produces more accurate source locations and increases the detection capability.

58 GEOSCIENCES↗