SEARCH · Engineering Papers

Results for “Geologic uncertainty”

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.

SourceResearch connectionFrom yearThrough year

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

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Impacts and tectonism in Earth and moon history of the past 3800 million years

The moon's surface, unlike the Earth's, displays a comparatively clear record of its past bombardment history for the last 3800 Myr, the time since active lunar tectonism under the massive premare bombardment ended. From Baldwin's (1987) tabulation of estimated ages for a representative sample of large lunar craters younger than 3800 Ma, six major cratering episodes can be discerned. These six bombardment episodes, which must have affected the Earth too, appear to match in time the six major episodes of orogenic tectonism on Earth, despite typical resolution errors of +/- 100 Myr and the great uncertainties of the two chronologies. Since more highly resolved events during the Cenozoic and Mesozoic Eras suggest the same correlation, it is possible that large impacts have influenced plate tectonics and other aspects of geologic history, perhaps by triggering flood basalt eruptions.

Stothers, Richard B.↗

DOE OSTI2021

Introducing the GeoRePORT Resource Size Tool: Reporting on Geothermal Resource Size Estimations Using the Geothermal Resource Portfolio Optimization and Reporting Technique (GeoRePORT): Preprint

The Geothermal Resource Portfolio Optimization and Reporting Technique (GeoRePORT) was developed with funding from the U.S. Department of Energy Geothermal Technologies Office to assist in identifying and pursuing long-term investment strategies through the development of a resource reporting protocol. The assessment protocols used in GeoRePORT allow for comparison of project attributes across locations and geological settings to understand the feasibility of geothermal development. This work introduces the Resource Size Tool, a new feature within the GeoRePORT package that compiles two independent methods for estimating geothermal resource size in terms of energy capacity in MW. Energy production potential for twenty-three case studies was estimated with the Resource Size Tool in order to 1) generate a reasonable range of resource size estimates for a particular geothermal field; 2) illustrate the advantages and limitations of each methodology (such as data input requirements, estimate accuracy and precision, and the appropriate circumstances of use); and 3) test the ability of the resource size tool to provide useful and accurate information for geothermal stakeholders. The tool employs two methods widely used in the geothermal industry: (1) USGS Volumetric and (2) Power Density. Results from our case studies show general overlap between these two methods in terms of resource size estimates; however, they also reveal key differences between the two approaches that should be considered when using such estimates to drive development. First, the two methods rely on different input parameters and therefore one method may be more appropriate and/or accurate for a given project than the other. Second, the Power Density method was found to generate wider ranges of resource size predictions, more consistently aligning with actual power production of the field but with larger scales of error; whereas the USGS Volumetric method predicts narrower ranges but tends to overestimate when compared to current MW production. Future work will refine variables used in the methods with input data from other sections of GeoRePORT and modify uncertainty levels based on the particular datasets used for a given project.

geological↗

DOE OSTI2023

Contributions of Climate Change and ENSO Variability to Future Precipitation Extremes Over California

The El Niño/Southern Oscillation (ENSO) affects the occurrence frequency and intensity of extreme precipitation through modulations of regional heat and moisture fluxes. California experiences particularly strong ENSO influences and models project different to its extreme precipitation. It remains unclear how diverse projections of future precipitation extremes relate to inter-model differences for those changing signals. Here, we use “large ensemble” simulations with multiple climate models along with the Coupled Model Intercomparison Project Phase 5 to investigate the range of precipitation extreme changes over California and the influences from ENSO-related teleconnections. We found that precipitation amount increases are much larger during El Niño relative to La Niña years, mainly caused by the differences in frequency of extreme events during different phases. The ENSO-driven effect is even larger than the overall climate change signal for the most extreme events, implying uncertainties from inter-model differences in ENSO-related SST variability for extreme precipitation changes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

River Deltas and Sea-Level Rise

Future sea-level rise poses an existential threat for many river deltas, yet quantifying the effect of sea-level changes on these coastal landforms remains a challenge. Sea-level changes have been slow compared to other coastal processes during the instrumental record, such that our knowledge comes primarily from models, experiments, and the geologic record. Here we review the current state of science on river delta response to sea-level change, including models and observations from the Holocene until 2300 CE. We report on improvements in the detection and modeling of past and future regional sea-level change, including a better understanding of the underlying processes and sources of uncertainty. We also see significant improvements in morphodynamic delta models. Still, substantial uncertainties remain, notably on present and future subsidence rates in and near deltas. Observations of delta submergence and land loss due to modern sea-level rise also remain elusive, posing major challenges to model validation. There are large differences in the initiation time and subsequent delta progradation during the Holocene, likely from different sea-level and sediment supply histories. Modern deltas are larger and will face faster sea-level rise than during their Holocene growth, making them susceptible to forced transgression. Regional sea-level projections have been much improved in the past decade and now also isolate dominant sources of uncertainty, such as the Antarctic ice sheet. Vertical land motion in deltas can be the dominant source of relative sea-level change and the dominant source of uncertainty; limited observations complicate projections. River deltas globally might lose 5% (~35,000 km 2 ) of their surface area by 2100 and 50% by 2300 due to relative sea-level rise under a high-emission scenario.

58 GEOSCIENCES↗

DOE OSTI2023

Endothermic physiology of extinct megatooth sharks

The evolution of the extinct megatooth shark, Otodus megalodon , and its close phylogenetic relatives remains enigmatic. A central question persists regarding the thermophysiological origins of these large predatory sharks through geologic time, including whether O. megalodon was ectothermic or endothermic (including regional endothermy), and whether its thermophysiology could help to explain the iconic shark’s gigantism and eventual demise during the Pliocene. To address these uncertainties, we present unique geochemical evidence for thermoregulation in O. megalodon from both clumped isotope paleothermometry and phosphate oxygen isotopes. Our results show that O. megalodon had an overall warmer body temperature compared with its ambient environment and other coexisting shark species, providing quantitative and experimental support for recent biophysical modeling studies that suggest endothermy was one of the key drivers for gigantism in O. megalodon and other lamniform sharks. The gigantic body size with high metabolic costs of having high body temperatures may have contributed to the vulnerability of Otodus species to extinction when compared to other sympatric sharks that survived the Pliocene epoch.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Combining Eddy Covariance Towers, Field Measurements, and the MEMS 2 Ecosystem Model Improves Confidence in the Climate Impacts of Bioenergy With Carbon Capture and Storage

ABSTRACT Carbon dioxide removal technologies such as bioenergy with carbon capture and storage (BECCS) are required if the effects of climate change are to be reversed over the next century. However, BECCS demands extensive land use change that may create positive or negative radiative forcing impacts upstream of the BECCS facility through changes to in situ greenhouse gas fluxes and land surface albedo. When quantifying these upstream climate impacts, even at a single site, different methods can give different estimates. Here we show how three common methods for estimating the net ecosystem carbon balance of bioenergy crops established on former grassland or former cropland can differ in their central estimates and uncertainty. We place these net ecosystem carbon balance forcings in the context of associated radiative forcings from changes to soil N 2 O and CH 4 fluxes, land surface albedo, embedded fossil fuel use, and geologically stored carbon. Results from long term eddy covariance measurements, a soil and plant carbon inventory, and the MEMS 2 process‐based ecosystem model all agree that establishing perennials such as switchgrass or mixed prairie on former cropland resulted in net negative radiative forcing (i.e., global cooling) of −26.5 to −39.6 fW m −2 over 100 years. Establishing these perennials on former grassland sites had similar climate mitigation impacts of −19.3 to −42.5 fW m −2 . However, the largest climate mitigation came from establishing corn for BECCS on former cropland or grassland, with radiative forcings from −38.4 to −50.5 fW m −2 , due to its higher plant productivity and therefore more geologically stored carbon. Our results highlight the strengths and limitations of each method for quantifying the field scale climate impacts of BECCS and show that utilizing multiple methods can increase confidence in the final radiative forcing estimates.

Falvo, Grant [Department of Plant, Soil and Microb↗

DOE OSTI2024

SCEC/USGS Community Stress-Drop Validation Study: How Spectral Fitting Approaches Influence Measured Source Parameters

Spectral source parameters used to estimate an earthquake’s stress drop (⁠Δσ⁠) can vary significantly across measurement approaches. The Statewide California Earthquake Center/U.S. Geological Survey Community Stress‐Drop Validation Study was initiated to compare source parameter estimates, focusing initially on a dataset from the 2019 Ridgecrest earthquake sequence. As part of that validation effort, here we focus on one potential source of uncertainty: whether spectral fitting approaches alone, applied to a common set of spectra from the 2019 Ridgecrest sequence result in different source parameter estimates. By using a common set of benchmark spectra analyzed across a consistent frequency band of 1–40 Hz, we eliminate many sources of variability. A subgroup of validation study participants volunteered to estimate the low‐frequency displacement (⁠Ω 0 ⁠) and corner frequency (⁠ƒ c ⁠) by fitting a smooth function to benchmark displacement spectra. Participants used linear‐ or log‐sampled spectra, assumed a Brune or Boatwright spectral model, and applied different misfit criteria. We compare 17 approaches used to estimate ⁠Ω 0 ⁠, ƒ c ⁠, and Δσ for 54 earthquake spectra. Our results reveal that 35% of events have Δσ estimates within a factor of two, whereas others exhibit variations exceeding an order of magnitude. The variability in and can largely be attributed to whether a spectrum is consistent with the smooth function of an idealized simple crack model. The trade‐off between Ω 0 and ƒ c may be more pronounced when using linearly sampled spectra, as higher frequency spectral bumps control the fits. As expected, methods that assumed a Boatwright model tended to have lower Ω 0 and somewhat higher ƒ c compared to those assuming a Brune model, although resulting Δσ estimates are similar. Finally, when compared to the overall validation study results, the fitting approach alone may account for between 5% and 90% (25% on average) of the total variability in spectral Δσ⁠.

58 GEOSCIENCES↗

NASA NTRS2010

Human Health and Performance Considerations for Exploration of Near-Earth Asteroids

This presentation will describe the human health and performance issues that are anticipated for the human exploration of near-Earth asteroids (NEA). Humans are considered a system in the design of any such deep-space exploration mission, and exploration of NEA presents unique challenges for the human system. Key factors that define the mission are those that are strongly affected by distance and duration. The most critical of these is deep-space radiation exposure without even the temporary shielding of a nearby large planetary body. The current space radiation permissible exposure limits (PEL) restrict mission duration to 3-10 months depending on age and gender of crewmembers and stage of the solar cycle. Factors that affect mission architecture include medical capability; countermeasures for bone, muscle, and cardiovascular atrophy during continuous weightlessness; restricted food supplies; and limited habitable volume. The design of a habitat that can maintain the physical and psychological health of the crew and support mission operations with limited intervention from Earth will require an integrated research and development effort by NASA s Human Research Program, engineering, and human factors groups. Limited abort and return options for an NEA mission are anticipated to have important effects on crew psychology as well as influence medical supplies and training requirements of the crew. Other important factors are those related to isolation, confinement, communication delays, autonomous operations, task design, small crew size, and even the unchanging view outside the windows for most of the mission. Geological properties of the NEA will influence design of sample handling and containment, and extravehicular activity capabilities including suit ports and tools. A robotic precursor mission that collects basic information on NEA surface properties would reduce uncertainty about these aspects of the mission as well as aid in design of mission architecture and exploration tasks.

Kundrot, Craig↗

NASA NTRS1984

Geophysical interpretation of satellite laser ranging measurements of crustal movement in California

As determined by satellite laser ranging the rate of contraction of a 900 kilometer baseline between sites located near Quincy in northern California and San Diego in southern California is about 61 to 65 mm/yr with a formal uncertainty of about 10 mm/yr. The measured changes in baseline length are a manifestation of the relative motion between the North America and Pacific tectonic plates. This long baseline result is compared to measurements made by more conventional means on shorter baselines. Additional information based on seismicity, geology, and theoretical modelling is also analyzed. Deformation lying within a few tens of kilometers about the major faults in southern California accounts for most, but not all of the observed motion. Further motion is attributable to a broader scale deformation in southern California. Data suggesting crustal movements north of the Garlock fault, in and near the southern Sierra Nevada and local motion at an observatory are also critically reviewed. The best estimates of overall motion indicated by ground observations lie between 40 and 60 mm/yr. This lies within one or two standard deviations of that deduced by satellite ranging but the possibility of some unresolved deficit cannot be dismissed. The long time scale RM2 plate tectonic model of Minster and Jordan predicts a contraction between 47 and 53 mm/yr depending on the extension rate of the Basin and Range. Thus the ground based observations, satellite laser ranging (SLR) results, and RM2 rates differ at about the 10 mm/yr level and are consistent with one another within the data and model uncertainties.

Cohen, S. C.↗

DOE OSTI2026

Denudation, solute export, landscape evolution modeling, and geographic information system data for the East River watershed, Colorado, USA (2020-2024)

This data package contains geographic information system (GIS) layers and tabular datasets associated with the study of lithologic controls on denudation, solute export, carbon-scaling relationships, and transient landscape evolution in the East River watershed near Crested Butte, Colorado, USA. The package includes GIS layers used to produce the Figure 2 map, including drainage, hillshade, lithology, sample locations, and basin polygons, together with comma-separated value (CSV) tables and matching CSV data dictionaries. One group of tables reports sample-level and catchment-level information for river-sediment samples analyzed for in situ-produced cosmogenic beryllium-10 (10Be), including sample names, outlet elevations, geographic coordinates, upstream drainage area, rock-type classes, production-rate scaling scheme, analyzed nuclide, catchment-averaged denudation rates, and associated lower and upper analytical uncertainties. Sample and catchment attributes provide the basis for comparing denudation rates across intrusive, shale, sedimentary, and mixed-lithology settings. A second group of tables reports supporting information for landscape-evolution modeling and the mapped geologic framework of the study area. Included files list parameter values and definitions for the two-phase landscape-evolution simulations, summarize full-domain model erosion fluxes and topographic metrics for different simulation configurations, provide a fixed-area carbon-model scaling table, and summarize mapped geologic units within the East River study domain, including geologic code, formation name, lithologic description, mapped area, and lithologic class grouping. Model outputs and geologic summaries support interpretation of transient landscape behavior and its relation to the mapped distribution of shale, intrusive, sedimentary, and surficial units. A third group of tables reports hydrologic and hydrochemical information used to quantify dissolved export from the watershed. Included files provide site-level values for drainage area, mean annual solute export, standard error of annual export, area-normalized solute yield, and equivalent weathering rate for five East River monitoring sites, along with metadata describing the number, sampling cadence, and date range of discharge records and partial and full total dissolved solids observations used in the solute-yield analyses. The package also contains a supplementary daily ion-load time series with daily mean discharge, discharge observation counts, dissolved concentrations, and daily loads for calcium, magnesium, sodium, potassium, chloride, sulfate, nitrate, fluoride, dissolved silica, charge-balance bicarbonate, and total dissolved solids. The package contains GIS files, comma-separated value files (.csv), CSV data dictionaries, a file-level metadata table, a package-tree text file, and a readme text file.

10Be↗

DOE OSTI2021

Atmospheric River Response to Arctic Sea Ice Loss in the Polar Amplification Model Intercomparison Project

The atmospheric river (AR) response to Arctic sea ice loss in the Northern hemisphere winter is investigated using simulations from the Polar Amplification Model Intercomparison Project. Results have shown that the midlatitude responses are dominated by dynamic effects. Poleward of around 60°N, the dynamic and thermodynamic effects cancel each other, resulting in relatively small responses. The response uncertainty can be characterized by leading uncertainty modes, with the responses over the Pacific and Atlantic projecting onto the northeastward extension and equatorward shift mode, respectively. In addition, the responses seem to be mean state-dependent: under the same forcing, models with more poleward-located climatological ARs tend to show stronger equatorward shifts over the Atlantic; over the Pacific, models with more westward-located climatological AR core tend to show stronger northeastward extensions. Finally, these relationships highlight the importance of improving the AR climatology representation on reducing the response uncertainty to Arctic sea ice loss.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1990

Carbon-isotopic analysis of dissolved acetate

Heating of dried, acetate-containing solids together with oxalic acid dihydrate conveniently releases acetic acid for purification by gas chromatography. For determination of the carbon-isotopic composition of total acetate, the acetate-containing zone of the chromatographic effluent can be routed directly to a combustion furnace coupled to a vacuum system allowing recovery, purification, and packaging of CO2 for mass-spectrometric analysis. For analysis of methyl carbon, acetic acid can be cryogenically trapped from the chromatographic effluent, then transferred to a tube containing excess NaOH. The tube is evacuated, sealed, and heated to 500 degrees C to produce methane by pyrolysis of sodium acetate. Subsequent combustion of the methane allows determination of the 13C content at the methyl position in the parent acetate. With typical blanks, the standard deviation of single analyses is less than 0.4% for acetate samples larger than 5 micromoles. A full treatment of uncertainties is outlined.

Non-NASA Center↗

DOE OSTI2020

Criticality Consideration for Geologic Disposal of DOE-Owned Spent Nuclear Fuels - 20255

Decisions must be made relative to packaging spent nuclear fuel (SNF) for storage, transportation, and eventual disposal. This paper evaluates potential criticality-related loading and packaging considerations for Department of Energy (DOE)-owned SNF that could be impacted favorably or unfavorably by disposal uncertainties and proposes approaches to minimize such impacts. Using Al-Based SNF as an example and drawing on analyses performed in support of the Yucca Mountain Repository License Application (LA), this paper explores the role of natural barriers, engineered barriers, and waste form characteristics relative to reducing the potential for postclosure criticality. Several parameters important to postclosure criticality are considered in order to determine which are useful discriminators to support packaging decisions. Using these parameters, a qualitative assessment is performed relative to three options for managing DOE-owned SNFs. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Initial Efforts Organizing WPNCS SG-8: Preservation of Expert Knowledge and Judgement Applied to Criticality Benchmarks

The Working Party on Nuclear Criticality Safety (WPNCS) under the guidance of the Organization for Economic Co-operation and Development (OECD) Nuclear Energy Agency (NEA) has over 20 years of experience addressing concerns related to static and transient configurations encountered within the nuclear fuel cycle: fuel fabrication, transportation, reprocessing, storage, and geological disposal. One of the cornerstone activities of the WPNCS is the International Criticality Safety Benchmark Evaluation Project (ICSBEP), which was established to identify a comprehensive set of criticality benchmark data, evaluate the data, including quantification of overall uncertainties; compile the data into a standardized format, perform sample calculations utilizing modern nuclear data sets and codes utilized in nuclear criticality safety, and formally document the work into a single source of verified benchmark data. Annually, members of the ICSBEP Technical Review Group (TRG) contribute evaluated benchmark data that undergoes comprehensive technical review prior to publication in the ICSBEP Handbook. In the years since the ICSBEP was established, there has been much work to prepare benchmark data to support validation activities in nuclear criticality safety. The 2020 edition of the ICSBEP Handbook contains acceptable benchmark specifications for 5,053 critical, subcritical, or near-critical configurations in 582 benchmark evaluations. Modern benchmark development benefits from decades of experienced international participants, a well-established handbook format, supplementary guides to deal with uncertainty quantification, and a comprehensive review process based upon independent reviews from international experts. The ICSBEP Handbook also contains 838 configurations deemed unacceptable to support criticality safety efforts. They are recorded, with the reasoning for their rejection, to preserve the experimental data, prevent reevaluation of data that are incomplete or contain known errors, and/or to potentially allow future reevaluation of the experiment pending the identification of sufficient data to resolve identified inconsistencies and errors. Users of the ICSBEP Handbook today might notice that the rigor and quality of modern criticality safety benchmarks is much greater than those prepared within the initial decade of the project. Benchmarks with 1s uncertainties in k eff greater than 1% were traditionally rejected unless they were identified as unique experiment types that encompassed materials, fuels, or designs not available from other benchmark experiments. However, benchmarks developed using modern experimental techniques and practices typically have uncertainties on the order of a few tenths of a percent. There have been ongoing efforts to improve the overall quality of previously published benchmark evaluations. Seventy-eight evaluations, containing approximately 600 configurations, have been revised just within the past decade. An additional eleven benchmarks are under revision for updated release in the 2020 edition of the ICSBEP Handbook. If some of the historic benchmarks were resubmitted in their current form to the TRG today, they would be rejected due to lack of data, missing components in the uncertainty analysis, or incomplete benchmark model development. The use of historic criticality safety benchmarks that underestimate the total uncertainty, lack properly quantified biases, or provide inadequate benchmark specifications do not sufficiently support modern criticality safety and nuclear data efforts. Although the ICSBEP Handbook is recognized by regulating bodies to support criticality safety, users are required to justify their reasons to ignore historic benchmark data and include additional safety margins within their designs. Discussions were held at the WPNCS 23rd Annual Meeting in September 2019 regarding the aforementioned issues. The resultant decision was to establish Subgroup 8 (SG-8): Preservation of Expert Knowledge and Judgement Applied to Criticality Benchmarks. The current activities of SG-8 are discussed herein.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

DEEPEN Leapfrog Geodata Model Cleaned and Reformatted Exploration Datasets from Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. As part of the DEEPEN 3D play fairway analysis (PFA) conducted at Newberry Volcano for multiple play types (conventional hydrothermal, superhot EGS, and supercritical), existing geoscientific exploration datasets needed to be acquired, cleaned, reformatted, and assembled in Leapfrog Geothermal. This GDR submission includes all of the cleaned and reformatted (X (m), Y (m), elevation (m), processed data values) datasets used to build the Leapfrog Geodata model. Existing datasets were acquired from the GDR, from AltaRock, and from other sources. This yielded the following datasets: - Digital elevation model produced from LiDAR data by Ramsey and Bard, 2016 - MT surveys from 2006, 2011, 2014, and 2017 (including single inversions) - Gravity surveys from 2006, 2007, and 2011 (including single) - Earthquake catalogs from PNSN, LLNL, and the Newberry EGS Demonstration project - Seismic velocity model from Templeton et al., 2014 - The Frone, 2015 temperature model and a new one produced through extrapolating downhole temperature measurements and the SMU temperature at depth maps. Two versions of the new model are provided: 250 m spacing and 500 m spacing - EarthVision geologic model with alteration from Moser et al., 2016 - Well data from EGS well 55-29, deep geothermal wells, coreholes (GEO N-2 through 5) and several thermal gradient holes - "Newberry Well Data:" Location, simple lithology, directional survey data, and temperature data for the 34 wells and coreholes used in the Newberry PFA Although there are additional 2D datasets available in the area, such as aeromagnetic surveys, these were not included in the analysis. While it may be possible to project these datasets into three dimensions by assuming the surface measurements do not vary with depth, this method is associated with high uncertainty. Preexisting inversions of these data were unavailable, and inverting additional geophysical datasets is outside the scope of this project.

15 GEOTHERMAL ENERGY↗

DOE OSTI2022

Robust Carbon Dioxide Plume Imaging Using Joint Tomographic Inversion of Seismic Onset Time and Distributed Pressure and Temperature Measurements (Final Report)

We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.

02 PETROLEUM↗

DOE OSTI2022

Validation of earthquake ground-motion models in southern California, USA, using precariously balanced rocks

Accurate estimates of earthquake ground shaking rely on uncertain ground-motion models derived from limited instrumental recordings of historical earthquakes. A critical issue is that there is currently no method to empirically validate the resultant ground-motion estimates of these models at the timescale of rare, large earthquakes; this lack of validation causes great uncertainty in ground-motion estimates. Here, we address this issue and validate ground-motion estimates for southern California utilizing the unexceeded ground motions recorded by 20 precariously balanced rocks. We used cosmogenic 10Be exposure dating to model the age of the precariously balanced rocks, which ranged from ca. 1 ka to ca. 50 ka, and calculated their probability of toppling at different ground-motion levels. With this rock data, we then validated the earthquake ground motions estimated by the Uniform California Earthquake Rupture Forecast, Version 3 (UCERF3) seismic-source characterization and the Next Generation Attenuation (NGA)-West2 ground-motion models. We found that no ground-motion model estimated levels of earthquake ground shaking consistent with the observed continued existence of all 20 precariously balanced rocks. The ground-motion model I14 estimated ground-motion levels that were inconsistent with the most rocks; therefore, I14 was invalidated and removed. At a 2475 year mean return period, the removal of this invalid ground-motion model resulted in a 2−7% reduction in the mean and a 10−36% reduction in the 5th−95th fractile uncertainty of the ground-motion estimates. Our findings demonstrate the value of empirical data from precariously balanced rocks as a validation tool for removing invalid ground-motion models and, in turn, reducing the uncertainty in earthquake ground-motion estimates.

Geology↗

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↗