Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geologic Modeling”

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

Enhanced Multi-Dimensional Inversion Through Target-Specific Inversion Parameter Bounds With an Application to Crosswell Electromagnetic for Sequestration Monitoring

In geophysical inversions, lower and upper model parameter bounds are a means of solution stabilization. Further, constraints that intend to let only geologically plausible inverse solutions pass are amenable to lower and upper bounds. Reliable prior information is paramount to construct such bound constraints. It is common practice to narrow and widen bound intervals for regions of, respectively, more and less certain prior information. Contrary to this practice, we experiment with widened bound intervals in zones that are poorly resolved by a given survey configuration but where prior information would suggest structural anomalies of interest. The purpose of enlarged parameter bounds that correlate spatially with predefined targets is to let the inversion explore a larger solution space, thus increasing the potential to resolve otherwise hidden anomalies. Application of the method is based on a carbon-sequestration baseline (pre-injection) crosswell electromagnetic (EM) field survey at the Containment and Monitoring Institute Field Research Station (Alberta, Canada), where impeded measurements led to generally reduced sensitivities for the interwell region. Synthetic-data proofs of concept use augmented bounds designed to boost the resolution of artificial plume targets, indicating an enhanced illumination compared to constant bounds. Comparative field data inversions with spatially variable bounds constructed from prior resistivity and velocity information highlight non-horizontal baseline structures.

3-D electrical resistivity imaging↗
DOE OSTI2025

Learning the factors controlling mineral dissolution in three-dimensional fracture networks: applications in geologic carbon sequestration

We perform a set of high-fidelity simulations of geochemical reactions within three-dimensional discrete fracture networks (DFN) and use various machine learning techniques to determine the primary factors controlling mineral dissolution. The DFN are partially filled with quartz that gradually dissolves until quasi-steady state conditions are reached. At this point, we measure the quartz remaining in each fracture within the domain as our primary quantity of interest. We observe that a primary sub-network of fractures exists, where the quartz has been fully dissolved out. This reduction in resistance to flow leads to increased flow channelization and reduced solute travel times. However, depending on the DFN topology and the rate of dissolution, we observe substantial variability in the volume of quartz remaining within fractures outside of the primary subnetwork. This variability indicates an interplay between the fracture network structure and geochemical reactions. We characterize the features controlling these processes by developing a machine learning framework to extract their relevant impact. Specifically, we use a combination of high-fidelity simulations with a graph-based approach to study geochemical reactive transport in a complex fracture network to determine the key features that control dissolution. We consider topological, geometric and hydrological features of the fracture network to predict the remaining quartz in quasi-steady state. We found that the dissolution reaction rate constant of quartz and the distance to the primary sub-network in the fracture network are the two most important features controlling the amount of quartz remaining. This study is a first step towards characterizing the parameters that control carbon mineralization using an approach with integrates computational physics and machine learning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Radionuclide Transport Simulations Supporting Proposed Borehole Waste Disposal in Israel

A scientific collaboration between the U.S. and Israel is underway to assess the suitability of a potential site for subsurface radioactive waste disposal in the Negev Desert, Israel. The Negev Desert has several favorable attributes for geologic disposal, including an arid climate, a deep vadose zone, interlayered low-permeability lithologies, and carbonate rocks with high uranium-sorption potential. These features may provide a robust natural barrier to radionuclide migration. Geologic and laboratory characterization data from the Negev Desert are incorporated into multiphase flow and transport models, solved using PFLOTRAN, to aid in site characterization and risk analysis that will support decision-making for waste disposal in an intermediate-depth borehole design. The lithology with the greatest uranium sorption potential at the site is phosphorite. We use modeling to evaluate the ability of this layer to impact uranium transport around a proposed disposal borehole. The current objective of the simulations is focused on characterizing hypothetical leakage from waste canisters and subsequent uranium migration under three infiltration scenarios. Here, we describe a hydrogeologic model based on data from a local exploratory borehole and present results for uranium flow and transport simulations under varying infiltration scenarios. We find that under the current climate conditions, it is likely that uranium will remain in the near-field of the borehole for thousands of years. However, under a hypothesized extreme climate scenario representing an increase in infiltration by a factor of 300x above present-day values, uranium may break through the phosphorite layer and exit the base of the model domain (~200 m above the water table) within 1000 years. Simulation results have direct implications for the planning of nuclear waste disposal in the Negev Desert, and specifically in intermediate-depth boreholes.

58 GEOSCIENCES↗
DOE OSTI2021

Hydraulic and Thermal Stimulation Program at Raft River Idaho, A DOE EGS

A Department of Energy Enhanced Geothermal System (EGS) stimulation program has injected over 254 million gallons of water into the well RRG-9 ST1 since the summer of 2013. Three major stimulations have been conducted during the program increasing injection flow rates from less than 20 gpm to 550 gpm. Geologic, water chemistry, microseismic activity, and borehole imaging data have been used to develop a conceptual model describing possible flow paths of this injected water. This model contains two major fracture zones one of which intersects the RRG-9 ST1 wellbore. Modified Hall and injectivity index plots constructed using injection flow rates, surface temperatures, and wellhead pressures show steady improvement in the injectivity of the well. Here, the injectivity index has risen from 0.15 gpm/psi to 2.0 gpm/ psi. A pressure falloff test conducted on April 28, 2015 indicates a reservoir permeability of 1,220 md and -5.38 skin factor. The well stimulation program was simulated numerically using an Idaho National Laboratory reservoir simulation code, FALCON. These simulations show a significant increase in the permeability of connecting fracture pathways after each stimulation event.

Enhanced geothermal system↗
DOE OSTI2024

Nitrogen Vacancy Centers in Diamond for Stress Sensing Applications: Theory and Experiments

Sensing of stress under high elevated environmental conditions with high resolution has critical importance for range of applications including earth’s subsurface scanning and geological CO2 storage monitoring. Using first principles density functional theory (DFT) approach combined with the theoretical modelling of low energy Hamiltonian, here we investigate a novel approach to detect unprecedented level of pressure by taking advantage of solid-state electron’s spin of Nitrogen vacancy (NV) centers in nanodiamond. We computationally explore the effect of strain on the defect band edges and band gaps by varying the lattice parameters of diamond supercell hosting a single NV center. A low energy Hamiltonian is developed that includes effect of stress on the energy level of ±1 spin manifold at the ground state. By quantifying the energy level shift and split, we predict the pressure sensing of up to 0.3 MPa/√Hz. We show superiority of the quantum sensing over traditional optical sensing techniques by discussing our results from DFT and theoretical modelling for the frequency shift per unit pressure. The proposed quantum stress sensing could be useful to measure earth’s subsurface vibrations. Our results open avenues for development of sensing technology with a high sensitivity and resolution under extreme pressure conditions that potentially has a wider applicability than existing pressure sensing technologies.

Paudel, Hari P.↗
DOE OSTI2023

Prediction of Long-Term Geochemical Change in Bentonite Based on the Interpretative THMC Model of the FEBEX In Situ Test

Since nuclear energy is crucial in the decarbonization of the energy supply, one hurdle to remove is the handling of high-level radioactive waste (HLW). Disposal of HLW in a deep geological repository has long been deemed a viable permanent option. In the design of a deep geological repository, compacted bentonite is the most commonly proposed buffer material. Predicting the long-term chemical evolution in bentonite, which is important for the safety assessment of a repository, has been challenging because of the complex coupled processes. Models for large-scale tests and predictions based on such models have been some of the best practices for such purposes. An 18-year-long in situ test with two dismantling events provided a unique set of chemical data that allowed for studying chemical changes in bentonite. In this paper, we first developed coupled thermal, hydrological, mechanical, and chemical (THMC) models to interpret the geochemical data collected in the in situ test and then extended the THMC model to 200 years to make long-term prediction of the geochemical evolution of bentonite. The interpretive coupled THMC model shows that the geochemical profiles were strongly affected by THM processes such as evaporation/condensation, porosity change caused by swelling, permeability change, and the shape of concentration profiles for major cations were largely controlled by transport processes, but concentration levels were regulated by chemical reactions, and the profiles of some species such as pH, bicarbonate, and sulfate were dominated by these reactions. The long-term THMC model showed that heating prolongs the time that bentonite becomes fully saturated in the area close to the heater/canister; however, once the bentonite becomes fully saturated, high concentrations of ions in bentonite near the heater, which was observed in the field test, will disappear; illitization continues for 50 years but will not proceed further.

Zheng, Liange↗
DOE OSTI2025

Streamflow Response to Glacier Mass Loss Varies With Basin Precipitation Across Alaska

Abstract Diminishing glaciers affect streamflow, and given the extent of glaciers in Alaska and adjacent Canada, continued glacier mass loss is likely to have profound effects on ecosystems sensitive to runoff. The effects of glacier mass loss on streamflow are likely to vary across the wide ranges of basin size, glacier cover, and precipitation in this region. In this study, we use U.S. Geological Survey (USGS) streamflow data with satellite‐based glacier volume change estimates to quantify how glacier mass loss subsidized streamflow over the 2000–2019 period for 116 glacierized basins. We examine interannual variability in that subsidy at three USGS‐monitored glaciers to explore the ability of the subsidy to buffer streamflow derived solely from precipitation. We found the relative importance of percent glacier cover on streamflow magnitude increases in drier basins. In the driest basins, glaciers produced 40 times greater percent glacier mass loss subsidies to streamflow for the percent glacier cover compared to the wettest basins. While the subsidy from glacier mass loss buffers interannual variability in streamflow to varying degrees, it can also increase streamflow variability. Smaller amounts of percent glacier cover are needed to produce summer‐melt‐dominated seasonal flow regimes in drier basins than in wetter basins. Decreasing glacier cover will eventually decrease summer streamflow, increasing spring streamflow in drier basins, and attenuating seasonality with increasing spring and autumnal streamflow in wetter basins. Quantifying the downstream effects of continued glacier mass loss without the computational expense of a hydrological model is broadly applicable in this changing climate.

Curran, Janet H. [U.S. Geological Survey Alaska Sc↗
DOE OSTI2023

Continental Scale Hydrostratigraphy: Basin‐Scale Testing of Alternative Data‐Driven Approaches

Abstract Integrated hydrological modeling is an effective method for understanding interactions between parts of the hydrologic cycle, quantifying water resources, and furthering knowledge of hydrologic processes. However, these models are dependent on robust and accurate datasets that physically represent spatial characteristics as model inputs. This study evaluates multiple data‐driven approaches for estimating hydraulic conductivity and subsurface properties at the continental‐scale, constructed from existing subsurface dataset components. Each subsurface configuration represents upper (unconfined) hydrogeology, lower (confined) hydrogeology, and the presence of a vertical flow barrier. Configurations are tested in two large‐scale U.S. watersheds using an integrated model. Model results are compared to observed streamflow and steady state water table depth (WTD). We provide model results for a range of configurations and show that both WTD and surface water partitioning are important indicators of performance. We also show that geology data source, total subsurface depth, anisotropy, and inclusion of a vertical flow barrier are the most important considerations for subsurface configurations. While a range of configurations proved viable, we provide a recommended Selected National Configuration 1 km resolution subsurface dataset for use in distributed large‐and continental‐scale hydrologic modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

User’s Manual for Seal_Flux: A Seal Barrier Reduced-Order Model (Update)

This report provides a brief description on the use of the Seal_Flux computer program developed as part of the effort to quantify the risk of geologic storage of carbon dioxide (CO 2 ) under the U.S. Department of Energy’s (DOE) National Risk Assessment Partnership (NRAP). The Seal_Flux code simulates the flow of CO 2 through a low permeability rock horizon or seal formation overlying the storage reservoir into which CO 2 is injected. A two-phase, relative permeability approach with Darcy’s law is used for one-dimensional (1D) flow computations of CO 2 through the horizon in the vertical direction. The code also allows the simulation of time-dependent processes that can influence such flow. However, as part of its design, Seal_Flux is what can be termed a “reduced-order model” (ROM) and is not intended as a full-functioning flow code. The theory and simulation in the code is streamlined and directed towards the implementation of Monte Carlo risk analyses of CO 2 transport or as termed in this context as “leakage.” While presented in this report as a stand-alone tool, the Seal_Flux code is intended to function in the future as one of several models as part of an integrated, systems-level model of CO 2 storage performance. Finally, the code is written in Python 3.10 to provide an open framework for further development by others and to assist in linking the code with other modules in an integrated assessment model.

58 GEOSCIENCES↗
DOE OSTI2023

Reservoir Thermal Energy Storage Benchmarking (Rev. 3)

A benchmarking analysis of RTES research funded by GTO through the Beyond Batteries projects was conducted against the ESGC to see where they fit within the identified ESGC Use Cases. The projects were found to advance knowledge in multiple ESGC use cases, either directly or in some cases, indirectly as enabling technologies. This analysis is helpful to understand where RTES and associated research fits into the larger discussion around energy storage technologies. Also, a retrospective analysis of the Beyond Batteries projects was conducted to evaluate what the projects learned and how the results can be applied to advance the value of RTES. Major results of each of the studies are summarized in Table 2. Additionally, a comparative metrics analysis for RTES was completed to understand where RTES lies within the energy storage industry. Metrics for evaluation of RTES and its comparison to other storage technologies were selected and ranges of their values compiled. The selected metrics – LCOE (levelized cost of energy), capital costs, roundtrip efficiency, energy storage capacity, and storage time – were chosen based on data availability and have a particularly strong influence on the potential deployment of a storage technology. Charts which compare the metrics are presented in section 4.3 and show ranges for each of the 10 selected technologies. However, due to a lack of domestic operational facilities, values for RTES and for portions of the remaining technologies are based on theoretical modeling and studies of best-case scenarios. LCOE estimates for RTES fall within the lower reaches of Figure 15, but nevertheless amount to 2 – 5 times the ESGC Roadmap goal for LCOE, for example in the Facilitating and Evolving Grid Use Case. Capital costs for RTES sit on the higher end (Figure 16) but are expected to decrease as new projects are developed and the technology is refined. The theoretical roundtrip efficiency reported for RTES varies from mid to high percentages (Figure 17) with efficiencies upwards of 93% in modeled scenarios in the Portland Basin (Bershaw et al.,2020). RTES is also expected to have the largest energy storage capacities and longest storage times, likely matched only by lower efficiency hydrogen storage. To better assess the role that RTES could play in energy storage we examined it’s potential in the U.S. The potential depends on many factors. Recently, many researchers have started looking at deep sedimentary basins, depleted oil and gas fields, and basalt formations as potential targets for RTES development. The United States Geological Survey (USGS) has analyzed various cities and shown substantial RTES potential in the cooling sector (Pepin et al., 2021). By modeling RTES in low-quality groundwater (e.g., brackish), it is shown to be favorable across the U.S. with particular suitability in the Illinois Basin, Coastal Plains, and Basin and Range regions. Seasonal RTES operations have also been modeled in the Portland Basin by those at the USGS and Portland State University to simulate an RTES system supplying heating loads needed for the Oregon Health and Science University. Simulations suggest that high conductive heat loss in the initial years exists but tends to decrease with increasing time and development of the resource due to self-insulating nature of the basalts (Burns et al., 2020). Other national laboratory efforts are taking a close look at many of the technical issues involved with RTES (McLing et al., 2019, McLing et al., 2022). These include difficulties in understanding geochemical, hydrogeological, mechanical, and microbiological changes at such elevated temperatures and operational scenarios. Major gaps in research are identified and suggested for future work. With this increased focus to understand how to make RTES successful in the U.S., this technology could be a potential solution to many of the nation’s energy storage problems. For the energy independence of this country, the DOE should prioritize de-risking this technology by making future investments in pilot-scale demonstrations to attract potential investors.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

FECM/NETL CO 2 Saline Storage Cost Model (2024): User’s Manual

The U.S. Department of Energy's (DOE) Office of Fossil Energy and Carbon Management (FECM), in collaboration with the National Energy Technology Laboratory (NETL), has developed the FECM/NETL CO 2 Saline Storage Cost Model (CO2_S_COM). This Excel-based tool provides a comprehensive framework for estimating the costs and breakeven prices associated with storing carbon dioxide (CO 2 ) in deep saline formations. Designed from the perspective of a CO 2 storage site owner, the CO2_S_COM incorporates four integrated modules—project management, financial analysis, activity cost estimation, and geological evaluation—to deliver fast, robust, and actionable insights for evaluating project finances. This is the user's manual for CO2_S_COM. The model may be accessed at this link: FECM/NETL CO2 Saline Storage Cost Model CO2_S_COM 2024 (v4) - Submissions - EDX

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

On the dissolution of a borosilicate glass with the use of isotopic tracing – Insights into the mechanism for the long-term dissolution rate

An understanding of the mechanisms responsible for controlling the long-term corrosion rate of nuclear waste glass is paramount if reliable glass dissolution models are to be used to calculate the controlled release of radionuclides from nuclear waste glass under geologic disposal conditions. Additionally, understanding silicate glass dissolution is also important for natural glasses to understand the role from the dissolution of these glasses has on the composition of natural aquatic systems. Two general mechanisms appear to be responsible for the elemental release from glass to the surrounding biosphere - ion exchange (release principally of alkalis) and matrix dissolution (release of the structural components of the glass). Key unknowns related to these mechanisms are the impact of surface layers on the altered glass and how these layers participate in corrosion. To better understand the impact of these layers on glass corrosion, we made two glasses of the same elemental composition, but with selected elements enriched in specific isotopes, e.g. more 10 B than 11 B. We reacted these glasses for 1 y then the solutions from each were exchanged and the chemistry of the solutions and the solids were followed for another 2 y and in some cases 3 y. This allowed us to follow changes in the glasses with respect to these elements and how these elements were transported through the alteration layers with time and, for at least one element, into the glass. Results from these experiments demonstrate that the release behavior of different elements is strongly dependent on their structural role in the glass (e.g. network formers and modifiers), but, more importantly, the role of water transport and subsequent ion exchange has in the long-term dissolution of the glass studied here. In this article, we highlight the behavior of four elements: lithium, sodium, silicon, and boron. Lithium and sodium, network modifying elements, have similar chemistries. The behavior of Li could be tracked in more detail by following both 6 Li and 7 Li; Li penetrates through the gel layer in both directions without hinderance and into and released from the glass deeper than sodium. Silicon, a network forming element, reacts with the silica-rich alteration layer. Boron, a network forming element, does not accumulate in the gel or the pristine glass and has a very sharp elemental profile between the pristine glass and alteration layer. Boron, lithium, and sodium elemental profiles suggest that there is little if any transport control within the alteration layer. The correlation between the Li and Na profiles and that of a water species suggest that the limiting release is a steady state between the diffusion of a water species and matrix dissolution that results in a low release rate for alkali and this provides a steady driving force for the long-term dissolution of glass.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Deep learning-accelerated 3D carbon storage reservoir pressure forecasting based on data assimilation using surface displacement from InSAR

Fast forecasting of the reservoir pressure distribution during geologic carbon storage (GCS) by assimilating monitoring data is a challenging problem. Due to high drilling cost, GCS projects usually have spatially sparse measurements from few wells, leading to high uncertainties in reservoir pressure prediction. To address this challenge, we use low-cost Interferometric Synthetic-Aperture Radar (InSAR) data as monitoring data to infer reservoir pressure build up. We develop a deep learning-accelerated workflow to assimilate surface displacement maps interpreted from InSAR and to forecast dynamic reservoir pressure. Employing an Ensemble Smoother Multiple Data Assimilation (ES-MDA) framework, the workflow updates three-dimensional (3D) geologic properties and predicts reservoir pressure with quantified uncertainties. We use a synthetic commercial-scale GCS model with bimodally distributed permeability and porosity to demonstrate the efficacy of the workflow. A two-step CNN-PCA approach is employed to parameterize the bimodal fields. The computational efficiency of the workflow is boosted by two residual U-Net based surrogate models for surface displacement and reservoir pressure predictions, respectively. The workflow can complete data assimilation and reservoir pressure forecasting in half an hour on a personal computer.

25 ENERGY STORAGE↗
DOE OSTI2020

Integration of Soft Data Into Geostatistical Simulation of Categorical Variables

Uncertain or indirect “soft” data, such as geologic interpretation, driller’s logs, geophysical logs or imaging, offer potential constraints or “soft conditioning” to stochastic models of discrete categorical subsurface variables in hydrogeology such as hydrofacies. Previous bivariate geostatistical simulation algorithms have not fully addressed the impact of data uncertainty in formulation of the (co) kriging equations and the objective function in simulated annealing (or quenching). This paper introduces the geostatistical simulation code tsim-s, which accounts for categorical data uncertainty through a data “hardness” parameter. In generating geostatistical realizations with tsim-s, the uncertainty inherent to soft conditioning is factored into both 1) the data declustering and spatial correlation functions in cokriging and 2) the acceptance probability for change of category in simulated quenching. The degree or sensitivity to which soft data conditions a realization as a function of hardness can be quantified by mapping category probabilities derived from multiple realizations. In addition to point or borehole data, arrays of data (e.g., as derived from a depth-dependency function, probability map, or “prior realization”) can be used as soft conditioning. The tsim-s algorithm provides a theoretically sound and general framework for integrating datasets of variable location, resolution, and uncertainty into geostatistical simulation of categorical variables. A practical example shows how tsim-s is capable of generating a large-scale three-dimensional simulation including curvilinear features.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Improving Building Footprint Extraction Using NAIP and 3DEP Lidar Derived Features with Deep Learning

Accurate building footprint extraction is critical for applications ranging from population estimation to disaster management. Although optical imagery provides detailed spectral information, it often struggles with shadows, occlusions, and background clutter in dense urban environments. Lidar data, by contrast, offer precise elevation and structural attributes but face challenges such as variable point density and noise. This study integrates multispectral imagery from the U.S. Department of Agriculture (USDA) National Agriculture Imagery Program (NAIP) with lidar-derived feature height and intensity from the U.S. Geological Survey (USGS) 3D Elevation Program (3DEP) to improve footprint extraction using a U-Net–based deep learning model. A six-band input stack (RGB, near-infrared, height, intensity) was developed, normalized, and tiled for training and evaluation against Microsoft Global Building Footprints (GBF). Results from the Houston, TX test site show that the six-band model achieved a precision of 0.86, recall of 0.88, F1 score of 0.87, and Intersection-over-Union (IoU) of 0.76, consistently outperforming four-band baselines by reducing false positives while maintaining sensitivity. Predictions on withheld Houston tiles confirmed strong within-region generalization, yielded a precision of 0.78, recall of 0.81, F1 score of 0.79, and IoU of 0.66. Qualitative analysis further revealed limitations stemming from both training label quality and vegetation–building confusion. These findings demonstrate the complementary value of integrating spectral and structural information for robust building footprint extraction and how domain adaptation strategies can be used to enhance cross-regional transferability.

Liu, Jung Kuan [United States Geological Survey (U↗
DOE OSTI2024

Modeling injection-induced fault slip using long short-term memory networks

Stress changes due to changes in fluid pressure and temperature in a faulted formation may lead to the opening/shearing of the fault. This can be due to subsurface (geo)engineering activities such as fluid injections and geologic disposal of nuclear waste. Such activities are expected to rise in the future making it necessary to assess their short- and long-term safety. Here, a new machine learning (ML) approach to model pore pressure and fault displacements in response to high-pressure fluid injection cycles is developed. The focus is on fault behavior near the injection borehole. To capture the temporal dependencies in the data, long short-term memory (LSTM) networks are utilized. To prevent error accumulation within the forecast window, four critical measures to train a robust LSTM model for predicting fault response are highlighted: (i) setting an appropriate value of LSTM lag, (ii) calibrating the LSTM cell dimension, (iii) learning rate reduction during weight optimization, and (iv) not adopting an independent injection cycle as a validation set. Several numerical experiments were conducted, which demonstrated that the ML model can capture peaks in pressure and associated fault displacement that accompany an increase in fluid injection. The model also captured the decay in pressure and displacement during the injection shut-in period. Further, the ability of an ML model to highlight key changes in fault hydromechanical activation processes was investigated, which shows that ML can be used to monitor risk of fault activation and leakage during high pressure fluid injections.

58 GEOSCIENCES↗
DOE OSTI2024

Ground-motions site and event specificity: Insights from assessing a suite of simulated ground motions in the San Francisco Bay Area

This article presents the results of a research that is part of a larger collaborative effort between the Lawrence Berkeley National Laboratory and the Pacific Earthquake Engineering Research Center, funded by the US Department of Energy Office of Cybersecurity, Energy Security and Emergency Response. The main objective of this study is to assess a suite of near and far-field simulated ground motions obtained from 20 realizations of an M7 Hayward Fault earthquake in the San Francisco Bay Area, California USA, and inform the selection of rupture simulation parameters leading to strong motions. To this aim, comparisons are conducted with NGA-W2 and directivity ground-motion models and a selected population of records. An archetypal steel moment-resisting frame is utilized to assess infrastructure response distributions. The analyses carried out for each simulated event and subdomain with consistent properties in terms of shallow shear-wave velocity proved to be instrumental for better interpreting the differences between simulated motions and empirical models. The main reasons identified for variances between simulations and empirical relationships included (1) directivity effects fully captured by the simulations across the full breadth of rupture models; (2) site vicinity to ruptures that incorporate large-slip patches, particularly if these are in the forward-directivity direction; and (3) presence of geologic structures that can “trap” seismic waves and produce ground motions with large amplitude and long signal duration. The analyses carried out in this work provide a path for interpreting ground-motion site and event specificity obtained from a suite of physics-based simulations, differing only in the rupture model characterization, to inform the selection of simulation scenarios for site-specific engineering analyses under strong excitations. Evidence from this work points to the possibility that current hazard models may underestimate ground-motion intensities in areas where the combined effect of directivity and site conditions results in large ground-motion amplitudes.

58 GEOSCIENCES↗
DOE OSTI2024

3D Quantification of Pore Networks and Anthropogenic Carbon Mineralization in Stacked Basalt Reservoirs

Basalt formations are promising candidates for the geologic storage of anthropogenic CO 2 due to their storage capacity, porosity, permeability, and reactive geochemical trapping ability. The Wallula Basalt Carbon Storage Pilot Project demonstrated that supercritical CO 2 injected into >800 m deep Columbia River Basalt Group stacked reservoir flow tops mineralizes to ankerite-siderite-aragonite on month-year timescales, with 60% of the 977 metric tonnes of CO 2 converted within two years. The potential impacts of mineral precipitation and consequent changes on rock porosity, pore structure, pore size, and pore size distributions have likely been underestimated hitherto. Herein, we address these knowledge gaps using X-ray Micro Computed Tomography (XMT) to evaluate the pore network architecture of sidewall cores recovered two years after CO 2 injection. In the present study, we performed a detailed quantitative analysis of the CO 2 -reacted basalt cores by XMT imaging. Reconstructed 3D images were analyzed to study the distribution and volumetric details of porosity and carbonate nodules in the cores along with the various other phases, providing insight into paragenesis and carbonate growth mechanisms, including mineralogic/chemical zonation. Finally, these findings are being used to parametrize multiphase reactive transport models to predict the fate and transport of subsurface CO 2 , enabling scale-up to commercial-scale geologic carbon storage in basalts and other reactive mafic-ultramafic formations.

54 ENVIRONMENTAL SCIENCES↗