Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geologic Modeling”

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

An Improved Model For Determining Salinity Recharge Time for Deep Borehole Disposal - 20405

In a previous publication (K. P. Travis, D. Burley and F. G. F. Gibb, WM2017 Conference, Phoenix Arizona paper no. 17480) we introduced a numerical model to estimate the time it would take for a pressure perturbation to subside following its creation from the construction of a deep geological borehole in water-saturated rock. The model was based on the notion of a point source of momentum - a solution of the time-dependent pressure-diffusion equation. Using the model, we estimated that it would take on the order of 10 k years for physical equilibrium to become re-established following the sinking of a 5 km borehole in granite. The significance of such a model is that it places an upper bound on the lifetime required of a borehole sealing system - a necessary input to a borehole post-closure safety case assessment. No engineered sealing system has ever been devised which is capable of retaining its sealing properties for 10 half-lives of long-lived radioisotopes in spent fuel or high-level waste. One of the key advantages of disposing of nuclear waste via Deep Borehole Disposal (DBD) is the natural sealing provided by density stratified groundwater. Upon sinking of a borehole, and subsequent filling with fresh water, brine or drilling mud, this natural barrier may be temporarily damaged. Over time, fresh brine from the far-field will flow towards or away from the hole (driven by a pressure gradient) and will eventually re-establish the original salinity gradient. It follows that the engineered seals need last only as long as the time required for this salinity gradient to reset itself. One of the limitations of our previous model was the use of a static boundary condition on the borehole wall. The model used a boundary pressure which varied quadratically with depth (arising from differences between the pressure of a column of fresh water in the borehole and that of a column of brine in the host rock). However as brine replaces fresh water in the borehole (driven by a pressure gradient), the boundary function must change with time. We now introduce an improved model which takes this time dependent boundary condition into account. The model also takes into account the time taken for a mixture of brine and fresh water to re-establish chemical equilibrium through the process of diffusion. The paper contains the mathematical details of our new iterative model as well as results showing the time taken to reach steady state, and concentration profiles for the components of the brine-filled borehole as a function of time together with a discussion on the implications of the results for developing a post-closure safety assessment for DBD. (authors)

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

A Solar-luminosity Model and Climate

Although the mechanisms of climatic change are not completely understood, the potential causes include changes in the Sun's luminosity. Solar activity in the form of sunspots, flares, proton events, and radiation fluctuations has displayed periodic tendencies. Two types of proxy climatic data that can be related to periodic solar activity are varved geologic formations and freshwater diatom deposits. A model for solar luminosity was developed by using the geometric progression of harmonic cycles that is evident in solar and geophysical data. The model assumes that variation in global energy input is a result of many periods of individual solar-luminosity variations. The 0.1-percent variation of the solar constant measured during the last sunspot cycle provided the basis for determining the amplitude of each luminosity cycle. Model output is a summation of the amplitudes of each cycle of a geometric progression of harmonic sine waves that are referenced to the 11-year average solar cycle. When the last eight cycles in Emiliani's oxygen-18 variations from deep-sea cores were standardized to the average length of glaciations during the Pleistocene (88,000 years), correlation coefficients with the model output ranged from 0.48 to 0.76. In order to calibrate the model to real time, model output was graphically compared to indirect records of glacial advances and retreats during the last 24,000 years and with sea-level rises during the Holocene. Carbon-14 production during the last millenium and elevations of the Great Salt Lake for the last 140 years demonstrate significant correlations with modeled luminosity. Major solar flares during the last 90 years match well with the time-calibrated model.

Perry, Charles A.↗
DOE OSTI2025

Coupled geomechanical investigation of depletion-induced fault reactivation

Fault reactivation during subsurface fluid production pose significant challenges to safe and sustainable resource extraction. Here, this study presents a three-dimensional coupled geomechanical framework to investigate the processes driving fault reactivation, capturing the interactions between reservoir dynamics and geomechanical responses. Verification against theoretical estimations based on linear poroelasticity confirms the model's capacity in representing reservoir background stress responses. However, the study reveals that relying solely on background stress states can underestimate or overestimate fault reactivation potential, emphasizing the importance of including localized stress perturbations such as differential compaction and stress redistribution. Applied to a fault (M1) inspired by the geological characteristics of the Groningen field, the model shows slip initiation at 2965 m depth with 16.0 MPa depletion, aligning with field observations where seismicity occurred at approximately 3 km depth after 15.8 MPa depletion. Parametric studies reveal: (1) inelastic reservoir compaction delays fault reactivation and mitigates fault slip by reducing stress concentration, (2) higher intermediate in-situ stress magnitudes decrease the Coulomb Failure Stress (CFS) increase rate and reduce fault slip, (3) larger fault offsets amplify shear stress near the offset zone, promoting earlier reactivation and longer rupture propagation, and (4) fault permeability significantly influences pressure diffusion, with low-permeability faults leading to sharper stress changes and earlier fault destabilization. These insights highlight the critical role of geological and mechanical parameters in fault reactivation and provide a predictive framework for mitigating induced seismicity risks.

Coupled geomechanical modeling↗
DOE OSTI2023

A phased workflow to define permit‐ready locations for large volume CO 2 injection and storage

Abstract To‐date, only two UIC Class VI permits have been issued by the US Environmental Protection Agency. We illustrate a four‐phase workflow to first identify regional storage resources and then down‐select sites to yield permit‐ready locations that can accept and store large volumes of CO 2 . Specific permit requirements should guide objectives and define deliverables of respective workflow phases. In the first phase we used available regional data and screened structure and injection zones to locate resources that match CO 2 volumes planned to be captured. Available data were also used to assess presence and depth of usable groundwater, the key resource being protected via permitting. We then used advanced, closed‐form, analytical solutions (EASiTool) to estimate CO 2 injectivity into each hydrologically connected injection compartment. In the second phase we acquired and conditioned additional wireline logs and leased available seismic datasets. We interpreted the depositional systems from wireline well‐log character and mapped sandbody geometry to interpolate injection and confining‐zone distribution. Using available data, we mapped faults and locations of freshwater and overpressure (or other capacity‐limiting geologic parameters) in more detail. In the third phase, we used the augmented geologic data to develop a static model for the selected area, extracted the areas of highest interest, and generated and ran dynamic (flow) models. In a fourth phase, we reduced major uncertainties identified in earlier phases. Our case study indicates that to complete preparation of a permit application requires (1) improved lithologic characterization information (thicknesses and horizontal and vertical connectivity) and (2) better definition of poorly defined local faults. © 2023 The Authors. Greenhouse Gases: Science and Technology published by Society of Chemical Industry and John Wiley & Sons Ltd.

58 GEOSCIENCES↗
DOE OSTI2024

Probabilistic Programming for Transportable Source Characterization and Uncertainty Quantification of the North Korean Nuclear Tests 2006–2017

Here, we introduce a transportable technique to determine the yield and depth of burial (DOB) from seismic source spectra of underground nuclear explosions. We demonstrate this technique on the six declared North Korean nuclear tests. This approach derives source spectra in absolute units from regional phase (Pg) amplitudes by correcting the observations for geometric spreading, attenuation, and site amplification. We couple the source spectra and explosion source models with a probabilistic programming framework that integrates deep learning techniques and Bayesian modeling. This approach permits the exchange of information across various data categories to quantify both the data and model uncertainty. This technique stands out as an innovative use of broad‐area propagation models, making it transportable across various geologic settings. This method proves to be effective in scenarios with diverse and/or limited observational data, even when the source depth is unknown. We present new independent estimates of absolute yield and DOB that are consistent with the prior assessments, underscoring the potential of this method in enhancing transportable nuclear explosion monitoring capabilities.

58 GEOSCIENCES↗
DOE OSTI2024

Utilizing 3D Mechanical Earth Models for Calibration and Validation in a Large-Scale Carbon Capture and Storage Project in North Dakota

Conference paper presented at 17th International Conference on Greenhouse Gas Control Technologies (GHGT-17), Calgary, Alberta, Canada, October 20–24, 2024. Three-dimensional (3D) mechanical earth models (MEMs) are pivotal in assessing geological sites for carbon dioxide (CO 2 ) storage and mitigating potential risks associated with storage and injection. The Energy & Environmental Research Center is exploring innovative carbon storage monitoring techniques at a CO 2 storage site near Beulah, North Dakota. These methods aim to provide lower-impact, faster feedback for commercial carbon capture and storage (CCS) projects. The research includes monitoring CO 2 injection using various approaches: 1) an automated, integrated modular monitoring station; 2) time-lapse electromagnetic (EM) field surveys; 3) drone-based surveillance; 4) time-lapse seismic techniques; and 5) advanced wellbore monitoring.

02 PETROLEUM↗
DOE OSTI2020

Technology Summary L-21044 Tunnel Sensor Hydrogeology Evaluations

Illicit tunneling along the U.S. Mexican border is a persistent national security problem. Several types of sensors could be to used to detect tunnels, but the hydrogeology of the subsurface makes some technologies ineffective. This effort used available borehole and surface geology data to develop 3D hydrogeologic models, which could then provide the input for optimizing subsurface sensor technologies to detect tunnels. Depending on the specific type of sensor, the following physical properties of the subsurface are important: longitudinal speed (velocity-p), shear speed (velocity-s), acoustic attenuation, electrical permittivity, electrical conductivity, density and elastic properties. LLNL developed a procedure which uses hydrogeological datasets to determine 3D hydrogeology cross sections. The physical properties of the rocks are generally not known precisely, thus ranges of properties would need to be considered. This cross-section data can then be utilized to enhance modeling studies to optimize sensor types/systems.

58 GEOSCIENCES↗
DOE OSTI2024

Machine Learning Applications in Analyzing the Role of Shale Barriers and Baffles for CO2 Storage

This study uses machine learning to analyze microseismic data from the Illinois Basin Decatur Project (IBDP) and quantify CO₂ plume extents. By leveraging well logs, microseismic records, and CO₂ injection metrics, the research predicts subsurface CO₂ plume dynamics. Findings show vertical clustering of microseismic events near the injection well, with CO₂ periodically breaching barriers due to buoyancy. K-Means clustering performed best, achieving the highest Silhouette Score and lowest Davies-Bouldin Index. This capability is crucial for real-time monitoring and management of CO₂ sequestration sites, validated against physical models and IBDP data, reinforcing CO₂ geological sequestration's viability and enhancing management tools.

Carr, Timothy↗
DOE OSTI2025

Optimizing Alabama’s CO 2 Storage in Shelby County (Project OASIS) - Deliverable Task 7: Risk Assessment

The “Optimizing Alabama’s CO 2 Storage in Shelby County: Project OASIS” CarbonSAFE Phase II Project seeks to build on regional data sets that demonstrate that the subsurface within Shelby County, Alabama has the potential to store commercial volumes of CO 2 safely, permanently, and economically. The primary target CO 2 Storage Complex is the deep Ketona Dolomite located within a 140 square mile area of the Valley and Ridge Region of Alabama. This deep saline reservoir is beneath a confining system encompassing at least 6,500 ft of shales and other low permeability sediments. Project OASIS drilled a deep stratigraphic test well to confirm the geological properties of the confining system and storage reservoir(s) within the Storage Complex. The geological data was incorporated into numerical models to establish the areal extent of the CO 2 plume and help design the storage site and its monitoring system. The project is managed by the Southern States Energy Board (SSEB), an interstate compact organization consisting of governors and state legislative leaders from sixteen southern states, Puerto Rico, and the U.S. Virgin Islands, as well as an appointee by the President of the United States. The organizational compact provides it with access to state government organizations and legislatures. The Board also maintains an Associate Members program comprised of energy resource companies, utilities, trade groups, academic R&D science and technology experts and energy consultants. Further, the SSEB staff is experienced in managing and coordinating complex energy and environmental programs, from research programs to full-scale design and demonstrations of new and innovative technologies. Project OASIS is a public-private partnership of six entities with multiple principal investigators (PIs). SSEB’s Lead PI and Co-PI are responsible for all aspects of project performance in accordance with the DOE-NETL Cooperative Agreement. SSEB has issued subgrants to Advanced Resources International, Inc., Alabama A&M University, Auburn University, Crescent Resource Innovation, and Oklahoma State University. Advanced Resources International, Inc., issued subgrants to Baker Hughes and Loudon Technical Services for field services

20 FOSSIL-FUELED POWER PLANTS↗
NASA NTRS1972

Thermal evolution of the moon.

The thermal history and current state of the lunar interior are investigated using constraints imposed by recent geological and physical data. Theoretical temperature models are computed taking into account different initial conditions, heat sources, differentiation and simulated convection. To account for the early formation of the lunar highlands, the time duration of magmatism and present-day temperatures estimated from lunar electrical conductivity profiles, it is necessary to restrict initial temperatures and abundances of radioactive elements. Successful models require that the outer half of the moon initially heated to melting temperatures, probably due to rapid accretion. Differentiation of radioactive heat sources toward the lunar surface occurred during the first 1.6 billion years. Temperatures in the outer 500 km are currently low, while the deep interior (radius less than 700 to 1000 km) is warmer than 1000 C, and is of primordial material.

Toksoz, M. N.↗
NASA NTRS1974

Solar luminosity variations and the climate of Mars

Attempts to resolve the solar neutrino flux problem have led to suggestions of large scale oscillations in the solar luminosity on a geological time scale. A simple climatological model of Mars indicates that its climate may be much more sensitive to luminosity changes than the earth's because of strong positive feedback mechanisms at work on Mars. Mariner-9 photographs of Mars show an abundance of large sinuous channels that point to an epoch of higher atmospheric pressures and abundant liquid water. Such an epoch could have been the result of large-scale, solar luminosity variations. However, our climatological model suggests that other less controversial mechanisms, such as obliquity or polar albedo changes, also could have led to such an epoch. As more becomes known about Mars, it may prove possible to formulate a history of Martian climate. By discovering effects that cannot be due to other mechanisms one may be able to form a chronology of solar luminosity variations to compare with data from the earth.

Toon, O. B.↗
NASA NTRS1979

Origin of Valles Marineris

The purpose of the paper is to propose an origin of Valles Marineris and contiguous features such as Noctis Labyrinthus, chaotic terrain, and the large channels of Mars. These are found to: (1) be continuous to each other, (2) have gradual transition between the morphologies, (3) be related by tectonic fracturing, and associated with features of possible magmatic origin. A model is presented in which the geological process is a magmatic process, and after initial uplift and fracturing, Valles Marineris, Noctis Labyrinthus and some chaotic terrain were produced by subsidence caused by lava withdrawal.

Schonfeld, E.↗
NASA NTRS1987

Morphologic contrasts between Nirgal and Auqakuh Valles, Mars: Evidence of different crustal properties

Photoclinometric measurements were made of sidewall slopes in Nirgal and Auqakuh Valles and these results were interpreted in terms of the geologic setting and a simple geomorphic model to provide insights into the physical properties of crustal materials in these areas. Nirgal was interpreted to be a runoff channel and Auqakuh to be a fretted channel. Geomorphologic arguments for the sapping origin of Nirgal and Auqakuh Valles were presented. The morphologies of the channels, however, differ greatly: the tributaries of Nirgal end abruptly in theater-headed canyons, whereas the heads of tributaries of Auqakuh shallow gradually. The plateau surface surrounding both channels appears to be covered by smooth materials, presumably lava flows; they are continuous and uneroded in the Nirgal area, but at Auqakuh they are largely eroded and several layers are exposed that total about 200 m in thickness. For Nirgal Valles, the measurements show that sidewalls in the ralatively shallow upper reaches of the channel have average slopes near 30 degrees and, in the lower reaches, sidewall slopes exceed 50 degrees. Auqakuh, on the other hand, has maximum sidewall slopes of 14 degrees and an approximate maximum depth of 1000 m. Faint, horizontal layering in portions of the lower reaches of Nirgal may indicate inhomogeneity in either composition or topography.

Mackinnon, David J.↗
NASA NTRS1989

Kinematics at the Intersection of the Garlock and Death Valley Fault Zones, California: Integration of TM Data and Field Studies

Kinematic relationships at the intersection of the southern Death Valley and Garlock fault zones were examined to identify and delineate the eastern structural boundary between the Mojave and the Basin and Range geologic terrains, and to construct a model for the evolution of this boundary through time. In order to accomplish this, satellite imagery was combined with field investigations to study six areas in the vicinity of the intersection, or possible extensions, of the fault zones. The information gathered from these areas allows the test of various hypotheses that were proposed to explain the interaction between the Death Valley and Garlock fault zones.

Verosub, Kenneth L.↗
NASA NTRS1991

Geology and structure of Beta Regio, Venus - Results from Arecibo radar imaging

Arecibo radar images of a portion of the equatorial region of Venus provide the first high resolution (1.5- to 2.0-km) synoptic coverage of Beta Regio. Within this area, tessera, a complex deformed terrain, is identified as a major geologic unit with the largest region corresponding to a plateau on the east flank of the highland. Three models are proposed to explain the origin and evolution of Beta Regio and are identified as Mantle Plume/Passive Crust, Mantle Plume/Active Crust, and Mantle Plume/Crustal Spreading. The Mantle Plume/Passive Crust model appears to be the most consistent with the geology in this region and suggests that a plume disrupts a preexisting region of tessera.

Senske, D. A.↗
NASA NTRS1993

Workshop on the Martian Northern Plains: Sedimentological, Periglacial, and Paleoclimatic Evolution

Papers that have been accepted for presentation at the Workshop on the Martian Northern Plains: Sedimentological, Periglacial, and Paleoclimatic Evolution, on 12-14 Aug. 1993 in Fairbanks, Alaska are included. Topics covered include: hydrological consequences of ponded water on Mars; morphological and morphometric studies of impact craters in the Northern Plains of Mars; a wet-geology and cold-climate Mars model: punctuation of a slow dynamics approach to equilibrium; the distribution of ground ice on Mars; and stratigraphy of the Martian Northern Plains.

Kargel, Jeffrey S.↗
NASA NTRS1988

The role of extraterrestrial phenomena in extinction

In the several years since the Alvarez report of anomalously high iridium concentrations at the Cretaceous-Tertiary boundary, evidence for the involvement of meteorite impacts in biological extinction has increased dramatically. Much more research will be needed, however, before meteorite impact is established as a general causal factor in extinction. Of ever greater long-term interest is the possibility that other extraterrestrial forces have had important influences on the evolution of life. To recognize the effects of such forces, it will be necessary to coordinate the research of astronomy and paleontology so that testable predictions can be formulated. It is possible that known, systematic changes in the Solar System or Galaxy have had effects on global biology and that these effects have been preserved in the paleontological record.

Non-NASA Center↗
NASA NTRS1988

The cryptoendolithic microbial environment in the Ross Desert of Antarctica: mathematical models of the thermal regime

Microbial activity in the Antarctic cryptoendolithic habitat is regulated primarily by temperature. Previous field studies have provided some information on the thermal regime in this habitat, but this type of information is limited by the remoteness of the site and the harsh climatic conditions. Therefore, a mathematical model of the endolithic thermal regime was constructed to augment the field data. This model enabled the parameters affecting the horizontal and altitudinal distribution of the community to be examined. The model predicts that colonization should be possible on surfaces with zenith angle less than 15 degrees. At greater zenith angles, colonization should be restricted to surfaces with azimuth angles less than 135 degrees or greater than 225 degrees. The upper elevational limit of the community should be less than 2,500 m. The thermal regime probably does not influence the zonation of the community within a rock.

Non-NASA Center↗