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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Relaxation of the Martian Crustal Dichotomy Boundary in the Ismenius Region

The origin of the Martian crustal dichotomy remains a puzzle that when solved can provide an insight to the geological and geophysical evolution of Mars. In this study we model crustal relaxation in order to better constrain the original topographic shape, rheology, and temperature of the Martian crust. Our approach is to model the detailed geologic history of the Ismenius region of Mars, including slope, strain, and timing of faulting [1]. This region may contain the best preserved section of the dichotomy boundary as it is relatively unaffected by large impacts and erosion. So far the only study Martian crustal relaxation [2] suggests that the original topographic shape of the dichotomy is preserved. However, in this area strain from faulting implies at least some relaxation [1].

Guest, A.↗
DOE OSTI2019

Smart Proxy Modeling of SACROC CO2-EOR

Large CO2-enhanced oil recovery (EOR) projects usually contain an abundance of geological and good performance data. While this volume of data leads to robust models, it often results in difficult to manage, slow-running numerical flow models. To dramatically reduce the numerical run-times associated with the traditional simulation techniques, this work investigated the feasibility of using artificial intelligence and machine learning technologies to develop a smart proxy model of the Scurry Area Canyon Reef Operators Committee (SACROC) oilfield, located in the Permian Basin, TX, USA. Smart proxy models can be used to facilitate injection-production optimization for CO2-EOR projects. The use of a coupled grid-based, and well-based surrogate reservoir model (SRM) (also known as smart proxy modeling) was investigated as the base of the optimization. A fit-for-purpose coupled SRM, which executes in seconds, was built based on high-resolution numerical reservoir simulation models of the northern platform of the SACROC oilfield. This study is unique as it is the first application of coupled SRM at a large oilfield. The developed SRM was able to identify the dynamic reservoir properties (pressure, saturations, and component mole-fraction) at every grid-block, along with the production characteristics (pressure and rate) at each well. Recent attempts to use machine learning and pattern recognition to build proxy models have been simplistic, with limited predictive capabilities. The geological model used in this study is comprised of more than nine million grid blocks. The high correlation between the actual component and SRM, which can be visualized by mapping the properties, along with the fast footprint of the developed model demonstrate the successful application of this methodology.

Vida, Gholami↗
DOE OSTI2014

Modeling CO 2 -Brine-Rock Interaction Including Mercury and H 2 S Impurities in the Context of CO 2 Geologic Storage

This study uses modeling and simulation approaches to investigate the impacts on injectivity of trace amounts of mercury (Hg) in a carbon dioxide (CO 2 ) stream injected for geologic carbon sequestration in a sandstone reservoir at ~2.5 km depth. At the range of Hg concentrations expected (7-190 ppbV, or ~ 0.06-1.6 mg/std.m 3 CO 2 ), the total volumetric plugging that could occur due to complete condensation of Hg, or due to complete precipitation of Hg as cinnabar, results in a very small porosity change. In addition, Hg concentration much higher than the concentrations considered here would be required for Hg condensation to even occur. Concentration of aqueous Hg by water evaporation into CO 2 is also unlikely because the higher volatility of Hg relative to H 2 O at reservoir conditions prevents the Hg concentration from increasing in groundwater as dry CO 2 sweeps through, volatilizing both H 2 O and Hg. Using a model-derived aqueous solution to represent the formation water, batch reactive geochemical modeling show that the reaction of the formation water with the CO 2 -Hg mixture causes the pH to drop to about 4.7 and then become buffered near 5.2 upon reaction with the sediments, with a negligible net volume change from mineral dissolution and precipitation. Cinnabar (HgS(s)) is found to be thermodynamically stable as soon as the Hg-bearing CO 2 reacts with the formation water which contains small amounts of dissolved sulfide. Liquid mercury (Hg(l)) is not found to be thermodynamically stable at any point during the simulation. Two-dimensional radial reactive transport simulations of CO 2 injection at a rate of 14.8 kg/s into a 400 m-thick formation at isothermal conditions of 106°C and average pressure near 215 bar, with varying amounts of Hg and H 2 S trace gases, show generally that porosity changes only by about ±0.05% (absolute, i.e., new porosity = initial porosity ±0.0005) with Hg predicted to readily precipitate from the CO 2 as cinnabar in a zone mostly matching the single-phase CO 2 plume. The precipitation of minerals other than cinnabar, however, dominates the evolution of porosity. Main reactions include the replacement of primarily Fe-chlorite by siderite, of calcite by dolomite, and of K-feldspar by muscovite. Chalcedony is also predicted to precipitate from the dissolution of feldspars and quartz. Although the range of predicted porosity change is quite small, the amount of dissolution and precipitation predicted for these individual minerals is not negligible. These reactive transport simulations assume that Hg gas behaves ideally. To examine effects of non-ideality on these simulations, approximate calculations of the fugacity coefficient of Hg in CO 2 were made. Results suggest that Hg condensation could be significantly overestimated when assuming ideal gas behavior, making our simulation results conservative with respect to impacts on injectivity. The effect of pressure on Henry’s constant for Hg is estimated to yield Hg solubilities about 10% lower than when this effect is not considered, a change that is considered too small to affect the conclusions of this report. Although all results in this study are based on relatively mature data and modeling approaches, in the absence of experimental data and more detailed site-specific information, it is not possible to fully validate the results and conclusions.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1994

Venus: Vertical accretion of crust and depleted mantle and implications for geological history and processes

Models for the vertical accretion of a basaltic crust and depleted mantle layer on Venus over geologic time predict the eventual development of a net negatively buoyant depleted mantle layer, its foundering and its remixing with the underlying mantle. The consequences of the development of this layer, its loss, and the aftermath are investigated and compared to the geologic record of Venus revealed by Magellan. The young average age of the surface of Venus (several hundred million years), the formation of the heavily deformed tessera regions, the subsequent emplacement of widespread volcanic plains, the presently low rate of volcanic activity, and impact crater population that cannot be distinguished from a completely spatially random distribution, and the small number of impact craters embayed by volcanism, are all consistent with the development of a depleted mantle layer, its relatively rapid loss followed by large-scale volcanic flooding, and its subsequent reestablishment. We outline a 'catastrophic' tectonic resurfacing model in which the foundering of the depleted mantle layer several hundred million years ago caused globally extensive tectonic deformation and obliteration of the cratering record, accompanied by upwelling of warm fertile mantle and its pressure-release melting to produce extensive surface volcanism in the following period. Venus presently appears to be characterized by a relatively thick depleted mantle layer and lithosphere reestablished over the last several hundred million years following the previous instability event inferred to have produced the tessera terrain.

Head, James W.↗
DOE OSTI2015

Numerical simulations of vertical growth of hydraulic fractures and brine migration in geological formations above the Marcellus shale

One of the critical environmental questions about hydraulic fracturing in shales is the potential for contamination of ground and surface water. There are two specific concerns arising from hydraulic treatments: 1) whether hydraulic fractures extend upward through overlying strata to reach overlying aquifers containing drinking water, and 2) whether injected fluids push native fluids upward into these overlying aquifers. In this work, the extent of likely fracture growth through overlying layers during hydraulic treatment of the Marcellus shale was estimated using a hydraulic fracture model. A wide range of material and fluid flow properties in a multi-layered geologic model was considered. The model was based on conditions and characteristics applicable to the Marcellus shale in that part of the Appalachian basin within southwestern Pennsylvania. Predictions of vertical termination frequencies for hydraulic fractures were used in a multi-layer model of the strata and natural fractures for studying brine migration through the natural and induced fracture network. NFFLOW, the software for explicitly modeling flow within networks of fractures, was utilized to compute transient flow rates according to the schedule of injected fluid during hydraulic fracturing. To aid our analysis, the modeled sequence of geologic strata was capped with a fictitious unfractured, but moderately-permeable layer, which serves as a monitoring zone. The analysis assumes one well lateral was placed in the middle of the Marcellus shale with hydraulic fractures penetrating layers in the model. The newly-developed geomechanical module within NFFLOW was used to represent stress-sensitivity of the fractures. This allows the opening and closing of fracture apertures with changes in fluid pressures within fracture segments. Pressure increases in the formations overlying the Tully limestone, indicating fluid flow, was observed due to the hydraulic stimulation; and the impact of these increased pressures on brine migration towards the surface was considered.

03 NATURAL GAS↗
DOE OSTI2023

Temperature uncertainty modelling with proxy structural data as geostatistical constraints for well siting: an example applied to Granite Springs Valley, NV, USA

Utilizing existing temperature and structural geology information around Granite Springs Valley, Nevada, we build 3D stochastic temperature models with the aims of evaluating the 3D uncertainty of temperature and choosing between candidate exploration well locations. The data used to support the modelling are measured temperatures and structural proxies from 3D geologic modelling (distance to fault, distance to fault intersections and terminations, Coulomb stress change and dilation tendency), the latter considered ‘secondary’ data. Two stochastic geostatistical techniques are explored for incorporating the structural proxies: cosimulation and local varying mean. With both the cosimulation and local varying mean methods, many equally-likely temperature models (i.e. realizations) are produced, from which temperature probability profiles are calculated at candidate well locations. To aid in choosing between the candidate locations, two quantities summarize the temperature probabilities: V prior and entropy. V prior quantifies the likelihood for economic temperatures at each candidate location, whereas entropy identifies where new information has the most potential to reduce uncertainty. In general, the cosimulation realizations have smoother spatial structure, and extrapolate high temperatures at candidate locations that are located along the direction of the longest spatial correlation, which are down dip from existing temperature logs. The smooth realizations result in tight temperature probability profiles that are easier to interpret, but they have unrealistic temperature reversals in some locations because of the dipping ellipsoid shape created and that the cosimulation technique does not enforce a conductive geothermal gradient as a baseline (i.e. linearly increasing temperature with depth). The local varying mean results produce realizations with more realistic geothermal gradients, with temperatures increasing downward since a depth-temperature relationship is included. However, because they have much noisier spatial nature compared to cosimulation, it is harder to interpret the temperature probability profiles. The different local varying mean results allow the geologist to determine which proxy (e.g. dilation v. distance to fault termination) should be used given the specific geothermal system. In general, V prior from local varying mean results identify locations that are close to high values for the structural proxies: areas with higher probabilities for higher temperatures. The entropy results identify where uncertainty is greatest and therefore new drilling information could be most useful. Though these techniques provide useful information, even when applied to areas of sparse data, our comparison of these two techniques demonstrates the need for new geothermal geostatistics techniques that combine the advantages of these two methods and that are tailored to the spatial uncertainty issues inherent in geothermal exploration.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

An Integrated Multiscale Experimental-Numerical Analysis on Reconsolidation of Salt-Clay Mixture for Disposal of Heat-Generating Waste (Final NEUP Technical Report)

The overall purpose of this research is to improve understanding of THMC coupling effect on the reconsolidation of granular (or crushed) salt-clay mixture used for seal systems of shafts and drifts in salt repositories. This proposed work is partially motivated by the recent work on the Waste Isolation Pilot Plant (WIPP) that shows the promising sealing capability of clay-salt mixture compared to crushed salt. In particular, primary emphasis is to develop a fully integrated multiscale experiment-numerical study to determine and explain what leads to the superior sealing ability of the clay-salt mixture. These research activities are designed to seek further understanding of (1) why clay additives may enhance the fluid trapping and (2) whether this flow barrier effect may prevail under different combinations of temperature, confining pressure, deviatoric stress and other foreseeable environmental factors. If successful, this enhanced flow trapping ability of the seal provides significant improvement to the seal and repository performance and therefore make the repository safer in the long-term. The experiment component includes microstructural investigation and macroscopic tests on a reconsolidated salt-clay mixture. In the microstructural study, the goal is to (1) characterize microscopic distributions of distinct phases (e.g., clay, salt crystal boundaries, trapped brine, and pore) to examine the connectivity of the pore network inside the salt-clay mixture with different amounts of clay additive and moisture content and (2) analyze multiscale imaging data to reconstruct the polycrystalline microstructures for numerical simulations. Meanwhile, macroscopic tests are performed to analyze how clay alters the failure/creep mechanisms in the salt-clay mixture. Microscopic and macroscopic experimental observations will both be used to calibrate and validate a multiscale model that explicitly simulates the capillary and multiphase flow in the connected pores and the deformation due to the presence of intra-crystalline brine at the pore scale via a new polyhedral discrete element–lattice Boltzmann method (DEM-LBM) coupling model. The pore-scale simulations are homogenized via an upscaling procedure that converts pore-scale information (e.g. force exerted on grain boundary, sliding, pressure-solution) to continuum measures (e.g. Cauchy stress, Darcy’s flow) at each integration point in the macroscopic multiphase TMHC model. This multiscale scheme will allow coupling be- tween high-fidelity simulations of brine-salt-clay interaction and the macroscopic TMHC model. The multiscale model helps the understanding of how the trapped brine inclusion affects the pressure-solution mechanism with the presence of clay and moisture. This work brings new insight into the sealing capacity of salt-clay mixture under elevated temperature over a long period of time - a key to evaluating the potential of salt-clay mixture usage for salt repositories.

42 ENGINEERING↗
DOE OSTI2015

SECARB Presentation/Publication of Annual Results [Slides]

Contents include the following: 1) SECARB’s Anthropogenic Test: Phase III Update, March 9, 2010; 2) SECARB Anthropogenic Test Overview: Integrated CO 2 capture, transport, and storage, March 10, 2011; 3) SECARB Phase III Anthropogenic Test–Risk Management Through Detailed Geologic Characterization and Modeling, March 9-10, 2011; 4) SECARB Stakeholders Meeting Plant Barry CCS Demo, March 8, 2012; 5) SECARB Plant Barry to Citronelle 4" CO 2 Pipeline, March 8, 2012; 6) Geologic Characterization and Baseline Injection Modeling, March 7-8, 2012; 7) Permitting and Injection Status, March 7-8, 2012; 8) Southern Company/MHI Ltd. 500 TPD CCS Demonstration; 9) Advanced Monitoring Techniques, March 13, 2013; 10) SECARB Citronelle CO 2 Storage & MVA Status, March 13, 2013; 11) SECARB Plant Barry to Citronelle, March 13, 2013; 12) SECARB Awards Anthropogenic Test, March 4, 2014; 13) SECARB Anthropogenic Test Project Assessment, March 4-5, 2014; 14) Southern Company/MHI Ltd. 500 TPD CCS Demonstration; 15) Southern Company CCS R&D Plant Barry CCS Demo, March 12, 2015; and 16) SECARB Phase III Citronelle Project (Anthropogenic Test) in Alabama, March 12, 2015.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS1976

Basalt models for the Mars penetrator mission: Geology of the Amboy Lava Field, California

Amboy lava field (San Bernardino County, California) is a Holocene basalt flow selected as a test site for potential Mars Penetrators. A discussion is presented of (1) the general relations of basalt flow features and textures to styles of eruptions on earth, (2) the types of basalt flows likely to be encountered on Mars and the rationale for selection of the Amboy lava field as a test site, (3) the general geology of the Amboy lava field, and (4) detailed descriptions of the target sites at Amboy lava field.

Greeley, R.↗
DOE OSTI2024

Performance evaluation of the USGS velocity model for the San Francisco Bay Area

In this study, we evaluated the performance of the United States Geological Survey velocity model developed for the San Francisco Bay Area (SFBA), version 21.1. The evaluation was performed through high-resolution three-dimensional physics-based ground motion simulations of seven small-magnitude earthquakes (ranging from magnitude 3.8 to 4.4) that occurred on the eastern side of the San Francisco Bay. The simulations were performed in the frequency range from 0 to 5 Hz with a minimum shear-wave velocity of 250 m/s, which allowed the capture of wave propagation effects of the near-surface soft materials that characterize local basins. Based on the direct comparison of Fourier amplitude spectra between recorded and simulated ground motions for more than 250 stations, we found that the velocity model generally performs well in the frequency range of 0.2–5 Hz. The median value of the Fourier amplitude residuals was found to be near zero for all seven earthquakes. The slight over-prediction of 0.2 log-natural units at frequencies above 3 Hz in our simulations was attributed to the potentially inaccurate representation of the source radiation pattern by a double-couple point source model, and simple representation of shallow small-scale underground structural complexity in the velocity model. Maps of spectral amplitude differences between the simulated and recorded data were used to identify areas responsible for systematic ground motion over-predictions or under-predictions. For example, while some sub-domains over soft sediments show over-prediction patterns, the block east of the Hayward fault is prone to exhibit patterns of under-prediction. These maps can be used to guide future refinements of the SFBA velocity model. Since our simulation methodology allows for the decoupling of the source and wave propagation effects, the ground motion data generated by our simulations can also be used to quantify the epistemic uncertainty due to the velocity model, in empirically based ground motion estimates for the SFBA.

58 GEOSCIENCES↗
DOE OSTI2024

CarbonSAFE Phase II: Optimizing Alabama’s CO 2 Storage in Shelby County, Alabama (Project OASIS), Milestone M3 - Site Specific Drilling Report

The Phase II Storage Complex Feasibility project, entitled "Optimizing Alabama's CO 2 Storage in Shelby County, Alabama (Project OASIS)," is a part of the DOE/NETL's CarbonSAFE initiative. The Project is managed by the Southern States Energy Board (SSEB), and includes participation from Advanced Resources International, Inc. (ARI), Crescent Resource Innovation, Southern Company, Alabama A&M University, Auburn University, and Oklahoma State University. Project OASIS is working to establish the foundation for a commercial-scale geologic storage complex for CO 2 captured from Plant Gaston (home of the National Carbon Capture Center) and surrounding industrial sources of CO 2 located in Shelby County, Alabama. The Project objectives are: • Demonstrate that the subsurface saline formations at the storage complex can store commercial volumes of CO 2 safely and permanently. • Develop a comprehensive Community Benefits Plan. • Develop the infrastructure framework for a CO 2 storage hub. • Develop a rigorous risk registry and to conduct a comprehensive risk assessment. • Develop a monitoring plan. • Develop a comprehensive site characterization plan to support an Underground Injection Control Class VI Permit in a future Phase III program. • Evaluate the commercial viability of the project. Project OASIS is about 30 miles southeast of Birmingham, Alabama within a geologic province called the Valley and Ridge (Figure 1.1). The Valley and Ridge Province comprises a sequence of Paleozoic carbonate and clastic rocks that underwent structural deformation during the Alleghanian Orogeny. Storage prospects occur in relatively flat lying structural panels located between thrust faults. Available geologic studies related to hydrocarbon exploration suggest that Cambro-Ordovician carbonates and Cambrian clastic units offer multiple potential storage intervals, and that regional confining systems are present, such as the tectonically thickened Floyd-Parkwood Shale. The surface property is owned by a timber and land stewardship company, The Westervelt Company, Inc., who worked with the Project Team to select and prepare adequate sites for geologic assessment. The purpose of drilling the Westover Stratigraphic Test Well #2 was to collect geologic data to model the feasibility of commercial scale CO 2 injection and storage. This includes geological and geophysical evaluations, reservoir engineering analyses, and risk assessments.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2013

Preliminary Model Comparison Results From the Sim-SEQ Project Using TOUGH2, STOMP, Eclipse, and VESA Approach

Sim-SEQ is an international initiative on model comparison for geologic carbon sequestration (GCS), with an objective to understand and, if possible, quantify model uncertainties. Model comparison efforts in Sim-SEQ are initially limited to one specific field test site, hereafter referred to as the Sim-SEQ Study site (or S-3 site). Within Sim-SEQ, different modeling teams are developing conceptual models of CO2 injection at the S-3 site. One of the conceptual models, developed by the LBNL team, is based on TOUGH2/EOS7C. In this paper, we present some preliminary model predictions of the S-3 site using the TOUGH2/EOS7C simulator. We also compare the predictions of the TOUGH2 simulator with three other conceptual models, developed by three different organizations, of the S-3 site.

carbon sequestration↗
DOE OSTI2024

Continuous integration data-driven platform of industrial-scale subsurface storage for real-time analytics

This project helped address the growing need for efficient and scalable models to support geological carbon and energy storage, which are crucial for achieving net-zero emissions. Traditionally accurate high-fidelity numerical models have been used to simulate relevant storage processes under a handful of processes, however such models are computationally demanding, making uncertainty quantification impractical. Consequently, we first developed a machine learning framework, based on Graph Neural Operators (GNOs), to improving the accuracy of model predictions for a fixed computational budget. We then developed an Ensemble of Improved Neural Operators (ENO), which uses bagging and Monte Carlo dropout techniques, to further improve prediction accuracy. Lastly, we developed the way to explain progressive transfer learning methods to reduce the amount of training data and computational cost of training (i.e., reduce trainable parameters) when using our models for multiple storage sites. Our numerical investigation, which used real-world case studies, demonstrated that our framework can significantly improve the safety and efficiency of geological storage operations, with potential applications in other domains such as geothermal reservoirs and climate modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

The Importance of Interflow to Groundwater Recharge in a Snowmelt-Dominated, Alpine Watershed, Geophysical Research Letters: Modeling and Data Package

This data package contains hydrologic modeling input and output files for the East River (ER_PRMS.zip) and Copper Creek (CC_GSFLOW.zip) for water years 1987 to 2018. These data were used to generate simulated results published in Geophysical Research Letters (Carroll et al., 2019) describing the importance of interflow in generating groundwater recharge in the upper subalpine. The East River and Copper Creek are part of the Watershed Function Scientific Focus Area (WFSFA) located in the Upper Colorado River Basin, United States. The hydrologic model uses the U.S. Geological Survey coupled groundwater and surface water flow model (GSFLOW). Input and output files are askii format, while the model executable (gsflow.exe), source code, example problems and documentation describing model files is provided at the following U.S. Geological Survey website. https://www.usgs.gov/software/coupled-ground-water-and-surface-water-flow-model-gsflow. GSFLOW executables are compiled for either Microsoft Windows or Linux operating systems.

54 ENVIRONMENTAL SCIENCES↗