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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Multiscale Modeling of CO 2 Migration and Trapping in Fractured Reservoirs with Validation by Model Comparison and Real-Site Applications (Final Technical Report)

This report documents the accomplishments achieved during the project titled “Multiscale Modeling of CO 2 Migration and Trapping in Fractured Reservoirs with Validation by Model Comparison and Real-Site Applications” funded by the US Department of Energy, Office of Fossil Energy. The objectives of the project were to develop modeling capabilities for predicting CO 2 and brine migration in fractured reservoirs during geologic carbon storage and to investigate the feasibility of carbon storage in fractured reservoirs using the newly developed modeling capabilities. To achieve these objectives new mass transfer functions were developed to adapt existing dual-continuum approaches – a commonly used approach in hydrocarbon reservoir modeling – to the CO 2 -brine system. In the dual-continuum approach the undisturbed rock matrix and the fractures are modeled as two separate continua, which are coupled through mass exchange of CO 2 and brine between the two continua. Mass transfer functions represent the mass exchange between fractures and the rock matrix, without requiring highly resolved numerical grids of individual fractures, thus allowing large modeling domains relevant to geologic carbon storage questions (tens of meters to hundreds of kilometers). Mass transfer functions for the exchange of both free-phase CO 2 and brine between fractures and rock matrix were developed for both the injection phase and the post-injection phase. During the injection phase the fractures are filled with CO 2 , while the rock matrix is initially filled with brine. Therefore, a mass transfer function was developed that is based on vertical displacement of brine by CO 2 in the rock matrix due to gravitational forces. Once CO 2 injection ceases and CO 2 migrates away from the direct vicinity of the injection site, the fractures will begin to re-fill with brine, while the rock matrix will have high CO 2 saturations. For these conditions a mass transfer functions was developed that is based on spontaneous imbibition, as the pore spaces in the rock matrix are much smaller than in the fractures. For formations with low permeability or high capillary entry pressure of the rock matrix, free-phase CO 2 may not be able to enter the rock matrix. However, free-phase CO 2 at the interface between fractures and rock matrix can dissolve into brine resident in the rock matrix. Therefore, additional mass transfer functions were developed based on aqueous-phase diffusion of CO 2 through the rock matrix. These mass transfer functions are based on existing solutions for diffusion and were expanded to be applicable to a variety of block shapes relevant to dual-continuum modeling. Implementations of the newly developed mass transfer functions showed that they give accurate results when compared to highly resolved models. In addition to developing new transfer functions, a new vertically-integrated dual-continuum modeling framework was developed. This development takes advantage of fast vertical segregation of CO 2 and brine in the fractures due to the density difference between CO 2 and brine and the high permeability of the fractures. Therefore, a vertically-integrated vertical equilibrium model for the fractures is coupled with a vertically-distributed model for the rock matrix. The newly developed approach significantly reduces the computational effort, but is only accurate for fracture permeabilities of 100 mD or higher. Finally, the newly developed models were applied to real and hypothetical geologic carbon storage sites to investigate the feasibility of using fractured reservoirs for carbon storage. The general results were, that storage of CO 2 – both as free-phase and dissolved – is feasible, but that storage capacities are likely lower than in unfractured systems. While the high permeability of the fractures reduces the injection pressure (or permits higher injection rates), fast migration of CO 2 may limit the overall storage volume, because CO 2 may escape laterally with not enough of the rock matrix being swept for effective storage.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2022

Process Interactions Can Change Process Ranking in a Coupled Complex System Under Process Model and Parametric Uncertainty

For a complex hydrologic system with multiple processes and process interactions, global sensitivity analysis is often used to identify important or influential parameters for model development and improvement. The identification is complicated by process model uncertainty, when a system process can be represented by multiple process models. This study develops a new total-effect process sensitivity index to identify influential processes under model uncertainty. This is done by extending Sobol's total-effect parameter sensitivity index for one system model to total-effect process sensitivity index for multiple system models to account for uncertainty in process models and model parameters. The total-effect process sensitivity index includes not only the first-order process sensitivity index for measuring the importance of individual processes but also higher-order indices that account for process interactions. The total-effect process sensitivity index can identify an influential process that itself and its interactions with other processes influence a model output. Here, the total-effect process sensitivity index is applied to two numerical examples: (a) Sobol's G*-functions with analytical solutions of first-order and total-effect process sensitivity indices, and (b) groundwater flow models with interactions between recharge, geology, and snowmelt processes. The second evaluation shows that, due to second-order and higher-order process interactions, the first-order and total-effect process sensitivity indices give different process ranking. It is thus necessary to estimate both first-order and total-effect process sensitivity indices to appreciate the difference between the first-order impact of a process alone and the overall total-effect impact of the process itself and its interactions with other processes on a model output.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Robust Ceramic Coaxial Cable Down-Hole Sensors for Long-Term In Situ Monitoring of Geologic CO 2 Injection and Storage

Energy generation by use of fossil fuels produces large volumes of CO 2 and other greenhouse gases, whose accumulation in the atmosphere is widely seen as undesirable. CO 2 Capture followed by sequestration has been identified as the solution. Subsurface geologic formations offer a potential location for long-term storage of CO 2 because of their requisite size. Unfortunately, the inaccessibility and complexity of the subsurface, the wide range of scales of variability, and the coupled nonlinear processes, impose tremendous challenges to determine the transport and predict the fate of the stored CO 2 . Among the various monitoring approaches, in situ down-hole monitoring of the various state parameters provides critical and direct data points that can be used to validate the models, optimize the injection, detect leakage and track the CO 2 plume. However, down-hole sensors that can withstand the harsh conditions and operate over decades of the project lifecycle remain unavailable. Given that the widespread of carbon capture and storage will be the necessity and reality in the future, fundamental and applied research is required to address the significant challenges and technological gaps in lack of long-term reliable down-hole sensors This project focused on the development and demonstration of a novel, low-cost, distributed, robust ceramic coaxial cable sensor platform for in situ down-hole monitoring of geologic CO 2 injection and storage with high spatial and temporal resolutions. The coaxial cable Fabry-Perot interferometer (CCFPI) has been studied as a general sensor platform for in situ, long-term, measurement of temperature, pressure and strain, which are critical to CO 2 injection and storage. A novel signal processing scheme has been developed and demonstrated for dense multiplexing of the sensors for low-cost distributed sensing with high spatial resolution. The developed temperature, pressure and strain sensors have been extensively tested under laboratory conditions that are similar to the downhole CO 2 storage environment, showing excellent capability for in situ monitoring the various parameters that are important to model, optimize the injection, detect leakage and track the CO 2 plume. In addition, the interactions between the sensor datum and the geological models have been investigated in details for the purposes of model validation, guiding sensor installation/placement, enhancement of model prediction capability and optimization of the injection processes. This project has resulted in the successful development of new ceramic coaxial cable based sensor systems that can monitor directly the changes in pressure, temperature, and strain caused by increased reservoir pressure and reduced reservoir temperature due to the supercritical CO 2 injection. Integrated with geological models, the sensors and measurement data can improve the possibility to identify plume movement and leakage in the cap rock and wells with higher precision and more accuracy. The low cost, ease of deployment, small size and dense multiplexing features of the new sensing technology will allow a large number of sensors to be deployed to address the objective to demonstrate that 99% of the CO 2 remains in the injection zone.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2016

Numerical Study of Critical Role of Rock Heterogeneity in Hydraulic Fracture Propagation

Log and seismic data indicate that most shale formations have strong heterogeneity. Conventional analytical and semi-analytical fracture models are not enough to simulate the complex fracture propagation in these highly heterogeneous formation. Without considering the intrinsic heterogeneity, predicted morphology of hydraulic fracture may be biased and misleading in optimizing the completion strategy. In this paper, a fully coupling fluid flow and geomechanics hydraulic fracture simulator based on dual-lattice Discrete Element Method (DEM) is used to predict the hydraulic fracture propagation in heterogeneous reservoir. The heterogeneity of rock is simulated by assigning different material force constant and critical strain to different particles and is adjusted by conditioning to the measured data and observed geological features. Based on proposed model, the effects of heterogeneity at different scale on micromechanical behavior and induced macroscopic fractures are examined. From the numerical results, the microcrack will be more inclined to form at the grain weaker interface. The conventional simulator with homogeneous assumption is not applicable for highly heterogeneous shale formation.

Dual Lattice Discrete Element Method↗
DOE OSTI2017

Improved image of intrusive bodies at Newberry Volcano, Oregon, based on 3D gravity modelling

Beneath Newberry Volcano is one of the largest geothermal heat reservoirs in the western United States and it has been extensively studied for the last 40 years. Several magmatic intrusions have been recognized at depths between 2.5 and 8 km and some of them identified as suitable targets for enhanced geothermal energy and tested during two previous EGS campaigns. These subsurface structures have been intersected by three deep wells and imaged by various geophysical methods including seismic tomography and magnetotellurics. Although three high quality gravity surveys were completed between 2006 and 2010 as part of various projects, a complete synthesis and interpretation of the gravity data has not yet been performed. Regional gravity data also exist in the vicinity of the Newberry volcano and have been added to these surveys to constitute a dataset with a total of 1418 gravity measurements. When coupled with existing geologic and geophysical data and models, this new gravity dataset provides important constraints on the depth and contours of the magmatic bodies previously identified by other methods and thus greatly contributing to facilitate any future drilling and stimulation works. Using the initial structures discovered by seismic tomography, inversion of gravity data has been performed. Shape, density values and depths of various bodies were allowed to vary and three main bodies have been identified. Densities of the middle and lower intrusive bodies (~2.6-2.7 g/cm3) are consistent with rhyolite, basalt or granites. Modeled density of the near-surface caldera body match that of a low density tephra material and the density of the shallow ring structures contained in the upper kilometer correspond to that of welded tuff or low-density rhyolites. Modeled bodies are in reality a composite of thin layers; however, average densities of the modeled gravity bodies are in good agreement with the density log obtained in one well located on the western flank (well 55-29). Final gravity data residuals show that most of the observed gravity anomalies at the surface can be explained by the modeled gravity bodies and are consistent with other site characterization information.

Cascade Range↗
NASA NTRS1993

Yana ring structure, North-Eastern Siberia: A possible counterpart of coronae on Venus

An arch of Verkhoyansky Range in NE Siberia, Russia is considered as W and S parts of the rim of the large ring structure, referred hereafter as Yana Ring Structure (YRS). It has general appearance similar to coronae on Venus and could be considered as the terrestrial counter part of coronae. The structure is located between 63 and 70 deg N, 125 and 140 deg E. The outer diameter of YRS (between the rim feet) is 680-700 km. The rim crest diameter is 500-520 km. The inner diameters (between the inner feet of the rim) is 320-350 km. The rim of YRS consist of the following topographic features: W segment is Orulgan Range (N part of Verkhoyansky Range), SW and S segments consist of S part of Verkhoyansky Range, SE segment is Elginskoye Ploskogorye (Tableland), E segment is Chersky Range (its part W of Indigirka River), NE segment is Selennyakhsky Range, and NW segment is Kular Range. Inside this ring of mountains Yanskoye Ploskogorye (Tableland) is located. The circular outlines of YRS are marked with the large river valleys located outside: Lena River (lower of Aldan River mouth) is located W from YRS, Aldan River (lower of Amga River mouth) is from S, Indigirka River (between 64 and 67 deg N lat) is from E, and Selennyakh River (the part running N to S) is from NE. The ring patterns inside the mountaineous rim of YRS are marked with the tributaries of Yana River. The largest ring here, 300 km in diameter, consists of Bytangay (with Bellyakh) River at W, Nelgese and Derbeke Rivers at S/SE, and Tuostakh River at E. The typical altitudes of YRS mountain rim are 1000-2000 m. The maximal ones are: 2400 m in W part of the rim, 2300 m in S, 1600 m in SE, 2200-2500 m in E, and 1300-1500 m in N. The central part of YRS have typical altitudes of 400-600 m. There are lowlands with altitudes 50-200 m to the S, W, N, and NW outside YRS Oymyakonskoye Nagorye (Highland) is located SE of YRS at 500-1500 m altitudes. It should be recognized from this description YRS have a general appearance of circular highland with its central plateau located 1000-2000 m lower than surrounding mountain rim, but 300-400 m higher than plains of lowlands located outside the rim. Such topographic patterns are typical for coronae on Venus. The diameter of YRS is equal to the largest coronae on Venus. There are only 3 or 4 structures of such or larger size on the whole Venus. This calls to compare morphology and geologic evolution of YRS with morphology and models of coronae evolutions especially the largest ones, to search the similarities and/or differences in geologic histories of the structures on the Earth and Venus.

Burba, G. A.↗
DOE OSTI2016

Development of an Integrated Hydrological Model for Simulation of Surface Runoff and Stream Flow in Tims Branch Watershed, Savannah River Site, South Carolina - 16203

This research is part of continued efforts to develop an integrated transport model to predict the distribution of tin within the overland flow in Tims Branch Watershed (TBW) at Savannah River Site (SRS), South Carolina. Tin was introduced into TBW during the application of an innovative mercury remediation technology implemented by the U.S. Department of Energy's Office of Environmental Management, which involved the injection of stannous (tin) chloride into mercury contaminated groundwater. Understanding the fate and transport of tin and its compounds is of primary importance due to their potential impact on the environment. In this study, we developed an overland hydrology model (MIKE SHE) capable of simulating surface hydrology throughout the TBW. The modeling application used historical precipitation, groundwater levels, geological data, and river discharges that were retrieved from government databases and input to the model. The model was developed to simulate spatio-temporal distribution of flow discharges, flow duration, and water levels in the TBW. The preliminary simulation results indicated that the model was capable of predicting flow depth and velocity within the study area during extreme climate events. The developed hydrology model will later on be coupled with transport model to simulate tin distribution within Tims Branch stream during seasonal storm events. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2020

The formation of Enceladus' Tiger Stripe Fractures from eccentricity tides

Enceladus has a young, tectonically active south polar region, which is erupting material from a prominent set of fractures called Tiger Stripes. No comparable activity is observed at the north pole, which is heavily cratered with limited tectonism. Given the many lines of evidence supporting a global ocean under Enceladus' icy shell, the reason for the dichotomy in geologic activity is unclear. We model the formation of the Tiger Stripes as tidally-driven fractures and examine the magnitudes of tidal stresses with different ice shell structures in order to explore whether and how tidal stress might explain Enceladus' distribution of tectonic activity. We find that eccentricity-driven tidal stresses would produce fractures of nearly identical orientations to the observed Tiger Stripe Fractures and that a 10-km difference in ice shell thickness between the north and south poles can result in substantially different tidal stress magnitudes, providing a natural explanation for the hemispheric dichotomy in tectonic activity on Enceladus. Finally, we synthesize these results with Enceladus' global geologic record to offer insight into the evolution of this enigmatic moon.

Alyssa RoseRhoden↗
NASA NTRS2002

Investigation of Error Patterns in Geographical Databases

The objective of the research conducted in this project is to develop a methodology to investigate the accuracy of Airport Safety Modeling Data (ASMD) using statistical, visualization, and Artificial Neural Network (ANN) techniques. Such a methodology can contribute to answering the following research questions: Over a representative sampling of ASMD databases, can statistical error analysis techniques be accurately learned and replicated by ANN modeling techniques? This representative ASMD sample should include numerous airports and a variety of terrain characterizations. Is it possible to identify and automate the recognition of patterns of error related to geographical features? Do such patterns of error relate to specific geographical features, such as elevation or terrain slope? Is it possible to combine the errors in small regions into an error prediction for a larger region? What are the data density reduction implications of this work? ASMD may be used as the source of terrain data for a synthetic visual system to be used in the cockpit of aircraft when visual reference to ground features is not possible during conditions of marginal weather or reduced visibility. In this research, United States Geologic Survey (USGS) digital elevation model (DEM) data has been selected as the benchmark. Artificial Neural Networks (ANNS) have been used and tested as alternate methods in place of the statistical methods in similar problems. They often perform better in pattern recognition, prediction and classification and categorization problems. Many studies show that when the data is complex and noisy, the accuracy of ANN models is generally higher than those of comparable traditional methods.

Dryer, David↗
DOE OSTI2023

An adjoint-based optimization method for jointly inverting heterogeneous material properties and fault slip from earthquake surface deformation data

SUMMARY Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson’s ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth’s structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M9 event, or distributed deformation along plate boundaries as constrained from InSAR.

Geochemistry & Geophysics↗
DOE OSTI2024

GDSA Repository Systems Analysis Investigations in FY 2024

The Disposal Research and Development (R&D) Program of the US Department of Energy (DOE) office of Nuclear Energy (NE-8) Spent Fuel and Waste Science and Technology (SFWST) Campaign is to conduct R&D on disposal of spent nuclear fuel (SNF) and high-level waste (HLW). The goal of the Geologic Disposal Safety Assessment (GDSA) within this project is to develop a disposal system modeling and analysis capability that supports the integrated modeling of coupled processes controlling disposal system performance of deep geologic repositories, including uncertainty. This report describes specific activities in the Fiscal Year (FY) 2024 associated with the GDSA Repository Systems Analysis (RSA) work package. The overall objective of the GDSA RSA work package is to develop generic deep geologic repository concepts and repository system performance models in crystalline, argillite, salt, and unsaturated alluvium potential host-rock environments, and to simulate and analyze these generic repository concepts and models using GDSA Framework toolkit, and other tools as needed.

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

Nevada national security site environmental remediation Large-Scale Multi-Well Aquifer Tests Used to Define Flow Systems - 15589

Over the past 5 years, 10 new groundwater wells have been drilled and hydraulically tested on Pahute Mesa Corrective Action Unit (CAU) under the Underground Test Area (UGTA) Activity at the Nevada National Security Site (NNSS). The new wells were drilled as part of a Phase II Corrective Action Investigation to gain a better understanding of the geology and hydrology of the system and to support development of a groundwater flow and transport model used for forecasting contaminant boundaries. Multi-well aquifer testing provides the most integrated assessment of hydraulic properties and connectivity within complex geologic systems. By providing hydraulic connectivity, these tests can be used to generate insight on potential contaminant flow paths and provide constraint on overall groundwater flow pathways. Drawdown in observation wells provides direct evidence of a hydraulic connection. During well drilling and testing activities, a total of 20 aquifer tests were conducted with as many as 36 individual completion zones. A signal-processing method developed by the U.S. Geological Survey was used to improve pumping drawdown detection. Using this approach, drawdown as far away as 3-4 kilometers has been detected. In addition to the typical analysis for hydraulic property values, the observed hydraulic connectivity are qualitatively and quantitatively used to inform the hydrogeologic framework model which is updated using other permissible geologic interpretations, when necessary, to incorporate the aquifer-testing results. These results will be used for developing and building confidence in the conceptual and numerical models for the Pahute Mesa CAU. (authors)

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

Heterogeneity, pore pressure, and injectate chemistry: Control measures for geologic carbon storage

Desirable outcomes for geologic carbon storage include maximizing storage efficiency, preserving injectivity, and avoiding unwanted consequences such as caprock or wellbore leakage or induced seismicity during and post injection. Here, to achieve these outcomes, three control measures are evident including pore pressure, injectate chemistry, and knowledge and prudent use of geologic heterogeneity. Field, experimental, and modeling examples are presented that demonstrate controllable GCS via these three measures. Observed changes in reservoir response accompanying CO 2 injection at the Cranfield (Mississippi, USA) site, along with lab testing, show potential for use of injectate chemistry as a means to alter fracture permeability (with concomitant improvements for sweep and storage efficiency). Further control of reservoir sweep attends brine extraction from reservoirs, with benefit for pressure control, mitigation of reservoir and wellbore damage, and water use. State-of-the-art validated models predict the extent of damage and deformation associated with pore pressure hazards in reservoirs, timing and location of networks of fractures, and development of localized leakage pathways. Experimentally validated geomechanics models show where wellbore failure is likely to occur during injection, and efficiency of repair methods. Use of heterogeneity as a control measure includes where best to inject, and where to avoid attempts at storage. Lastly, an example is use of waste zones or leaky seals to both reduce pore pressure hazards and enhance residual CO 2 trapping.

58 GEOSCIENCES↗
NASA NTRS1983

Two-dimensional Finite Element Modeling for Modeling Tectonic Stress and Strain

Techniques of finite element analysis in two dimensional plane strain were applied to problems of geophysics and tectonics. More specifically, the flexibility of the finite element method was employed to address problems involving geological complexity and fault interactions. The modeling of effective anisotropy in material elastic properties proved useful in describing the deformation of faulted crustal blocks. The applications of this modeling work to problems of actual tectonics in southern California was explored. Preliminary models show encouraging agreement with measured tectonic strain in this region, and modeling work was done to gain an understanding of the stress state in a locked fault region with future seismic potential.

Lyzenga, G. A.↗