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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Integrated Modeling Report for CO 2 Storage with Enhanced Oil Recovery in Northern Michigan

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. MRCSP has been investigating various reservoir characterization and modeling technologies related to Carbon Capture, Utilization, and Storage (CCUS) in conjunction with carbon dioxide enhanced oil recovery (CO 2 -EOR) operations in multiple Silurian-age (Niagaran), oil-bearing, carbonate pinnacle reefs in northern Michigan, USA. This report provides a comprehensive discussion of reservoir modeling studies that were conducted for tracking oil production, forecasting CO 2 plume migration, and estimating associated storage in a number of reefs that were at different stages of their CO 2 -EOR life cycle. The modeling process for simulating oil production, CO 2 injection, and associated storage in these reefs entails two phases. The first phase, geologic framework modeling, integrates all pertinent geological and geophysical data (from logs, cores and seismic surveys) about reservoir structure, geometry, rock types, and property distributions (porosity, permeability, water saturation) into a 3-D distributed grid-based static earth model (SEM). The second phase, dynamic reservoir modeling, uses the SEM as a platform to simulate the movement of oil, gas, water, and CO 2 within the reservoir during primary hydrocarbon production, as well as during subsequent phases such as CO 2 -injection assisted EOR, plume migration, and associated storage. In addition, an assessment of coupled process effects is also carried out, where the impacts of geochemical and geo-mechanical processed induced by CO 2 injection are studied.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2015

Predictive Modeling of Terrestrial Radiation Exposure from Geologic Materials

Aerial gamma ray surveys are important for those working in nuclear security and industry for determining locations of both anthropogenic radiological sources and natural occurrences of radionuclides. During an aerial gamma ray survey, a low flying aircraft, such as a helicopter, flies in a linear pattern across the survey area while measuring the gamma emissions with a sodium iodide (NaI) detector. Currently, if a gamma ray survey is being flown in an area, the only way to correct for geologic sources of gamma rays is to have flown the area previously. This is prohibitively expensive and would require complete national coverage. This project’s goal is to model the geologic contribution to radiological backgrounds using published geochemical data, GIS software, remote sensing, calculations, and modeling software. K, U and Th are the three major gamma emitters in geologic material. U and Th are assumed to be in secular equilibrium with their daughter isotopes. If K, U, and Th abundance values are known for a given geologic unit the expected gamma ray exposure rate can be calculated using the Grasty equation or by modeling software. Monte Carlo N-Particle Transport software (MCNP), developed by Los Alamos National Laboratory, is modeling software designed to simulate particles and their interactions with matter. Using this software, models have been created that represent various lithologies. These simulations randomly generate gamma ray photons at energy levels expected from natural radiologic sources. The photons take a random path through the simulated geologic media and deposit their energy at the end of their track. A series of nested spheres have been created and filled with simulated atmosphere to record energy deposition. Energies deposited are binned in the same manner as the NaI detectors used during an aerial survey. These models are used in place of the simplistic Grasty equation as they take into account absorption properties of the lithology which the simplistic equation ignores.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Pahute Mesa-Oasis Valley Hydrostratigraphic Framework Model for Corrective Action Units 101 and 102: Central and Western Pahute Mesa, Nye County, Nevada (Rev. 0)

A new, revised three-dimensional (3-D) hydrostratigraphic framework model (HFM) for the Central and Western Pahute Mesa Corrective Action Units 101 and 102 was completed in 2019. The initial (Phase I) model was completed in 2002 and was documented in A Hydrostratigraphic Model and Alternatives for the Groundwater Flow and Contaminant Transport Model of Corrective Action Units 101 and 102: Central and Western Pahute Mesa, Nye County, Nevada (BN, 2002). Subsequent flow and transport modeling revealed the need to collect additional data to reduce uncertainties. The Phase II data collection and characterization effort included extensive well drilling and testing activities, and rebuilding of the Phase I 3-D HFM, resulting in the Phase II HFM (NSTec, 2014). Further drilling (Well ER-20-12), data collection, and interpretive work (e.g., Thirsty Canyon Lineament [Wurtz and Day, 2018]) were incorporated into the Pahute Mesa-Oasis Valley (PM-OV) HFM. The PM-OV HFM incorporates the area of interest from the Phase II HFM and adds additional areas, primarily to the north and west with minor additions to the south and east. The new, expanded area matches the PM-OV groundwater basin as defined in Fenelon et al. (2016) and includes Pahute Mesa, a former nuclear testing area at the Nevada National Security Site; and Oasis Valley, a groundwater discharge area downgradient from contaminant source areas on Pahute Mesa. The data and information necessary to build the PM-OV HFM are provided in Navarro (2019). The PM-OV HFM will be used by hydrologic modelers who are tasked with developing a model to determine how contaminants are transported by groundwater flow in an area of complex geology. The model area is large (more than 6,250 square kilometers) and geologically complex, including Paleozoic- to Mesozoic-age structures (e.g., thrust faults, Thirsty Canyon Lineament), at least seven Tertiary-age calderas, several intrusive bodies, and many relatively recent basin-and-range normal faults. Investigators from the U.S. Geological Survey; National Security Technologies, LLC; National Laboratories; and U.S. Department of Energy contractors have organized the more than 300 volcanic units and other rocks in the model area into 77 hydrostratigraphic units (HSUs). The volcanic rocks are subdivided into 42 aquifers, 11 confining units, and 9 composite units (containing both aquifer and confining unit rocks). The underlying pre-Tertiary rocks were divided into 5 HSUs, comprising 2 aquifers and 3 confining units. The model also includes 1 alluvial aquifer unit, 8 intrusive confining units, and 1 granite confining unit. The model depicts the thickness, extent, and geometric relationships of these HSUs (“layers” in the model) along with all the major structural features that control them, including calderas and faults. The most substantive differences between the Phase II HFM and the PM-OVHFM include (1) addition of the new data and resulting revised interpretations from the additional area encompassed by the PM-OV groundwater basin; (2) incorporation of 16 new faults; (3) 12 boreholes from the northern extension; (4) 2 additional caldera-collapse collars; and (5) correlation of the simulated faulting and hydrostratigraphy in the eastern part of the model domain with the delineations represented in the neighboring Rainier Mesa-Shoshone Mountain HFM-Prime and Yucca Flat-Climax Mine HFM. All PM-OV HFM delineations were based on published and unpublished data from geophysical investigations, surface geologic maps, and drill-hole records. Phase II approached model uncertainty using the model-building software to test model configurations and hydrologic consequences using a range of values for parameters such as fault offset amounts and HSU elevations. The PM-OV HFM expands on this approach to encompass the entire PM-OV groundwater basin.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Model Package Report: Geoframework Model of the Hanford Site 100 Area

The purpose of the 100 Area Geoframework Model (GFM) is to provide a reasonable, consistent, and defensible three dimensional representation of the hydrostratigraphic units below the River Corridor at the Hanford Site to support contaminant fate and transport models. The GFM is a three dimensional representation of the subsurface geologic structure. From this, three dimensional geologic model-exported results, in the form of points or surfaces, are used as inputs to populate and assemble the various numerical model architectures. The objective of this report is to define the process used to produce a hydrostratigraphic model for the hydrostratigraphic units beneath the Hanford Site 100 Area. The GFM may support several other CH2M HILL Plateau Remediation Company project needs and objectives, including providing geologic information to support remedial investigations and actions and as a tool to present River Corridor geology and contaminant extents. The GFM was constructed based on information through 2019 available in the Integrated Data Management System, Hanford Environmental Information System, and the Hanford Site geologic contacts (GeoContacts) datasets. Revisions to the 100 Area GFM may occur to incorporate new data and information. Each version will be maintained in configuration control using date stamps for identifying supporting databases, figures, and interpretations. Supporting data include the final three dimensional geoframework surfaces, two dimensional structure and isopach maps, and all the geologic contact inputs and interpreted geologic log data. These final products provide a supporting set of geologic information that together allow the creation of the GFM. This information is managed, updated, and maintained via CH2M HILL Plateau Remediation Company under the Environmental Model Management Archive.

58 GEOSCIENCES↗
DOE OSTI2024

CO 2 rock physics modeling for reliable monitoring of geologic carbon storage

Monitoring, verification, and accounting (MVA) are crucial to ensure safe and long-term geologic carbon storage. Seismic monitoring is a key MVA technique that utilizes seismic data to infer elastic properties of CO 2 -saturated rocks. Reliable accounting of CO 2 in subsurface storage reservoirs and potential leakage zones requires an accurate rock physics model. However, the widely used CO 2 rock physics model based on the conventional Biot-Gassmann equation can substantially underestimate the influence of CO 2 saturation on seismic waves, leading to inaccurate accounting. We develop an accurate CO 2 rock physics model by accounting for both effects of the stress dependence of seismic velocities in porous rocks and CO 2 weakening on the rock framework. We validate our CO 2 rock physics model using the Kimberlina-1.2 model (a previously proposed geologic carbon storage site in California) and create time-lapse elastic property models with our new rock physics method. We compare the results with those obtained using the conventional Biot-Gassmann equation. Our innovative approach produces larger changes in elastic properties than the Biot-Gassmann results. Using our CO 2 rock physics model can replicate shear-wave speed reductions observed in the laboratory. Our rock physics model enhances the accuracy of time-lapse elastic-wave modeling and enables reliable CO 2 accounting using seismic monitoring.

58 GEOSCIENCES↗
DOE OSTI2013

Evaluation of experimentally measured and model-calculated pH for rock-brine-CO2 systems under geologic CO2 sequestration conditions

pH is an essential parameter for understanding the geochemical reactions that occur in rock-brine-CO2 systems when CO2 is injected into deep geologic formations for long-term storage. Due to a lack of reliable experimental methods, most laboratory studies conducted under geological CO2 sequestration (GCS) conditions have relied on thermodynamic modeling to estimate pH. The accuracy of these model predictions is typically uncertain. In our previous work, we have developed a method for pH determination by in-situ spectrophotometry. In the present work, we expanded the applicable pH range for this method and measured the pH of several rock-brine-CO2 systems at GCS conditions for five rock samples collected from ongoing GCS demonstration projects. Experimental measurements were compared with pH values calculated using several geochemical modeling approaches. The effect of different thermodynamic databases on the accuracy of model prediction was evaluated. Results indicate that the accuracy of model calculations is rock-dependent. For rocks comprised of carbonate and sandstone, model results generally agreed well with experimentally measured pH; however, for basalt, significant differences were observed. These discrepancies may be due to the models’ failure to fully account for certain reaction occurring between the basalt minerals the CO2-saturated brine solutions.

pH, carbon sequestration↗
DOE OSTI2026

Mountain Basin Controls on the Snow-to-Streamflow Signal: An AIC-Weighted Multiple Linear Regression Framework

A regression-based analysis quantifies how basin characteristics modulate the snow-to-streamflow signal. First, we use the ERA5-Land reanalysis gridded product (European Centre for Medium Range Weather Forecasts reanalysis 5 -Land component) for 4,655 hydrologic unit code - 10 (HUC10) mountain basins across the western United States (US) for water years 1987–2024. Linear regressions are performed for peak snow water equivalent (SWE) and annual streamflow for each mountain basin. Models use ordinary least squares in Python’s statsmodels package. After which, an Akaike Information Criterion (AIC)–weighted ensemble multiple linear regression (MLR) framework with 47 watershed traits is used to predict the linear regression coefficient of determination (r-squared) defining the ability of peak SWE to predict annual streamflow across all mountain basin. Predictor sets are constrained to avoid multicollinearity by excluding models with variance inflation factors (VIF) greater than 5. Mountain basin traits included in the MLR include seasonal climate, topography, vegetation type and structure, and bedrock geology. Accepted models are considered if their AIC is within 2.0 of the model with the minimum AIC, or best model. To compare predictor influence across acceptable models, we computed standardized regression coefficients. To evaluate structural redundancy among models, we constructed binary inclusion vectors for each acceptable model, denoting whether a predictor was present (1) or absent (0). Core predictor variables are defined as occurring in at least 67% of the acceptable models. For this regional analysis, only one model was found acceptable, with higher snow-to-streamflow translation (higher r-squared) occurring in colder mountain basins with higher relative winter precipitation, more snow accumulation and a lower fraction of annual precipitation that falls in the spring and summer. The second component of the data package uses previously published, high-resolution output from an integrated hydrological model of the East River watershed using the U.S. Geological Survey Groundwater and Surface water Flow model (GSFLOW, doi:10.15485/1998576). East River MLR expands upon the approach described above to explore the response of five streamflow metrics—annual streamflow, runoff efficiency, 7-day minimum flow, low-flow duration, and non-perennial stream fraction to snow system indicators including peak SWE, snow-covered area, snow disappearance date, and the fraction of basin area characterized by low-to-no snow, as well as seasonal precipitation and temperature, and annual hydrologic variables representing soil moisture, evapotranspiration (ET), the partitioning of incoming precipitation to evapotranspiration (ET/P), groundwater storage, and groundwater inflow to streams. MLR was done on all water years (P0: 1987-2024) and for each period as determined in the split analysis using pooled regression techniques (P1: 1987-2011 and P2: 2012-2024) to evaluate shifting predictor variable emphasis on streamflow generation. Results indicate that since 2012, peak SWE has lost statistical strength in its prediction of annual streamflow and runoff efficiency, and the indirect influence of spring temperature has emerged as critically important. Low-flow metrics remain largely influenced by soil moisture, vegetation water use and groundwater inflows with summer precipitation becoming a direct influence on minimum summer flow. Together, these data and Python-based analysis tools provide a framework for identifying the key watershed characteristics that control how streamflow responds to snow from year to year. The package also helps quantify uncertainty in statistical models and assess how snow–streamflow relationships vary across regions and over time. This dataset contains comma-separated values files (.csv), text files (.txt), python code files (.py), figure files (.png), and shapefiles (.cpg, .dbf, .prj, .sbn, .sbx, .shp, .xml). Further details on file contents and MLR execution can be found in the readme file and the FLMD files. Work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Collaborative Seismic Earth Model: Generation 2

Geological interpretations, earthquake source inversions and ground motion modeling, among other applications, require models that jointly resolve crustal and mantle structure. With the second generation of the Collaborative Seismic Earth Model (CSEM2), we present a global multi-resolution tomographic Earth model that serves this purpose. The model evolves through successive regional- and global-scale refinements. While the first generation aggregated regional models, with this study, we ensure consistency between all individual submodels, resulting in a model that accurately explains wave propagation across scales. Recent regional tomographic models were incorporated, comprising continental-scale inversions for Asia and Africa, as well as regional inversions for the Western US, Central Andes, Iran, and Southeast Asia. Across all regional refinements, over 793,000 source-receiver pairs contributed. Moreover, the long-wavelength Earth model (LOWE) introduces large-scale structures outside of pre-existing local refinements. A full-waveform inversion for global anisotropic P-and S-wave speed structure over a total of 194 iterations with a minimum period of 50 s on a large data set of 1 hr of waveform data from 2,423 earthquakes and over 6 million source-receiver pairs ensures that regional updates in the crust and uppermost mantle translate into updates of deeper, global-scale structure. To test the performance of CSEM2, we evaluate waveform fits between observed and synthetic seismograms at 50 s for an independent data set on the global scale, and on the regional scale for lower periods. We accurately simulate waveforms within and across regional refinements, maintaining the original resolution of the submodels embedded in the global framework.

58 GEOSCIENCES↗
DOE OSTI2016

Co-Optimization of CO 2 -EOR and Storage Processes in Mature Oil Reservoirs

This article presents an optimization methodology for CO 2 enhanced oil recovery in partially depleted reservoirs. A field-scale compositional reservoir flow model was developed for assessing the performance history of an active CO 2 flood and for optimizing both oil production and CO 2 storage in the Farnsworth Unit (FWU), Ochiltree County, Texas. A geological framework model constructed from geophysical, geological, and engineering data acquired from the FWU was the basis for all reservoir simulations and the optimization method. An equation of state was calibrated with laboratory fluid analyses and subsequently used to predict the thermodynamic minimum miscible pressure (MMP). Initial history calibrations of primary, secondary and tertiary recovery were conducted as the basis for the study. After a good match was achieved, an optimization approach consisting of a proxy or surrogate model was constructed with a polynomial response surface method (PRSM). The PRSM utilized an objective function that maximized both oil recovery and CO 2 storage. Experimental design was used to link uncertain parameters to the objective function. Control variables considered in this study included: water alternating gas cycle and ratio, production rates and bottom-hole pressure of injectors and producers. Other key parameters considered in the modeling process were CO 2 purchase, gas recycle and addition of infill wells and/or patterns. The PRSM proxy model was ‘trained’ or calibrated with a series of training simulations. This involved an iterative process until the surrogate model reached a specific validation criterion. A sensitivity analysis was first conducted to ascertain which of these control variables to retain in the surrogate model. A genetic algorithm with a mixed-integer capability optimization approach was employed to determine the optimum developmental strategy to maximize both oil recovery and CO 2 storage. The proxy model reduced the computational cost significantly. The validation criteria of the reduced order model ensured accuracy in the dynamic modeling results. The prediction outcome suggested robustness and reliability of the genetic algorithm for optimizing both oil recovery and CO 2 storage. The reservoir modeling approach used in this study illustrates an improved approach to optimizing oil production and CO 2 storage within partially depleted oil reservoirs such as FWU. Lastly, this study may serve as a benchmark for potential CO 2 –EOR projects in the Anadarko basin and/or geologically similar basins throughout the world.

02 PETROLEUM↗
DOE OSTI2021

Agreement of Analytical and Simulation-Based Estimates of the Required Land Depth in Climate Models

Previous analytical and simulation-based analyses suggest that deeper land surface models are needed to realistically simulate the terrestrial thermal state in climate models, with implications for land-atmosphere interactions. Analytical approaches mainly focused on the subsurface propagation of harmonics such as the annual temperature signal, and a direct comparison with climate-change model output has been elusive. Here, this study addresses the propagation of a harmonic pulse fitted to represent the timescale and amplitude of anthropogenic warming. Its comparison to land model simulations with stepwise increased bottom boundary depth leads to an agreement between the simulation-based and analytical frameworks for long-term climate trends. Any depth increase gradually decreases the relative error in the subsurface thermodynamics, and a minimum depth of 170 m is recommended to simulate the ground climate adequately. The approach provides an accurate estimate of the required land-model depth for climate-change simulations and assesses the relative bias in insufficiently deep land models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1992

Sudbury project (University of Muenster-Ontario Geological Survey): Summary of results - an updated impact model

In 1984 the Ontario Geological Survey initiated a research project on the Sudbury structure (SS) in cooperation with the University of Muenster. The project included field mapping (1984-1989) and petrographic, chemical, and isotope analyses of the major stratigraphic units of the SS. Four diploma theses and four doctoral theses were performed during the project (1984-1992). Specific results of the various investigations are reported. Selected areas of the SS were mapped and sampled: Footwall rocks; Footwall breccia and parts of the sublayer and lower section of the Sudbury Igneous Complex (SIC); Onaping Formation and the upper section of the SIC; and Sudbury breccia and adjacent Footwall rocks along extended profiles up to 55 km from the SIC. All these stratigraphic units of the SS were studied in substantial detail by previous workers. The most important characteristic of the previous research is that it was based either on a volcanic model or on a mixed volcanic-impact model for the origin of the SS. The present project was clearly directed toward a test of the impact origin of the SS without invoking an endogenic component. In general, our results confirm the most widely accepted stratigraphic division of the SS. However, our interpretation of some of the major stratigraphic units is different from most views expressed. The stratigraphy of the SS and its new interpretation is given as a basis for discussion.

Avermann, M.↗