SEARCH · Engineering Papers
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.
Presentation at the 17th Greenhouse Gas Technologies Conference in Calgary, Alberta, Canada - describing research progress and computational tools developed for geologic carbon storage risk assessment under the National Risk Assessment Partnership.
Effective seismic monitoring is essential for verifying CO₂ containment, detecting potential leakage, and optimizing operational decisions in geologic carbon storage. Here, this study presents a time-adaptive, elastic-wave sensitivity-guided framework for designing cost-effective seismic monitoring layouts for tracking CO₂ plume migration. The method is based on elastic-wave sensitivity analysis, which quantifies how variations in subsurface properties impact seismic wavefields. Two complementary design strategies are developed: one based on selecting a fixed number of seismic sources (Method A), and the other based on selecting source–receiver pairs contributing to a fixed fraction of cumulative elastic-wave sensitivity energy (Method B). The optimization workflow to identify source–receiver configurations with the highest detection potential is demonstrated using a hypothetical GCS scenario at the Kimberlina site in California using simulations of elastic-wave sensitivity data at multiple post-injection timesteps. Results show that both strategies adapt to evolving plume geometries and wavefield sensitivities, with Method B offering broader spatial coverage and Method A ensuring simpler deployment. This framework enables site-specific, cost-effective, and risk-informed seismic survey designs, enhancing the ability to monitor CO₂ migration over time in evolving geological environments
Data collected via seismic and infrasound array deployments are leveraged in the geosciences to detect and characterize a myriad of natural and anthropogenic sources. These deployments consist of numerous sensors placed in a predetermined configuration to amplify signal strength and improve the efficacy of array processing techniques used to measure signal directionality and waveform coherence. High‐fidelity feature extraction is often predicated on interstation distance as well as the frequency content and wavelength of an incident signal. Numerous array processing softwares analyze data in sequential frequency bands to obtain a more detailed characterization of a signal. However, current algorithms are limited in their ability to determine optimal array configuration for each band. We introduce an open‐source Python code, called Cardinal, to process seismic and infrasound array data in discretized time–frequency space with the option of applying an adaptive array design to determine optimal subarray configuration for each frequency band. To reduce computational time, the array processing step can be run in parallel using multithreading. Furthermore, the software has the capability to aggregate array processing results from different time–frequency pixels to produce separate sets of detections, or families, with added utility via the application of an adaptive semblance threshold, which aids in isolating signals‐of‐interest from coherent background noise. Upon appropriate configuration, Cardinal exhibits the potential to combine distinct seismic and infrasound phases into separate families.
Ambitious international greenhouse gas emissions reduction targets demand a rapid transformation to a low-carbon economy. This transformation includes the accelerated adoption of carbon dioxide (CO2) capture and storage (CCS) technology. However, as with any large-scale engineering enterprise, the widespread commercial-scale deployment of geologic carbon storage (GCS) raises important questions about technology and cost-effectiveness, safety, environmental risk, and long-term liability. Effectively assessing and managing risks and liability associated with GCS projects is a key technical need throughout the project life cycle-from site selection and permitting to monitoring design, operational risk management, and post-operational site closure. This presentation highlights recent advancements in tools for quantitative risk assessment, being developed by the National Risk Assessment Partnership (NRAP). NRAP is a multi-year, multinational laboratory research collaboration sponsored by the U.S. Department of Energy's Office of Fossil Energy and Carbon Management. Our focus will be on these tools' applications in addressing critical stakeholder questions related to supporting permitting to ensure secure and environmentally protective storage; designing effective and efficient monitoring plans; evaluating the effectiveness of remedial actions and risk management alternatives; and informing liability assessment and investment decisions. This paper will detail the key functionality of NRAP’s Open-Source Integrated Assessment Model (NRAP-Open-IAM), a computational framework for assessing leakage risk and containment assurance. This model features streamlined workflows for calculating leakage risk profiles, delineating risk-based area of review, and assessing contingency plans and post-injection site care requirements. ORION is an open-source, observation-based ensemble forecasting toolkit to help operators assess the seismic hazard at a carbon storage site. The State of Stress Analysis Tool (SOSAT), designed to assess subsurface stress conditions and evaluate geomechanical risk resulting from CO2 injection in an area of interest will also be presented. We will also introduce a prototype model to evaluate storage project costs and liability associated with risk management. The Technoeconomic and Liability Evaluation for Storage (TALES) model uses results from forecasts of leakage and induced seismicity risk to estimate the lifecycle cost of managing risk. Finally, a preliminary example of how the NRAP Risk-based Adaptive Monitoring Plan (RAMP) tool can be used to design efficient and effective site monitoring plans and estimate the detectability of fluid leakage will be provided. The relevance of these tools for addressing key stakeholder questions amidst uncertainty will be emphasized.
PySolate is a Python‐based toolset that implements the continuous wavelet transform and nonlinear thresholding operations to denoise or designal seismic data, following Langston and Mousavi (2019). This filtering approach can remove microseismic noise to isolate intermediate‐period seismic signals that are key to enabling full‐waveform modeling and analysis of smaller‐magnitude regional events. This approach is best for the application to signals with frequency or time separation of signal and noise, in contrast to Fourier analysis, which is effective when signal and noise are separated in frequency. We demonstrate the Python toolset using the six announced Democratic People’s Republic of Korea declared nuclear tests, showing the effectiveness of isolating the seismic signal compared to standard bandpass filtering. In conclusion, we also demonstrate the ease of using the toolset with any Python processing tools.
Site characterization for underground injection and storage of gigatonne-scale CO₂ requires reliable and cost-effective methods to detect and characterize faults and fractures and to assess their stress state and fault activation potential. This is critical, as wastewater injection and disposal have been shown to activate faults and induce earthquakes, and CO₂ leakage remains a key concern for long-term storage. In this project, we developed seismic methods to detect and characterize large-scale sedimentary and crystalline basement faults and associated small-scale fractures below conventional seismic imaging resolution using multicomponent (9C) surface seismic data. Machine learning was used to automatically interpret large-scale faults, providing key information for estimating the maximum magnitude of potential induced earthquakes. High-fidelity imaging was achieved by exploiting redundancy across multiple elastic wave modes, where independent images from different modes and frequencies cross-validate each other. We also used our nonlinear signal comparison (NLSC) method for ground roll removal, improving data quality in complex near-surface conditions. The methods were validated using field data acquired in central Montana. Results show that basement faults extend into the sedimentary section and that small-scale fractures are widespread above the basement. The inferred stress orientation is consistent with regional stress data, and the estimated maximum induced earthquake magnitude is small (Mw ~2.3). The developed workflow provides a practical approach for fault and fracture characterization and for assessing induced seismicity and leakage risk. It is directly applicable to CO₂ storage site selection and to other subsurface systems.
The author has identified the following significant results. Four previously unknown faults were discovered in basement terrane of the Peninsular Ranges. These have been named the San Ysidro Creek fault, Thing Valley fault, Canyon City fault, and Warren Canyon fault. In addition fault gouge and breccia were recognized along the San Diego River fault. Study of features on Skylab imagery and review of geologic and seismic data suggest that the risk of a damaging earthquake is greater along the northwestern portion of the Elsinore fault than along the southeastern portion. Physiographic indicators of active faulting along the Garlock fault identifiable in Skylab imagery include scarps, linear ridges, shutter ridges, faceted ridges, linear valleys, undrained depressions and offset drainage. The following previously unrecognized fault segments are postulated for the Salton Trough Area: (1) An extension of a previously known fault in the San Andreas fault set located southeast of the Salton Sea; (2) An extension of the active San Jacinto fault zone along a tonal change in cultivated fields across Mexicali Valley ( the tonal change may represent different soil conditions along opposite sides of a fault). For the Skylab and LANDSAT images studied, pseudocolor transformations offer no advantages over the original images in the recognition of faults in Skylab and LANDSAT images. Alluvial deposits of different ages, a marble unit and iron oxide gossans of the Mojave Mining District are more readily differentiated on images prepared from ratios of individual bands of the S-192 multispectral scanner data. The San Andreas fault was also made more distinct in the 8/2 and 9/2 band ratios by enhancement of vegetation differences on opposite sides of the fault. Preliminary analysis indicates a significant earth resources potential for the discrimination of soil and rock types, including mineral alteration zones. This application should be actively pursued.
The passive seismic network deployed on the Moon during the Apollo missions operated for eight years and allowed the observation of seismic activity. More than 12500 seismic events were registered, where 28 were classified as shallow moonquakes with moment magnitudes up to 4.1. Seismic events of this nature pose a significant risk to future long-term lunar habitats; thus, these events must be carefully studied and considered in the seismic design of these structures. This paper proposes a preliminary seismic hazard assessment imposed by shallow moonquakes. The hazard assessment is performed using the Probabilistic Seismic Hazard Analysis (PSHA) methodology, considering previous studies and theories regarding the seismic environment of the Moon. The study zone covers ~860 km2 of the Taurus-Littrow Valley, containing the Apollo 17 landing site and the Lee-Lincoln lobate scarp as the considered seismic source. The seismic hazard is quantified in terms of peak ground acceleration (PGA) and spectral acceleration (5% damped pseudo-acceleration, PSA). Seismic hazard deaggregation scenarios, Uniform Hazard Spectra (UHS) for different hazard levels, and a Conditional Mean Spectrum (CMS) for a target period of 0.2 s are obtained to quantify the seismic hazard on a specific site on the Moon. The developed seismic hazard assessment provides a preliminary approach for realistic scenarios to conduct structural designs that ensure the seismic performance of fully operational long-term lunar structures.
Presentation slides on “Fracture Network Quantification for Enhanced Reservoir Characterization in CO2 Storage Sites” for the CCUS 2025 conference held in Houston, Texas March 3-5, 2025. Accurate reservoir characterization is fundamental for the safe and efficient operation of carbon capture, utilization, and storage (CCUS) projects. This process is not only critical during initial site selection but also plays a key role in continuously updating reservoir geomodels throughout injection operations to effectively manage subsurface storage sites. Our primary objective is to develop a comprehensive understanding of fracture networks within CO2 injection reservoirs, which are critical for predicting the behavior of injected fluids, ensuring storage integrity post-injection, and mitigating risks such as leakage or induced seismicity.
Carbon Capture and Storage (CCS) is in transition from first-of-a kind projects and research-orientated pilots to commercially-motivated applications. Monitoring results from many newly developed and planned large scale commercial projects are limited; however, it is worthwhile to assess their evolution and consider new strategies as part of an effort to assess and document best practices. Commercial monitoring is targeted to activities that comply with regulatory drivers and de-risk investments. Commercial monitoring also supports accounting that storage has occurred and is tied to project financing. It deals with long time frames and large volumes injected into multiple wells and multiple projects in favorable areas. We see developing trends toward reproducible workflows that systematically reduce risks and clarify expectations for oversight and long-term surveillance. Monitoring techniques showing increasing trends include injection zone pressure as a history-matching and compliance tool. To reduce cost and environmental impact of time-lapse seismic data collection, deploying new approaches and tools, such as use of fibre and installed sources are increasingly applied. Concern over the risk of induced seismicity by regulatory bodies and the general public has increased, which has also resulted in increased monitoring. Some techniques used in the early research phases have been sidelined or used only in restricted applications. For example, geochemical analyses in the injection zone as well as the environment are now being deployed less than it was in research-oriented programs, except in the US where it is required by the permitting process. Expectations of frequent area-wide near surface monitoring have also decreased.
The project Characterization of Most Promising Sequestration Formations in the Rocky Mountain Region 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 Rocky Mountain Carbon Capture and Storage (RMCCS) project investigated multiple geologic formations and characterized a local site on the Colorado Plateau for future CCS opportunities. The RMCCS project focused on the Cretaceous Dakota, Jurassic Entrada, and Pennsylvanian Weber Sandstones, the three largest regional formations. All formations in this project are potential CO2 storage resources for future power plants, natural gas processing plants, cement plants, and oil shale development projects. The area adjacent to Craig, Colorado, (Sand Wash Basin) was the area selected for detailed geologic characterization on the RMCCS project. The basin was selected in part because the geology can be extrapolated to other sites on the Colorado Plateau. Field mapping and seismic surveys were conducted to identify and evaluate the basin's structural configuration. A 9,745-foot deep characterization well was drilled to collect 131 feet of core and a suite of geophysical well log data. Petrophysical tests on samples of core were used to calibrate geophysical log data, which can be used to obtain storage resource estimates and evaluate associated uncertainty as well as simulate the hydrologic behavior of injected CO2. A detailed analysis of the primary formations (Dakota, Entrada and Weber sandstones) yielded a more accurate CO2 storage resource assessment for these formations within the Colorado Plateau; RMCCS estimates indicate a total CO2 storage resource of more than 38,000 million metric tons. The characterization of the Sand Wash Basin (2-D seismic surveys, multiple well logs and lithological, petrophysical and geochemical analyses) allowed for a detailed 3-D model to be constructed. The model served as the framework for analyses ranging from CO2 storage resource, injectivity, and subsurface flow to uncertainty estimates to evaluation of risk.
Within the NIF (National Ignition Facility) sustainment project, LLNL is replacing all 1,728 blast shields with new ones. Before the new blast shield glass can be installed on NIF, it must pass through the AMOL (Amplifier Main Optics Loop). The AMOL process involves cleaning and coating the glass, assembling it into a Line Replaceable Unit (LRU), and performing testing and inspection. This work will be conducted in the OPFX cleanroom facility (Building 391, Room 1250) and requires several pieces of custom equipment. One key piece of equipment within the loop is the Handling Frame Assembly Table (see Figure 1.1). Its purpose is to assist operators in installing the blast shield glass into handling frames (1014442709) (see Figure 1.2). The table is also used to store handling frame assembly components. The primary components stored are the L-frames (1015208271), which are kept in a drawer underneath the table (see Figure 1.3). The purpose of this safety note is to ensure that the system does not pose a risk of injury to personnel during a seismic event. The system owner is OMST Engineering & Maintenance.
While probabilistic risk assessment (PRA) of nuclear facilities is expected to include internal and external hazards for a risk-informed and performance-based design, the current state of practice treats each hazard independently. However, such an independent treatment of hazards may not account for the correlations between different hazards and their response of and damage to the structures, systems, and components (SSCs) in a plant resulting in underestimating the overall risk. This project proposes to advance the multi-hazard PRA of nuclear facilities to more adequately evaluate concurrent hazards and contribute to an increased safety of nuclear plants. A framework for multi-hazard PRA will be developed by identifying concurrent hazard events (both internal and external) and event sequences that include interdependencies through the response of SSCs. An example application of the multi-hazard PRA framework will be demonstrated by considering a generic pressurized water reactor (PWR) subjected to seismic and internal flooding hazards. Computational models for the response of components will be developed to generated multi-hazard fragility surfaces under seismic and flooding loads. A PRA model consisting of event and fault trees will also be developed to quantify the multi-hazard risk profile and compare it with the independent hazard risk profile. Overall, by advancing the multi-hazard PRA of nuclear facilities, this project enhances nuclear safety and reduces costs by mitigating unforeseen consequences caused by correlations between concurrent hazards.
The pursuit of knowledge and curiosity are once again driving humanity to seek a human presence on extraterrestrial bodies. While numerous challenges must be overcome on the moon to achieve this goal, the moonquake hazard has often been overlooked. Using NASA’s Moon to Mars Architecture requirements for surface habitation and accounting for the midterm/long-term mission duration, a framework for vulnerability assessment is developed for essential nonstructural elements (NSEs) inside an inflatable habitat subjected to moonquakes. The aim of this study is to develop an approach to assess the vulnerability of typical NSEs that support essential equipment, such as the environmental control and life support system, under a paucity of information related to the seismic hazard. We also emphasize that the launch dynamic environment and lunar seismic environment exhibit notable differences in their characteristics, which can lead to higher lateral acceleration levels than initially expected. This acceleration may be linked to high-occurrence seismic events. However, due to the lack of information regarding the frequency of these events, the risk of NSEs being subjected to unforeseen loads increases. It is imperative to equip future lunar habitats with seismic mitigation measures until additional seismometers are deployed near upcoming surface human outposts.
Understanding creep is essential for assessing seismic hazards along the creeping Hayward Fault, which threatens the San Francisco area. Current 1-D monitoring methods, like creepmeters and theodolites, provide limited measurements with low spatial resolution. InSAR offers higher 2-D spatial sampling but lower temporal resolution. We introduce a novel method using coda wave interferometry (CWI) to monitor creep in 2-D with high resolution. By calculating changes in ambient seismic noise travel times, we relate these measurements to creep offsets, bridging the gap between existing methods. Applied to seismic data from 2021 around the Hayward Fault, our CWI-derived offsets aligned with creepmeter measurements, revealing continuous and episodic creep. A 2-D regionalization identified a spatial boundary between these regimes, corroborating previous findings. This approach could enhance seismic monitoring and improve understanding of fault mechanics and earthquake risks.
Poster on “FRACML: A Machine Learning Based Tool to Quantify Reservoir Scale Fracture Network for CO2 Storage” for the CCUS 2025 conference held in Houston, Texas March 3-5, 2025. The accurate characterization of subsurface fracture networks is essential for the secure operation of carbon capture, utilization, and storage (CCUS) projects. A thorough understanding of the spatial distribution of subsurface faults and fractures is crucial for predicting CO2 plume evolution and minimizing risks such as potential leakage into overlying formations or induced seismicity. In this context, robust fracture network quantification plays a pivotal role in reservoir management, providing the data necessary to fine-tune operational parameters, and ensure the environmental and economic viability of CCUS projects. As part of the U.S. Department of Energy’s SMART (Science-informed Machine Learning for Accelerating Real-time Decisions in Subsurface Applications) initiative, we focused on the development and application of a machine learning-based tool (FRACML) designed to quantify and map fracture networks using real-world (non-synthetic) data from an active CO2 injection site. Our objective is to demonstrate the utility of this tool in improving operational efficiency and safety across CCUS sites.
Enhanced Geothermal Systems (EGS) produce clean energy by circulating fluid through hot rock deep underground and bringing that heat to the surface to generate electricity. For this process to work reliably, fluids must be able to move efficiently through networks of natural or engineered fractures in the rock. Enhancing and maintaining subsurface permeability over time is essential for sustainable energy production. However, fluid injection changes the underground temperature, pressure, rock stress, and chemistry, which can alter permeability and sometimes trigger earthquakes. Predicting these interconnected processes remains a key challenge. To address this, we combined high-temperature laboratory experiments with high-fidelity simulations to better understand how fractures in geothermal reservoirs evolve over time. Our experiments measured how fractures respond to stress, slip, slip rate, and chemical reactions under geothermal conditions. These data were integrated into coupled thermal-hydrological-mechanical-chemical and earthquake (THMC+E) models tailored to the Utah FORGE site. The validated modeling framework improves predictions of reservoir performance and seismic response and helps guide operational decisions. This work reduces technical risk and strengthens the scientific foundation needed to make geothermal energy a reliable and scalable clean energy resource.
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.