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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Upscaling Methods Applied to a Fine-Scale Reservoir Model

This study was conducted as part of the Southwest Regional Partnership on Carbon Sequestration (SWP) project to evaluate how upscaling fine-scale simulation models to coarse-scale simulation models impacted the results. The focus was on the Farnsworth Unit (FWU) and its Morrow' B' Sandstone reservoir, specifically the west half of the field. Due to data limitations and the geologic characteristics of the surrounding area, the upscaling was limited to the west half of the FWU rather than a broader basinscale model. The primary aim was to explore how upscaling impacts numerical simulation models, particularly regarding CO 2 -enhanced oil recovery (EOR) and storage capacity predictions. Upscaling was necessary to reduce computational demands when transitioning from high-resolution geological models to coarser grids, as large-scale simulations with finer grids can be computationally prohibitive. This study expands on previous work by the SWP to understand how additional upscaling, applied to already fine-scale numerical simulation models, affects reservoir performance simulations (Ampomah, Balch, & Grigg, 2015). This is key to understanding how loss of resolution can affect coarsescale model results that may be used for large sensitivity analyses, uncertainty quantifications, and training data for machine learning applications.

02 PETROLEUM↗
DOE OSTI2024

Dynamic modeling studies of basin-scale pressure interference and CO 2 plume evolution in multi-well geologic CO 2 storage

This research employs numerical modeling to analyze how CO 2 plumes and pressure fronts evolve when CO 2 is injected into a single storage formation from multiple projects located in proximity. Here, this analysis also evaluates the extent to which injection well spacing alleviates pressure buildup in the absence of active pressure management tactics. The simulation approach was based on a single, homogenous saline aquifer in which CO 2 injection occurs under a one-injector baseline case and several multi-well cases where well spacing varies. Analysis results show that the extent of pressure buildup is in the range of tens or a few hundreds of kilometers and contingent upon the defining pressure buildup demarcating the front edge. For the geological setting evaluated in this paper, our analysis suggests that without active basin pressure management strategies, commercial-scale projects would likely need to be sited far apart to avoid pressure interference from one another. Analysis results show the radius of CO 2 plume varies approximately from 2 to 3 km from injection wells (each injecting 1 Mt/year for 30 years) depending on cases and modeling parameters assumed. Given the pressure interference, this paper thus draws attention to the importance of greater coordination among storage operators and regulatory stakeholders. Because this analysis assumes a very specific geologic setting, this exploratory analysis bears further investigations across other geologic scenarios.

42 ENGINEERING↗
NASA NTRS1993

Collisional and dynamical history of Gaspra

Interpretation of the impact record on Gaspra requires understanding of the effects of collisions on a target body of Gaspra's size and shape, recognition of impact features that may have different morphologies from craters on larger planets, and models of the geological processes that erase and modify impact features. Crater counts on the 140 sq km of Gaspra imaged at highest resolution by the Galileo spacecraft show a steep size-frequency distribution (cumulative power-law index near -3.5) from the smallest resolvable size (150 m diameter) up through the large feature (1.5 km diameter crater) of familiar crater-like morphology. In addition, there appear to be as many as eight roughly circular concavities with diameters greater than 3 km visible on the asteroid. If we restrict our crater counts to features with traditionally recognized crater morphologies, these concavities would not be included. However, if we define craters to include any concave structures that may represent local or regional damage at an impact size, then the larger features on Gaspra are candidates for consideration. Acceptance of the multi-km features as craters has been cautious for several reasons. First, scaling laws (the physically plausible algorithms for extrapolating from experimental data) indicate that Gaspra could not have sustained such large-crater-forming impacts without being disrupted; second, aside from concavity, the larger structures have no other features (e.g. rims) that can be identified with known impact craters; and third, extrapolation of the power-law size distribution for smaller craters predicts no craters larger than 3 km over the entire surface. On the other hand, recent hydrocode modeling of impacts shows that for given impact (albeit into a sphere), the crater size is much larger than given by scaling laws. Gaspra-size bodies can sustain formation of up to 8-km craters without disruption. Besides allowing larger impact craters, this result doubles the lifetime since the last catastrophic fragmentation event up to one billion years. Events that create multi-km craters also globally damage the material structure, such that regolith is produced, whether or not Gaspra 'initially' had a regolith, contrary to other models in which initial regolith is required in order to allow current regolith. Because the globally destructive shock wave precedes basin formation, crater size is closer to the large size extrapolated from gravity-scaling rather than the strength-scaling that had earlier been assumed for such small bodies. This mechanism may also help explain the existence of Stickney on Phobos. Moreover, rejection of the large concavities as craters based on unfamiliar morphology would be premature, because (aside from Stickney) we have no other data on such large impact structures on such a small, irregular body. The eight candidate concavities cover an area greater than that counted for smaller craters, because they are most apparent where small craters cannot be seen: on low resolution images and at the limb on high resolution images. We estimate that there are at least two with diameter greater than 4 km per 140 sq km, which would have to be accounted for in any model that claims these are impact craters.

Greenberg, R.↗
DOE OSTI2013

3D Model of the McGinness Hills Geothermal Area

The McGinness Hills geothermal system lies in a ~8.5 km wide, north-northeast trending accommodation zone defined by east-dipping normal faults bounding the Toiyabe Range to the west and west-dipping normal faults bounding the Simpson Park Mountains to the east. Within this broad accommodation zone lies a fault step-over defined by north-northeast striking, west-dipping normal faults which step to the left at roughly the latitude of the McGinness Hills geothermal system. The McGinness Hills 3D model consists of 9 geologic units and 41 faults. The basal geologic units are metasediments of the Ordovician Valmy and Vininni Formations (undifferentiated in the model) which are intruded by Jurassic granitic rocks. Unconformably overlying is a ~100s m-thick section of Tertiary andesitic lava flows and four Oligocene-to-Miocene ash-flow tuffs: The Rattlesnake Canyon Tuff, tuff of Sutcliffe, the Cambell Creek Tuff and the Nine Hill tuff. Overlying are sequences of pre-to-syn-extensional Quaternary alluvium and post-extensional Quaternary alluvium. 10-15 degrees eastward dip of the Tertiary stratigraphy is controlled by the predominant west-dipping fault set. Geothermal production comes from two west dipping normal faults in the northern limb of the step over. Injection is into west dipping faults in the southern limb of the step over. Production and injection sites are in hydrologic communication, but at a deep level, as the northwest striking fault that links the southern and northern limbs of the step-over has no permeability.

15 GEOTHERMAL ENERGY↗
DOE OSTI2019

Broadband (0–5 Hz) Fully Deterministic 3D Ground‐Motion Simulations of a Magnitude 7.0 Hayward Fault Earthquake: Comparison with Empirical Ground‐Motion Models and 3D Path and Site Effects from Source Normalized Intensities

We report on high‐performance computing (HPC) fully deterministic simulation of ground motions for a moment magnitude (⁠Mw⁠) 7.0 scenario earthquake on the Hayward fault resolved to 5 Hz using the SW4 finite‐difference code. We computed motions obeying physics‐based 3D wave propagation at a regional scale with an Mw 7.0 kinematic rupture model generated following Graves and Pitarka (2016). Both plane‐layered (1D) and 3D Earth models were considered, with 3D subsurface material properties and topography interpolated from a model of the U.S. Geological Survey (USGS). The resulting ground‐motion intensities cover a broader frequency range than typically considered in regional‐scale simulations, including higher frequencies relevant for engineering analysis of structures. Median intensities for sites across the domain are within the reported between‐event uncertainties (⁠τ⁠) of ground‐motion models (GMMs) across spectral periods 0.2–10 s (frequencies 0.1–5 Hz). The within‐event standard deviation ϕ of ground‐motion intensity measurement residuals range 0.2–0.5 natural log units with values consistently larger for the 3D model. Source‐normalized ratios of intensities (3D/1D) reveal patterns of path and site effects that are correlated with known geologic structure. These results demonstrate that earthquake simulations with fully deterministic wave propagation in 3D Earth models on HPC platforms produce broadband ground motions with median and within‐event aleatory variability consistent with empirical models. Systematic intensity variations for the 3D model caused by path and site effects suggest that these epistemic effects can be estimated and removed to reduce variation in site‐specific hazard estimates. This study motivates future work to evaluate the validity of the USGS 3D model and investigate the development of path and site corrections by running more scenarios.

Rodgers, Arthur J.↗
DOE OSTI2023

DEEPEN 3D PFA Index Models for Exploration Datasets at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. As part of the development of the DEEPEN 3D play fairway analysis (PFA) methodology for magmatic plays (conventional hydrothermal, superhot EGS, and supercritical), index models needed to be developed to map values in geoscientific exploration datasets to favorability index values. This GDR submission includes those index models. Index models were created by binning values in exploration datasets into chunks based on their favorability, and then applying a number between 0 and 5 to each chunk, where 0 represents very unfavorable data values and 5 represents very favorable data values. To account for differences in how exploration methods are used to detect each play component, separate index models are produced for each exploration method for each component of each play type. Index models were created using histograms of the distributions of each exploration dataset in combination with literature and input from experts about what combinations of geophysical, geological, and geochemical signatures are considered favorable at Newberry. This is in attempt to create similar sized bins based on the current understanding of how different anomalies map to favorable areas for the different types of geothermal plays (i.e., conventional hydrothermal, superhot EGS, and supercritical). For example, an area of partial melt would likely appear as an area of low density, high conductivity, low vp, and high vp/vs. This means that these target anomalies would be given high (4 or 5) index values for the purpose of imaging the heat source. To account for differences in how exploration methods are used to detect each play component, separate index models are produced for each exploration method for each component of each play type. Index models were produced for the following datasets: - Geologic model - Alteration model - vp/vs - vp - vs - Temperature model - Seismicity (density*magnitude) - Density - Resistivity - Fault distance - Earthquake cutoff depth model

15 GEOTHERMAL ENERGY↗
DOE OSTI2018

Global perturbation of the marine calcium cycle during the Permian-Triassic transition

A negative shift in the calcium isotopic composition of marine carbonate rocks spanning the end-Permian extinction horizon in South China has been used to argue for an ocean acidification event coincident with mass extinction. This interpretation has proven controversial, both because the excursion has not been demonstrated across multiple, widely separated localities, and because modeling results of coupled carbon and calcium isotope records illustrate that calcium cycle imbalances alone cannot account for the full magnitude of the isotope excursion. Here, we further test potential controls on the Permian-Triassic calcium isotope record by measuring calcium isotope ratios from shallow-marine carbonate successions spanning the Permian-Triassic boundary in Turkey, Italy, and Oman. All measured sections display negative shifts in δ 44/40 Ca of up to 0.6‰. Consistency in the direction, magnitude, and timing of the calcium isotope excursion across these widely separated localities implies a primary and global δ 44/40 Ca signature. Based on the results of a coupled box model of the geological carbon and calcium cycles, we interpret the excursion to reflect a series of consequences arising from volcanic CO 2 release, including a temporary decrease in seawater δ 44/40 Ca due to short-lived ocean acidification and a more protracted increase in calcium isotope fractionation associated with a shift toward more primary aragonite in the sediment and, potentially, subsequently elevated carbonate saturation states caused by the persistence of elevated CO 2 delivery from volcanism. Locally, changing balances between aragonite and calcite production are sufficient to account for the calcium isotope excursions, but this effect alone does not explain the globally observed negative excursion in the δ 13 C values of carbonate sediments and organic matter as well. Only a carbon release event and related geochemical consequences are consistent both with calcium and carbon isotope data. Furthermore, the carbon release scenario can also account for oxygen isotope evidence for dramatic and protracted global warming as well as paleontological evidence for the preferential extinction of marine animals most susceptible to acidification, warming, and anoxia.

58 GEOSCIENCES↗
DOE OSTI2015

Anisotropic Elastic-Waveform Modeling for Fracture Characterization in EGS Reservoirs

Enhanced geothermal systems (EGS) contain newly created fractures in addition to possible existing fractures. Accurate characterization and monitoring of EGS reservoirs are crucial for optimal placement of new wells and effective extraction of geothermal heat. The fractured reservoirs behave as anisotropic media where seismic waves propagate with different velocities along different directions. In addition, the anisotropic properties of fluid-filled fracture zones could be different from those of dry fracture zones. We develop an optimized rotated staggered-grid elastic-wave finite-difference method for simulating seismic-wave propagation in heterogeneous, anisotropic media. Our new method uses a few extra grid points and optimized finite-difference coefficients based on the space-time dispersion relation, and reduce numerical dispersion of the conventional rotated staggered-grid finite-difference scheme. We validate our new method using synthetic vertical-seismic-profiling (VSP) data for an anisotropic geophysical model built with geologic features found at the Raft River EGS reservoir. This improved and optimized rotated staggered-grid finite-difference method provides an essential tool for analyzing VSP data, reverse-time migration, and elastic-waveform inversion in anisotropic, fractured reservoirs.

58 GEOSCIENCES↗
DOE OSTI2016

Amplitude Models for Discrimination and Yield Estimation

This seminar presentation describes amplitude models and yield estimations that look at the data in order to inform legislation. The following points were brought forth in the summary: global models that will predict three-component amplitudes (R-T-Z) were produced; Q models match regional geology; corrected source spectra can be used for discrimination and yield estimation; three-component data increase coverage and reduce scatter in source spectral estimates; three-component efforts must include distance-dependent effects; a community effort on instrument calibration is needed.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

Reduced diffusion and enhanced retention of multiple radionuclides from pore structure characterization of barrier materials for enhanced repository performance

Fluid flow and chemical transport in porous media are the macroscopic consequences of pore structure, which integrates geometry (e.g., pore size and surface area, pore-size distribution) and topology (e.g., pore connectivity). Low-permeability geological media whose pores are poorly interconnected will exhibit the characteristics of anomalous diffusion and sample size-dependent effective porosity, which will strongly impact long-term net diffusion and retention of radionuclides in geological repository settings involving different host rocks and barrier materials. A suite of innovative and complementary experimental approaches is utilized to study the microscopic pore structure and macroscopic fluid flow & chemical transport for a range of host rocks and barrier materials, in addition to standard clay minerals and reference rocks. With a particular focus on quantifying the presence and magnitude of “isolated” pores for a reduced effective porosity in low-permeability geomedia, the integrated methodologies for basic properties and pore structure characterization of these geomedia include X-ray diffraction, thin section petrography, grain size distribution, water immersion porosimetry after vacuum-pulling for full saturation, mercury intrusion porosimetry, nitrogen physisorption, scanning electron microscopy, X-ray computed tomography, and (ultra-)small angle neutron (X-ray) scattering. In addition, custom-designed gas diffusion, tracer recipe involving a range of anionic and cationic chemicals with subsequent analyses by laser ablation and inductively coupled plasma-mass spectrometry, along with batch sorption, column transport, and imbibition tests were conducted for coupled effects of pore structure and chemical retention/transport. From the perspectives of pore structure in conjunction with multiple and complementary approaches to examining a range of sample sizes under different observational scales, we find that the poor pore connectivity is prevalent in low-permeability media (mudstone and crystalline rock) that is related to geological processes (e.g., compaction, diagenesis and thermal maturation). For example, the deep and organic matter-rich mudstones have a much smaller effective porosity than the total porosity (as a result of poor pore connectivity) and associated diffusion coefficient, and the effective porosity & diffusion coefficients are also dependent upon the sample sizes used in the measurement. Similarly, most of the pore space in the shallow mudstone is also controlled by pore-throat diameters in the 5-50 nm range of intergranular pore types from its fine-grained nature, but with an overall good pore connectivity. However, the nm-sized pore space (physically pore-network architecture) and strong sorption capacities (chemical retention from clay minerals) of both shallow and deep mudstones lead to the synergistic retention of cationic radionuclides and their utilities as effective host rocks and barrier materials. Our unique approaches of studying how the micro-scale pore structure affect macro-scale fluid flow, diffusion & retention, and chemical transport produce improved mechanistic understanding, and realistic quantification, of diffusion and retention of typical radionuclides in a range of generic host rocks and barrier materials (clay/shale, salt, crystalline rock, and tuff), with the overall results leading to scientifically-based understanding of enhanced isolation (from both diffusion and retention) of radionuclides and improved confidence on the long-term performance of geological repository to store high-level radioactive wastes. In addition to the training of 25 undergraduates, graduates, and postdocs of UTA, the scientists (organizations) involved in performing this work (e.g., discussion, sample sharing, and operation of SANS and SAXS instruments) include Ed Matteo, Yifeng Wang, and Kristopher Kuhlman (Sandia National Laboratories), Jens Birkholzer, Liange Zheng, Tim Kneafsey, and Sharon Borglin (Lawrence Berkeley National Laboratory), Mavrik Zavarin (Lawrence Livermore National Laboratory), Yukio Tachi and Yuta Fukatsu (Japan Atomic Energy Agency), Mieke de Craen (Euridice, Belgium), Markus Bleuel (NIST), Wei-Ren Chen, Gergely Nagy, Changwoo Do, William Heller, Larry Anovitz, and Kenneth Littrell (ORNL), as well as Jan Illvsky, Ivan Kuzmenko, Ju-Sang Park and Jon Almers (ANL). Key deliverables include a total of 13 peer-reviewed journal articles (nine published and three under review), 23 presentations at scientific conferences (AAPG, AAPG Southwest Section, AGU, Asian Clay Conference, GSA, GSA South-Central Section, IHLRWM, InterPore, International Conference on Chemistry and Migration Behavior of Actinides and Fission Products in the Geosphere, International Conference on Coupled Processes in Fractured Geological Media: Observation, Modeling and Application), and academic institutions (UTA, New Mexico State University; University of Poitiers, France; University of Helsinki, Finland; Uppsala University, Sweden; Istanbul Technical University, Turkey) and other organizations (Andra, France; Posiva Oy, Finland).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

​Understanding a Stratigraphic Hydrothermal Resource – Geophysical Imaging at Steptoe Valley, Nevada

Sandia National Laboratories partnered with a multi-disciplinary group of subject matter experts to evaluate a stratigraphic geothermal resource in Steptoe Valley, Nevada using both established and novel geophysical imaging techniques. The stratigraphic reservoir in northern Steptoe Valley was previously discovered during oil and gas exploration. Subsequent studies, such as the Nevada Play Fairway Analysis, included data which further highlighted potential resource targets in the basin. Geophysical surveys, complemented with refined geologic mapping and geochemical sampling, were deployed to further characterize the resource. The resulting 3D geologic interpretation, conceptual model refinements, and reservoir simulations suggest that a >100MWe power-capable reservoir is likely economically accessible using conventional well placement and stimulation techniques in the Paleozoic carbonates of the deep/central basin of northern Steptoe Valley. Additional geophysical characterization and exploration drilling efforts are recommended to calibrate interpretation and determine where/how to potentially develop the northern Steptoe resource. The geophysical tools, interpretations, lessons learned, and public data generated by this study establish an exploration methodology to inform decisions for characterization and development of northern Steptoe Valley and other stratigraphic geothermal reservoirs in the western U.S.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Carbonated Brine Injection for a Pilot Site as a Low-Risk Geologic Carbon Storage Strategy: A Simulation Study

This technical report presents a preliminary study for CBI implementation in the Pena Creek pilot site in Dimmit County, Texas, targeting the Edwards Limestone Formation for SWD, in collaboration with NGL Energy Partners. The objective of this study is to find the conservative CO 2 molality for safe storage and potential risks for CBI at this pilot site. Geomodel was performed to get site-specific reservoir geology information. Reservoir modeling was performed to investigate dissolved CO 2 molality and pressure changes at relevant injection conditions. Wellbore corrosion and compatibility was discussed with corrosion modeling. Dissolved CO 2 long-term fate in the reservoir was examined with reactive transport modeling. This plan would set the scientific foundation as the initial step for future implementation in the field.

58 GEOSCIENCES↗
NASA NTRS1971

Genetic implications of the shapes of Martian and lunar craters.

Discussion of the problem of lunar crater genesis in the light of some considerations prompted by lunar and Martian crater similarities. By analogy with the moon, Martian craters seem both to vary in initial shape according to the energy of the impact that formed them and to have been modified subsequently by endogenic and surface processes. A proposed model for the geologic development of large Martian and lunar craters outlines a time-dependent sequence of events. Craters which have undergone rapid isostatic adjustment on the moon have distinctive morphologies and occur preferentially along mare basin-upland margins.

Pike, R. J.↗
NASA NTRS1975

Some aspects of core formation in Mercury

An evaluation of existing theories on the existence of the planet's metallic core is presented. Topics considered are: (1) magnetic fields; (2) surface geology; (3) cosmochemical models.

Solomon, S. C.↗
NASA NTRS1977

Earth resources-regional transfer activity contracts review

A regional transfer activity contracts review held by the Earth Resources Office was summarized. Contracts in the earth resources field primarily directed toward applications of satellite data and technology in solution of state and regional problems were reviewed. A summary of the progress of each contract was given in order to share experiences of researchers across a seven state region. The region included Missouri, Kentucky, Tennessee, Mississippi, Alabama, Georgia, and North Carolina. Research in several earth science disciplines included forestry, limnology, water resources, land use, geology, and mathematical modeling. The use of computers for establishment of information retrieval systems was also emphasized.

Bensko, J., Jr.↗
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.↗