Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geologic model”

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 1,009 records · Page 56

Surface Modeling to Support Small-Body Spacecraft Exploration and Proximity Operations

In order to simulate physically plausible surfaces that represent geologically evolved surfaces, demonstrating demanding surface-relative guidance navigation and control (GN&C) actions, such surfaces must be made to mimic the geological processes themselves. A report describes how, using software and algorithms to model body surfaces as a series of digital terrain maps, a series of processes was put in place that evolve the surface from some assumed nominal starting condition. The physical processes modeled in this algorithmic technique include fractal regolith substrate texturing, fractally textured rocks (of empirically derived size and distribution power laws), cratering, and regolith migration under potential energy gradient. Starting with a global model that may be determined observationally or created ad hoc, the surface evolution is begun. First, material of some assumed strength is layered on the global model in a fractally random pattern. Then, rocks are distributed according to power laws measured on the Moon. Cratering then takes place in a temporal fashion, including modeling of ejecta blankets and taking into account the gravity of the object (which determines how much of the ejecta blanket falls back to the surface), and causing the observed phenomena of older craters being progressively buried by the ejecta of earlier impacts. Finally, regolith migration occurs which stratifies finer materials from coarser, as the fine material progressively migrates to regions of lower potential energy.

Riedel, Joseph E.↗

Feasibility of a Deep Direct-Use Geothermal System at the University of Illinois Urbana-Champaign

This study assesses the feasibility of using deep direct-use (DDU) geothermal energy in agricultural research facilities on the University of Illinois at Urbana-Champaign campus to exploit low-temperature sedimentary basins, such as the Illinois Basin. Subsurface components of the system include extraction and injection wells and downhole pumps. Surface equipment includes heat pumps/exchangers, and fluid transport and monitoring systems. Two geologic formations in the region exhibit a potential as sources for geothermal energy, based on pre initial temperatures and flow rates of fluids. The St. Peter and Mt. Simon Sandstones lie at depths of 634 and 1,280 m, respectively. Geocellular modeling is used to characterize the reservoirs. A St. Peter Sandstone model was made for an area south of the campus. Petrophysical and geothermal properties used are based on data from the closest wells penetrating the formations. Characterization of the Mt. Simon Sandstone is in progress and is not discussed here. Extraction and injection flows simulated with different wellbore configurations provide estimates of fluid flow out of and into the reservoir. The models are used to optimize flow rates, bottomhole pressure, and temperature of the produced fluid. Individual wellbore models simulate subsurface heat loss and gain, providing guidance on the optimal type and amount of insulation in the wellbore. Design of the surface facilities will address aspects of fluid delivery, heat exchange, capital operating costs, heat loss, and corrosion. Heat capacity and flow rates are assessed to estimate life-cycle costs and benefits, including the environmental benefits of reducing greenhouse gases and water use and increased energy efficiency. A preliminary analysis of surface configurations for the DDU system (including cascading applications) based on building heat loads is being conducted to identify multiple system designs that will maximize performance, energy efficiency, and cost recovery.

15 GEOTHERMAL ENERGY↗

Stratigraphy of the crater Copernicus

The stratigraphy of copernicus based on its olivine absorption bands is presented. Earth based spectral data are used to develop models that also employ cratering mechanics to devise theories for Copernican geomorphology. General geologic information, spectral information, upper and lower stratigraphic units and a chart for model comparison are included in the stratigraphic analysis.

Paquette, R.↗

International collaboration in disposal research: comparative modeling of coupled processes in the DECOVALEX project

Abstract. This presentation gives an overview of an international research collaboration for advancing the understanding and modeling of coupled thermo-hydro-mechanical-chemical (THMC) processes in geological systems. The creation of the international DECOVALEX project, now running for more than 25 years, was motivated by the recognition that prediction of these coupled effects is an essential part of the performance and safety assessment of geologic disposal systems for radioactive waste and spent nuclear fuel. DECOVALEX emphasizes joint analysis and comparative modeling of state-of-the-art field and laboratory experiments, across a range of host rock options and repository designs. Participating research teams are from radioactive waste management organizations, national research institutes, regulatory agencies, universities, and consulting groups, providing a wide range of perspectives and solutions to these complex problems. The presentation provides examples of the research contributions made collectively in past DECOVALEX tasks and also touches on the unique modeling challenges tackled in the ongoing project phase, referred to as DECOVALEX-2023. The current phase comprises 17 partner organizations, about 50 modeling teams, and 7 modeling tasks, which cover a broad portfolio from fundamental studies on gas migration to full-scale in situ heater experiments in different host rocks to performance assessment studies. Together, these examples illustrate that the insight and scientific knowledge gained within the DECOVALEX project would not have been possible if one group had studied these complex THMC modeling challenges alone rather than within a truly collaborative setting.

Birkholzer, Jens T.↗

Impacts of Different Operation Conditions and Geological Formation Characteristics on CO2 Sequestration in Citronelle Dome, Alabama

Major concerns of carbon dioxide (CO2) sequestration in subsurface formations are knowledge of the well injectivity and gas storage capacity of the formation, the CO2 pressure and saturation plume extensions during and after injection, and the risks associated with CO2 leakage and fault reactivation. Saline reservoirs are considered as one of the target formations for CO2 sequestration through structural, residual, dissolution, and mineral trapping mechanisms. The boundary condition of the saline reservoir dictates the pressure and saturation plume extension of the injected supercritical CO2 that could expand over large distances. This can lead to sources of risk, e.g., leakage and/or fault reactivation due to presence of wells, thief zones, and geological discontinuities. Therefore, there is a critical need to develop a model that describes how risk-related performance metrics (i.e., the CO2 saturation plume size, the pressure differential plume area, and the pressure differential at specific locations) vary as a function of the size of injection, time following injection, injection operations, and geologic environment. In this study, a systematic reservoir modeling studies of anthropogenic CO2 sequestration in Citronelle dome, Alabama, was performed where all relevant scenarios and conditions to address the questions of the saturation and pressure plume size in the area of review (AoR) and post-injection site care (PISC) are considered. The objective for this study was firstly to systematically simulate CO2 sequestration, i.e., saturation dynamics, and pressure behavior over a range of operational and geological conditions and to derive conclusions about the factors influencing saturation and pressure plume size, post-injection behavior, and the risk associated with them, by developing third-generation reduced order models (ROMs) for reservoir behavior. Finally, to assess the uncertainty associated with our studies, Latin Hypercube Sampling (LHS) together with an experimental design technique, i.e., Plackett–Burman design, was used. Application of Pareto charts and respond surfaces enabled us to determine the most important parameters impacting saturation and pressure plume sizes and to quantify the auto- and cross-correlation among different parameters in both history-matched and upscaled models.

58 GEOSCIENCES↗

Satellite Gravity Drilling the Earth

Analysis of satellite-measured gravity and topography can provide crust-to-core mass variation models for new insi@t on the geologic evolution of the Earth. The internal structure of the Earth is mostly constrained by seismic observations and geochemical considerations. We suggest that these constraints may be augmented by gravity drilling that interprets satellite altitude free-air gravity observations for boundary undulations of the internal density layers related to mass flow. The approach involves separating the free-air anomalies into terrain-correlated and -decorrelated components based on the correlation spectrum between the anomalies and the gravity effects of the terrain. The terrain-decorrelated gravity anomalies are largely devoid of the long wavelength interfering effects of the terrain gravity and thus provide enhanced constraints for modeling mass variations of the mantle and core. For the Earth, subcrustal interpretations of the terrain-decorrelated anomalies are constrained by radially stratified densities inferred from seismic observations. These anomalies, with frequencies that clearly decrease as the density contrasts deepen, facilitate mapping mass flow patterns related to the thermodynamic state and evolution of the Earth's interior.

vonFrese, R. R. B.↗

Investigation of representing hysteresis in macroscopic models of two-phase flow in porous media using intermediate scale experimental data

Incorporating hysteresis into models is important to accurately capture the two phase flow behavior when porous media systems undergo cycles of drainage and imbibition such as in the cases of injection and post-injection redistribution of CO 2 during geological CO 2 storage (GCS). In the traditional model of two-phase flow, existing constitutive models that parameterize the hysteresis associated with these processes are generally based on the empirical relationships. This manuscript presents development and testing of mathematical hysteretic capillary pressure—saturation—relative permeability models with the objective of more accurately representing the redistribution of the fluids after injection. The constitutive models are developed by relating macroscopic variables to basic physics of two-phase capillary displacements at pore-scale and void space distribution properties. The modeling approach with the developed constitutive models with and without hysteresis as input is tested against some intermediate-scale flow cell experiments to test the ability of the models to represent movement and capillary trapping of immiscible fluids under macroscopically homogeneous and heterogeneous conditions. The hysteretic two-phase flow model predicted the overall plume migration and distribution during and post injection reasonably well and represented the postinjection behavior of the plume more accurately than the nonhysteretic models. In conclusion, based on the results in this study, neglecting hysteresis in the constitutive models of the traditional two-phase flow theory can seriously overpredict or underpredict the injected fluid distribution during post-injection under both homogeneous and heterogeneous conditions, depending on the selected value of the residual saturation in the nonhysteretic models.

58 GEOSCIENCES↗

The Role of Capillary Hysteresis and Pore-Scale Heterogeneity in Limiting the Migration of Buoyant Immiscible Fluids in Porous Media

Understanding the main mechanisms affecting long-term migration and redistribution of injected CO 2 in geological carbon storage is needed for developing predictive models to assess environmental risks and designing monitoring schemes. Preparation of a postinjection site care plan is required for CO 2 injection wells, including monitoring of pressure changes and injected CO 2 plume. Knowledge gaps exist regarding assessment of postinjection monitoring timeframes for the CO 2 plume because the processes driving long-term CO 2 plume migration and trapping are not fully understood. In the postinjection stage of geological carbon storage, redistribution of CO 2 occurs mainly due to buoyancy and capillary forces. This work presents experimental and modeling studies to investigate processes contributing to postinjection plume distribution and stabilization. We conducted a flow cell experiments (0.5 m × 0.05 m × 0.01 m) with two immiscible fluid phases in a glass bead porous medium to study postinjection plume behavior. We employed a hysteretic macroscopic two-phase flow model to interpret the experimental results and to understand main processes leading to plume stabilization. Our findings show that capillary pressure hysteresis explains the experimentally observed plume shape and redistribution at early postinjection stages; however, the long-term plume migration and eventual plume stabilization can only be represented when in addition microscale heterogeneity is accounted for. Results also show that plume stabilization can be extremely slow and that the migration of the plume front can occur through multiple intermittent bursts over long times. Here, further studies are needed to understand implications of the results for more realistic porous media and large-scale storage reservoirs.

58 GEOSCIENCES↗

Play-fairway analysis for geothermal exploration: Examples from the Great Basin, western USA

Elevated permeability within fault systems provides pathways for circulation of geothermal fluids. Future geothermal development depends on precise and accurate location of such fluid flow pathways in order to both accurately assess geothermal resource potential and increase drilling success rates. The collocation of geologic characteristics that promote permeability in a given geothermal system define the geothermal ‘fairway’, the location(s) where upflow zones are probable and where exploration efforts including drilling should be focused. We define the geothermal fairway as the collocation of 1) fault zones that are ideally oriented for slip or dilation under ambient stress conditions, 2) areas with a high spatial density of fault intersections, and 3) lithologies capable of supporting dense interconnected fracture networks. Areas in which these characteristics are concomitant with both elevated temperature and fluids are probable upflow zones where economic-scale, sustainable temperatures and flow rates are most likely to occur. Employing a variety of surface and subsurface data sets, we test this ‘play-fairway’ exploration methodology on two Great Basin geothermal systems, the actively producing Brady’s geothermal system and a ‘greenfield’ geothermal prospect at Astor Pass, NV. These analyses, based on 3D structural and stratigraphic framework models, reveal subsurface characteristics about each system, well beyond the scope of standard exploration methods. At Brady’s, the geothermal fairways we define correlate well with successful production wells and pinpoint several drilling targets for maintaining or expanding production in the field. In addition, hot-dry wells within the Brady’s geothermal field lie outside our defined geothermal fairways. At Astor Pass, our play-fairway analysis provides for a data-based conceptual model of fluid flow within the geothermal system and indicates several targets for exploration drilling.

15 GEOTHERMAL ENERGY↗

Lunar and Planetary Science XXVIII

The present conference discusses such topics as density crossovers in lunar picrites, the geology of the Cassini impact basin, Mars, nanobacteria in carbonates, the properties of shocked aerogels, a chemical model of Comet Halley, lunar mascons, the impact evolution of icy regoliths, the geology of the Venera 8 landing site, the photogeologic mapping of northern Venus, HST observations of Mars, observational constraints on the rotational dynamics of Mars, and primordial magnetic field measurements from the moon. Also discussed are models of the S2 fluorescence spectra of comets, Martian crater ejecta, the heights of Venusian steep-sided domes, cloud-climate interactions on Venus, the Humorum basin geology from Clementine data, an early Amazonian lake in the Gale crater of Mars, nebular fractionations and Mn-Cr systematics, the Rock Chipper planetary surface sample collection, Mariner 10 stereo images of Mercury, remote and local stresses and Calderas on Mars, the electrostatic charging of saltating particles, SO2 detected on Callisto, the Mars Explorer Planetary Data System, an assessment of explosive venting on Europa, the sequential faulting history of the Mars Valles Marineris, a search for Martian sediments, the composition and internal structure of Europa, long-term and 'diurnal' tidal stresses on Europa, and episodic greenhouse climates on Mars.

Source record↗

Using Direct Numerical Simulation of Pore-Level Events to Improve Pore-Network Models for Prediction of Residual Trapping of CO2

Direct numerical simulation and pore-network modeling are common approaches to study the physics of two-phase flow through natural rocks. For assessment of the long-term performance of geological sequestration of CO 2 , it is important to model the full drainage-imbibition cycle to provide an accurate estimate of the trapped CO 2 . While direct numerical simulation using pore geometry from micro-CT rock images accurately models two-phase flow physics, it is computationally prohibitive for large rock volumes. On the other hand, pore-network modeling on networks extracted from micro-CT rock images is computationally efficient but utilizes simplified physics in idealized geometric pore elements. This study uses the lattice-Boltzmann method for direct numerical simulation of CO 2 -brine flow in idealized pore elements to develop a new set of pore-level flow models for the pore-body filling and snap-off events in pore-network modeling of imbibition. Lattice-Boltzmann simulations are conducted on typical idealized pore-network configurations, and the interface evolution and local capillary pressure are evaluated to develop modified equations of local threshold capillary pressure of pore elements as a function of shape factor and other geometrical parameters. The modified equations are then incorporated into a quasi-static pore-network flow solver. The modified model is applied on extracted pore-network of sandstone samples, and saturation of residual trapped CO 2 is computed for a drainage-imbibition cycle. The modified model yields different statistics of pore-level events compared with the original model; in particular, the occurrence of snap-off in pore-throats is reduced resulting in a more frontal displacement pattern along the main injection direction. Compared to the original model, the modified model is in closer agreement with the residual trapped CO 2 obtained from core flow experiments and direct numerical simulation.

Kohanpur, Amir H.↗

Conceptualization of the Predevelopment Groundwater Flow System and Transient Water-Level Responses in Yucca Flat, Nevada National Security Site, Nevada

Contaminants introduced into the subsurface of Yucca Flat, Nevada National Security Site, by underground nuclear testing are of concern to the U.S. Department of Energy and regulators responsible for protecting human health and safety. The potential for contaminant movement away from the underground test areas and into the accessible environment is greatest by groundwater transport. The primary hydrologic control on this transport is evaluated and examined through a set of contour maps developed to represent the hydraulic-head distribution within the two major aquifer systems underlying the area. Aquifers and confining units within these systems were identified and their extents delineated by merging and analyzing hydrostratigraphic framework models developed by other investigators from existing geologic information. Maps of the hydraulic-head distributions in the major aquifer systems were developed from a detailed evaluation and assessment of available water-level measurements. The maps, in conjunction with regional and detailed hydrogeologic cross sections, were used to conceptualize flow within and between aquifer systems.

58 GEOSCIENCES↗

SMART Deliverable 6.1.2a: Application of the ORION tool to the IBDP Carbon Storage Site

Forecasting and managing potential induced seismic activity is one of the challenges facing commercialscale geologic carbon sequestration (GCS), as well as other geologic energy extraction and byproduct disposal technologies. Historically, the process to develop robust, science-based forecasts of induced seismicity has required an integrated effort from experts in seismology, geomechanics, and reservoir engineering to manage data, develop and evaluate models of subsurface processes, and to calibrate and interpret the results from a range of models to understand site behavior relative to prescribed standards and in the context of uncertainty in geologic characterization data, forecasting models, and operational scenario uncertainty. The Operational Forecasting of Induced Seismicity (ORION) toolkit is an open-source, observation-based forecasting toolkit that is being co-developed by two U.S. DOE-funded initiatives: the National Risk Assessment Partnership (NRAP) and the Scienceinformed Machine Learning for Accelerating Real Time Decisions in Subsurface Applications (SMART) Initiative. ORION is designed to provide functionality to support decision making about seismic hazard analysis and risk management for GCS stakeholders ranging from the public to site operators to expert seismologists. The tool, which is written as open-source code in the Python programming language, is composed of a desktop graphical user interface (GUI) and an underlying forecasting engine. The forecasting engine uses available reservoir properties, well and fluid injection scenario details, and observed seismic catalog data as inputs to produce a set of temporal and spatio-temporal seismic forecasts.

58 GEOSCIENCES↗

An Overview of the Source Physics Experiment at the Nevada National Security Site (SPE-N)

Understanding the changes in seismic energy as it travels from the near field to the far field is the ultimate goal in monitoring for explosive events of interest. This requires a clear understanding of explosion phenomenology as it relates to seismic, infrasound, and acoustic signals. Although there has been much progress in modeling these phenomena, this has been primarily based in the empirical realm. As a result, the logical next step in advancing the seismic monitoring capability of the United States is to conduct field tests that can expand the predictive capability of the physics-based modeling currently under development. The Source Physics Experiment at the Nevada National Security Site (SPE-N) is the first step in this endeavor to link the empirically based with the physics-based modeling. This is a collaborative project between National Security Technologies (NSTec), Lawrence Livermore National Laboratory (LLNL), Los Alamos National Laboratory (LANL), Sandia National Laboratories (SNL), the Defense Threat Reduction Agency (DTRA), and the Air Force Technical Applications Center (AFTAC). The test series require both the simple and complex cases to fully characterize the problem, which is to understand the transition of seismic energy from the near field to the far field; to understand the development of S-waves in explosives sources; and how anisotropy controls seismic energy transmission and partitioning. The current series is being conducted in a granite body called the Climax Stock. This location was chosen for several reasons, including the fairly homogenous granite; the location of previous nuclear tests in the same rock body; and generally the geology has been well characterized. The simple geology series is planned for 7 shots using conventional explosives in the same shot hole surrounded by Continuous Reflectometry for Radius vs. Time Experiment (CORRTEX), Time of Arrival (TOA), Velocity of Detonation (VOD), down-hole accelerometers, surface accelerometers, infrasound, and a suite of seismic sensors of various frequency bands from the near field to the far field. This allows for the use of a single test bed in the simple geology case instead of multiple tests beds to obtain the same results. The shots are planned at various depths to obtain a Green’s function, scaled-depth of burial data, nominal depth of burial data and damage zone data. SPE1-N was conducted in May 2011 as a 220 lb (100 kg) TNT equivalent calibration shot at a depth of 180 ft (55 m). SPE2-N was conducted in October 2011 as a 2200 lb (1000 kg) TNT equivalent calibration shot at a depth of 150 ft (46 m). SPE3-N was conducted in July 2012 as a 2200 lb (1000 kg) TNT equivalent calibration shot at a depth of 150 ft (46 m) in the damaged zone. Over 400 data channels were recorded for each of these shots and data recovery was about 95% with high signal to noise ratio. Once the simple geology site data has been utilized, a new test bed will be developed in a complex geology site to test these physics based models. Ultimately, the results from this project will provide the next advances in the science of monitoring to enable a physics-based predicative capability.

Source record↗

The metabolic regimes of 356 rivers in the United States

A national-scale quantification of metabolic energy flow in streams and rivers can improve understanding of the temporal dynamics of in-stream activity, links between energy cycling and ecosystem services, and the effects of human activities on aquatic metabolism. The two dominant terms in aquatic metabolism, gross primary production (GPP) and aerobic respiration (ER), have recently become practical to estimate for many sites due to improved modeling approaches and the availability of requisite model inputs in public datasets. We assembled inputs from the U.S. Geological Survey and National Aeronautics and Space Administration for October 2007 to January 2017. We then ran models to estimate daily GPP, ER, and the gas exchange rate coefficient for 356 streams and rivers across the continental United States. We also gathered potential explanatory variables and spatial information for cross-referencing this dataset with other datasets of watershed characteristics. Furthermore, this dataset offers a first national assessment of many-day time series of metabolic rates for up to 9 years per site, with a total of 490,907 site-days of estimates.

54 ENVIRONMENTAL SCIENCES↗

Efficient Multi-Scale, Thermo-Hydro-Mechanical Numerical Model for Simulation of Long-Term Stability of Rooms for Nuclear Waste Disposal

This project demonstrated the feasibility of developing a numerical modeling tool that provides an efficient and unique approach to analyze the geomechanical performance of a deep geological repository (DGR) for high-level nuclear waste. The numerical modeling tool is physics-based and models the thermo-hydro-mechanical (THM) processes during the required time (up to one million years) for a DGR. The unique capability of the developed code is explicit representation of fracturing and a discrete fracture network (DFN) in the model that includes the entire repository layout for simulated time of one million years. Therefore, it can potentially address one of the challenges of analysis of a DGR, which is efficient numerical simulation of a multi-scale, physics-based models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Geologic application of thermal inertia imaging using HCMM data

During the July to September 1980 quarter the final tapes were received completing the order and preliminary processing was done. Thermal Inertia images for each of the three test sites, Death Valley and Pisgah Crater, California and Goldfield, Nevada were created using registered HCMM day/night pairs and the JPL model. A comprehensive study and analysis of the geologic application of all acquired HCMM data is in progress.

Paley, H. N.↗

Radiolysis Process Model

Assessing the performance of spent (used) nuclear fuel in geological repository requires quantification of time-dependent phenomena that may influence its behavior on a time-scale up to millions of years. A high-level waste repository environment will be a dynamic redox system because of the time-dependent generation of radiolytic oxidants and reductants and the corrosion of Fe-bearing canister materials. One major difference between used fuel and natural analogues, including unirradiated UO 2 , is the intense radiolytic field. The radiation emitted by used fuel can produce radiolysis products in the presence of water vapor or a thin-film of water (including OH• and H• radicals, O 2 -, eaq, H 2 O 2 , H 2 , and O 2 ) that may increase the waste form degradation rate and change radionuclide behavior. H 2 O 2 is the dominant oxidant for spent nuclear fuel in an O 2 depleted water environment, the most sensitive parameters have been identified with respect to predictions of a radiolysis model under typical conditions. As compared with the full model with about 100 reactions it was found that only 30-40 of the reactions are required to determine [H 2 O 2 ] to one part in 10–5 and to preserve most of the predictions for major species. This allows a systematic approach for model simplification and offers guidance in designing experiments for validation.

Hydrogen Peroxide↗