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 955 records · Page 53

Evaluation of Drilling Performance at The Geysers with Machine Learning Methods Using Geologic Data

A recent well, GDC-36, was drilled in The Geysers Geothermal Field served in a Department of Energy-industry to demonstrate improved drilling performance with polycrystalline diamond compact (PDC) bits. Both PDC and roller cone drill bits were used to drill this well. Key challenges encountered during drilling included lost circulation in the mud-drilled section, and bit damage interfacial severity in the deeper, air-drilled section. The objective of this study is to evaluate the drilling performance in relation to the local geological characteristics using machine learning methods. By applying K-clustering to the sonic log data, we were able to identify areas correlated with measured lost circulation. Also, the boundaries defined by clustering of the mineralogical and lithological data from the mud logs correlate well with interfacial severity during drilling. A random forest model was employed to build correlation between drilling data and rock strength. The confined compressive strength (CCS) of the rock in the training of the machine learning model was inferred from the dipole sonic log. The R-squared of the testing data is 0.78, and the RMSE (Root Mean Squared Error) is 0.06. The trained model was used to forecast rock strength for the section where sonic log data are not available. CCS could also be inferred from mud logs provided the relationship between mineralogy and rock strength is established through core testing data.

15 GEOTHERMAL ENERGY↗

Long-term sinking of nuclear waste canisters in salt formations by low-stress creep at high temperature

Abstract Rock salt has a self-sealing capacity, low permeability, and high thermal conductivity, making it a potential host for heat-generating nuclear waste. The feasibility of nuclear waste disposal within salt formations has been investigated mostly for small-sized canisters. Geologic disposal of larger-sized canisters originally designed for spent fuel storage and transportation has lately been examined as a cost-effective alternative. This raises questions about their long-term vertical movement due to their weight and high decay heat. Low-stress creep governs this movement; however, most salt constitutive models do not incorporate it. In this paper, the Norton and the WIPP creep models are compared with the Lux/Wolters/Lerche (LWL) model and a simpler model that combines linear and Norton creep laws (named combined creep model). The LWL and combined creep models consider pressure solution creep, though all incorporate dislocation creep. The models are first applied to creep tests under various stress levels. The LWL and the combined creep models results fit the experimental data well in both high and low stress ranges, whereas the Norton and WIPP models results only fit in higher-stress ranges. The different models are further applied for analyzing long-term canister movement. A sinking rate of $$-4.4\times 10^{-7}$$ - 4.4 × 10 - 7 mm/year was predicted using the Norton and WIPP models versus $$-2.1\times 10^{-2}$$ - 2.1 × 10 - 2 mm/year and $$-3.1\times 10^{-2}$$ - 3.1 × 10 - 2 mm/year using the LWL and the combined creep models, respectively. This comparative study confirms that creep models calibrated exclusively against high-deviatoric stress data might result in an inaccurate estimation of waste packages sinking rate in salt formations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Ultrafast radiographic imaging and tracking: An overview of instruments, methods, data, and applications

Ultrafast radiographic imaging and tracking (U-RadIT) use state-of-the-art ionizing particle and light sources to experimentally study sub-nanosecond transients or dynamic processes in physics, chemistry, biology, geology, materials science and other fields. These processes are fundamental to modern technologies and applications, such as nuclear fusion energy, advanced manufacturing, communication, and green transportation, which often involve one mole or more atoms and elementary particles, and thus are challenging to compute by using the first principles of quantum physics or other forward models. One of the central problems in U-RadIT is to optimize information yield through, e.g. high-luminosity X-ray and particle sources, efficient imaging and tracking detectors, novel methods to collect data, and large-bandwidth online and offline data processing, regulated by the underlying physics, statistics, and computing power. We review and highlight recent progress in: (a.) Detectors such as high-speed complementary metal-oxide semiconductor (CMOS) cameras, hybrid pixelated array detectors integrated with Timepix4 and other application-specific integrated circuits (ASICs), and digital photon detectors; (b.) U-RadIT modalities such as dynamic phase contrast imaging, dynamic diffractive imaging, and four-dimensional (4D) particle tracking; (c.) U-RadIT data and algorithms such as neural networks and machine learning, and (d.) Applications in ultrafast dynamic material science using XFELs, synchrotrons and laser-driven sources. Hardware-centric approaches to U-RadIT optimization are constrained by detector material properties, low signal-to-noise ratio, high cost and long development cycles of critical hardware components such as ASICs. Interpretation of experimental data, including comparisons with forward models, is frequently hindered by sparse measurements, model and measurement uncertainties, and noise. Alternatively, U-RadIT make increasing use of data science and machine learning algorithms, including experimental implementations of compressed sensing. Machine learning and artificial intelligence approaches, refined by physics and materials information, may also contribute significantly to data interpretation, uncertainty quantification and U-RadIT optimization.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Dislocation theory of steady and transient creep of crystalline solids: Predictions for olivine

In applications critical to the geological, materials, and engineering sciences, deformation occurs at strain rates too small to be accessible experimentally. Instead, extrapolations of empirical relationships are used, leading to epistemic uncertainties in predictions. To address these problems, we construct a theory of the fundamental processes affecting dislocations: storage and recovery. We then validate our theory for olivine deformation. This model explains the empirical relationships among strain rate, applied stress, and dislocation density in disparate laboratory regimes. It predicts the previously unexplained dependence of dislocation density on applied stress in olivine. The predictions of our model for Earth conditions differ from extrapolated empirical relationships. For example, it predicts rapid, transient deformation in the upper mantle, consistent with recent measurements of postseismic creep.

36 MATERIALS SCIENCE↗

Estimating perturbed stress from 3-D borehole displacements induced by fluid injection in fractured or faulted shales

SUMMARY Hydrofracturing stress measurements in fractured and anisotropic shales are notoriously difficult, because opening of existing geological features tends to prevent the creation of a pure hydraulic fracture perpendicular to the least compressive principal stress. Here we show how adding 3-D borehole-displacement measurements while conducting the hydraulic injection test helps to better constrain the principal stress orientations and magnitudes. We developed a 3-D fully coupled hydromechanical numerical model to analyse the displacement, fluid pressure and injection flow-rate data measured during an injection pressure-step-rate test conducted to activate a faulted borehole interval in the Mont Terri Opalinus Clay (Switzerland). We find that injected fluids can only penetrate the fault when it is at or above the Coulomb failure pressure. Borehole displacement orientations are sensitive to a ∼15° variation in the stress–tensor orientation and a 1 MPa stress magnitude variation. Although some dispersion occurs while rupture is propagating along the fault plane ∼4 m away from the borehole, the maximum density of displacement orientations consistently informs about the stress orientation. Thus, an extended injection step-rate approach coupled with an accurate in situ measurement of the borehole wall displacements can be used to better constrain the local stress field perturbations in fractured shales and in heterogeneous rock in general.

Guglielmi, Yves↗

Technoeconomic Analysis of Infrastructure Buildout Scenarios

The success of CCUS deployment in the southeast will depend on the optimized pipeline network from CO 2 point sources to geologic sinks for CO 2 storage. An optimal pipeline framework will reduce both environmental impacts and pipeline building cost. However, finding an optimal pipeline/transport infrastructure design is a non-trivial task and requires simultaneous usage of large volumes of data, computational resources, and a state-of-the-art simulator. Such a transport infrastructure model must consider point sources for CO 2 capture and associated volumes, sinks for CO 2 storage, and transportation from source to sink via pipeline networks. These considerations must be addressed simultaneously and systemically in an optimization process, in which a defined objective function (i.e., capital, variable CO 2 capture, transport, and storage cost function) is required to be minimized with the consideration of practical constraints (e.g., CO 2 flow through pipelines is less than the maximum capacity). Although optimization has been widely used in subsurface resources production and CO 2 storage, its application in large scale CCUS infrastructure design is rarely reported in the literature. SimCCS, developed by Los Alamos National Laboratory, is an open-source CCUS pipeline infrastructure design toolset that facilitates the optimization of pipeline infrastructure networks.

99 GENERAL AND MISCELLANEOUS↗

Revisiting ENSO and IOD Contributions to Australian Precipitation

Tropical modes of variability, such as El Niño–Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD), exert a strong influence on the interannual variability of Australian precipitation. Nevertheless, commonly used indices of ENSO and IOD variability display significant co-variability that prevents a robust quantification of the independent contribution of each mode to precipitation anomalies. This co-variability issue is often addressed by statistically removing ENSO or IOD variability from the precipitation field before calculating teleconnection patterns. However, by performing a suite of coupled and uncoupled modeling experiments in which either ENSO or IOD variability is physically removed, we show that ENSO-only-driven precipitation patterns computed by statistically removing the IOD influence significantly underestimate the impact of ENSO on Australian precipitation variability. Here, inspired by this, we propose a conceptual model that allows one to effectively separate the contribution of each mode to Australian precipitation variability.

54 ENVIRONMENTAL SCIENCES↗

Exploiting SMILEs and the CMIP5 Archive to Understand Arctic Climate Change Seasonality and Uncertainty

Arctic Amplification (AA) exhibits a distinct seasonal dependence; it is weakest in boreal summer and strongest in winter. Here, we analyze simulations from single-model initial-condition large ensembles and Coupled Model Intercomparison Project Phase 5 to decipher the seasonal evolution of Arctic climate change. Models agree that the annual maximum AA shifts from autumn into winter over the 21st century, accompanied by similar shifts in sea-ice loss and surface turbulent heat fluxes, whereas the maximum precipitation shifts only into late autumn. However, the exact seasonal timing and magnitude of these shifts are highly uncertain. Decomposing the uncertainty into model structural differences, emission scenarios, and internal variability reveals that model differences dominate the total uncertainty, which also undergo autumn-to-winter shifts. We also find that the scenario uncertainty is unimportant for projections of AA. These results highlight that understanding model differences is critical to reducing uncertainty in projected Arctic climate change.

54 ENVIRONMENTAL SCIENCES↗

The Roles of Tidal Marshes in the Estuarine Biochemical Processes: A Numerical Modeling Study

Observations suggest that the existence of tidal marsh can alter the oxygen and nutrient dynamics in adjacent water bodies, but assessing the impacts of large tidal marshes on an estuary is challenging. In this study, we use a modeling approach to investigate the roles of tidal marshes on the estuarine biochemical processes. The marsh model, which simulates the ecological functions of marshes at seasonal and annual time-scales, is embedded inside an unstructured-grid three-dimensional hydrodynamic and eutrophication model (SCHISM-ICM). This modeling system simulates the growth and metabolism of the tidal marshes and links biological processes to nutrient dynamics in the water column and sediment. This model dynamically simulates nutrient recycling and physical transport of the materials between marshes and open water through wetting-drying processes. This coupled model system is validated and successfully applied to the York River Estuary. Model results suggest that tidal marshes influence the local diurnal dissolved oxygen (DO) cycle by exporting dissolved organic carbon and high sediment oxygen demand in the marsh system through the tidal exchange. The high deposition rates of organics and diurnal DO cycle enhance the sediment release of phosphorus. On the other hand, marshes tend to decrease dissolved inorganic nitrogen in the water column by settling particulate nutrients and enhancing the denitrification process. The study demonstrates that tidal marshes exert substantial impacts on the estuarine biochemical processes. The developed tidal marsh model enhances eutrophication modeling and advances the understanding of the feedback effects between marsh biogeochemistry and estuarine eutrophication processes on a systemic scale.

54 ENVIRONMENTAL SCIENCES↗

FECM/NETL CO2 Saline Storage Cost Model CO2_S_COM 2024 (v4)

The U.S. Department of Energy's (DOE) Office of Fossil Energy and Carbon Management (FECM), in collaboration with the National Energy Technology Laboratory (NETL), has developed the FECM/NETL CO2 Saline Storage Cost Model (CO2_S_COM). This Excel-based tool provides a comprehensive framework for estimating the costs and breakeven prices associated with storing carbon dioxide (CO2) in deep saline formations. Designed from the perspective of a CO2 storage site owner, the CO2_S_COM incorporates four integrated modules—project management, financial analysis, activity cost estimation, and geological evaluation—to deliver fast, robust and actionable insights for screening project finances.

CO2 storage↗

Unlocking the Tight Oil Reservoirs of the Powder River Basin, Wyoming

The project focused on detailed geologic characterization, geomechanical studies, well completion optimization, stimulation monitoring, and field development strategies. A key aspect involved partnerships with industry and academic collaborators such as Occidental Petroleum, Southern Illinois University, Britt Rock Mechanics and Piri Technologies. Data acquisition included drilling, logging, coring, deployment of fiber optics, and microseismic monitoring. The project emphasized feedback loops for continuous model updating and integration of economic evaluations to guide development strategies.

02 PETROLEUM↗

Optimizing Deep Geothermal Drilling for Energy Sustainability in the Appalachian Basin

This study investigates the geological and geomechanical characteristics of the MIP 1S geothermal well in the Appalachian Basin to optimize drilling and address the wellbore stability issues encountered. Data from well logs, sidewall core analysis, and injection tests were used to derive elastic and rock strength properties, as well as stress and pore pressure profiles. A robust 1D-geomechanical model was developed and validated, correlating strongly with wellbore instability observations. This revealed significant wellbore breakout, widening the diameter from 12 ¼ inches to over 16 inches. Advanced technologies like Cerebro Force™ In-Bit Sensing were used to monitor drilling performance with high accuracy. This technology tracks critical metrics such as bit acceleration, vibration in the x, y, and z directions, Gyro RPM, stick-slip indicators, and bending on the bit. Cerebro Force™ readings identified hole drag caused by poor hole conditions, including friction between the drill string and wellbore walls and the presence of cuttings or debris. This led to higher torque and weight on bit (WOB) readings at the surface compared to downhole measurements, affecting drilling efficiency and wellbore stability. Optimal drilling parameters for future deep geothermal wells were determined based on these findings.

Environmental Sciences & Ecology↗

Feasibility of using reactive silicate particles with temperature-responsive coatings to enhance the security of geologic carbon storage

The large-scale deployment of geologic carbon storage will require minimizing the risk of CO 2 leakage. Here we report on a novel strategy to seal leakage pathways by injecting functionalized reactive calcium silicate (CaSiO 3 ) particles. The particles consist of a mineral silicate core within a temperature-sensitive Poly(N,N-dimethylaminoethyl methacrylate) coating that limits contact between the core and the surrounding CO 2 (aq) above a switching temperature making them relatively unreactive deep in the subsurface. If the CO 2 leaks to shallower depths and cooler temperatures, the coatings would extend and enable the silicate core to react with CO 2(aq) to form solid carbonate precipitates that reduce the permeability of the fracture or flow path. A combined laboratory- and simulation-based investigation was conducted to explore key early-stage research questions related to the feasibility of this leakage management strategy. Sand columns were injected with either coated or uncoated CaSiO 3 particles and the permeability was measured before and after exposure to CO 2 . Experiments were conducted at 35 °C and 70 °C, a range that spans the switching temperature for the coatings. Results show that at the cooler temperature the coating effectively enabled reaction and permeability reduction. At higher temperatures, the coated particles limited mineral carbonation and the permeability decrease was minimal compared to columns with uncoated particles. Here, the morphology and location of the precipitates provides insight about permeability evolution. Model simulations of 2D radial particle transport from an injection point provided insight into the effects of particle size, coating properties, and concentrations of particles in several representative rock formations.

54 ENVIRONMENTAL SCIENCES↗

Assessing Impacts on Pressure Stabilization and Leasing Acreage for CO 2 Storage Utilizing Oil Migration Concepts

Favorable geological storage for CO 2 has long been pictured as large anticlines with thick sandstones, similar to oil reservoirs in the petroleum system. Unlike oil, however, stored CO 2 does not need to be recoverable, which raises the possibility of using dissolution and residual trapping to augment the capacity of buoyant traps and tap more of the bulk rock volume. The work presented builds on that idea, asking the following question: If we inject CO 2 down to a syncline – analogous to the carrier bed in the petroleum system – how would this injection mechanism impact storage capacity and plume shape, migration, and stabilization? To address this question, we built a reservoir model, based on seismic interpretation of Middle Miocene strata, offshore Galveston, Texas. 3-D seismic and well logs were used to characterize key intervals. Reservoirs chosen for modeling are progradational-aggradational sands with mud intercalation. They have a higher degree of heterogeneity than the more conventional reservoirs commonly targeted for CO 2 storage. Modeling investigated how far the CO 2 plume would migrate under two scenarios: (1) injecting CO 2 at the base of the salt withdrawal basin (syncline scenario) and (2) injecting CO 2 at the base of the structural closure, similar to a common injection well location for EOR purposes (base scenario). For each scenario, we separately simulated injection of 30 MT of CO 2 and 60 MT of CO 2 continuously for 30 years and observe the plume and pressure evolution for 100 years after the injection stops. The simulation shows that injecting the CO 2 into a syncline limits the vertical migration of CO 2 , thus making synclinal injection more secure. In the syncline scenario, the geological layer around the injection point is more heterogeneous than the layer in the base scenario; thus, the CO 2 tends to migrate laterally. Additionally, in the syncline scenario, the plume does not even reach the upper part of the anticline, allowing us to safely store an additional amount of CO 2 into the reservoir. Furthermore, the simulation also shows that in the syncline scenario, the times needed for the reservoir to reach its stabilized pressure after the end of injections are faster. To summarize, CO 2 injection at the base of a syncline could provide additional storage, increase the safety of the project from the limited vertical plume migration, and expedite plume stabilization, which could result in the decrease of monitoring frequency as the project runs, thus lowering the operating cost of the project in the long run.

58 GEOSCIENCES↗

The importance of sub-meter-scale snow roughness on conductive heat flux of Arctic sea ice

Abstract The conductive heat flux through the snow and ice is a critical component of the mass and energy budgets in the Arctic sea ice system. We use high horizontal resolution (3–15 cm) measurements of snow topography to explore the impacts of sub-meter-scale snow surface roughness on heat flux as simulated by the Finite Element method. Simulating horizontal heat flux in a variable snow cover modestly increases the total simulated heat flux. With horizontal heat flux, as opposed to simple 1D-vertical heat flux modeling, the simulated heat flux is 10% greater than that for uniform snow with the same mean snow thickness for a 31.5 × 21 m region of sea ice (the largest region we studied). Vertical-only (1D) heat flux simulates just a 6% increase for the same region. However, this is highly dependent on observation resolution. Had we measured the snow cover at 1 m horizontal spacing or greater, simulating horizontal heat flux would not have changed the net heat flux from that simulated with vertical-only heat flux. These findings suggest that measuring and modeling snow roughness at sub-meter horizontal scales may be necessary to accurately represent horizontal heat flux on level Arctic sea ice.

Geology↗

River sinuosity describes a continuum between randomness and ordered growth

River channels are among the most common landscape features on Earth. An essential characteristic of channels is sinuosity: their tendency to take a circuitous path, which is quantified as along-stream length divided by straight-line length. River sinuosity is interpreted as a characteristic that either forms randomly at channel inception or develops over time as meander bends migrate. Studies tend to assume the latter and thus have used river sinuosity as a proxy for both modern and ancient environmental factors including climate, tectonics, vegetation, and geologic structure. But no quantitative criterion for planform expression has distinguished between random, initial sinuosity and that developed by ordered growth through channel migration. This ambiguity calls into question the utility of river sinuosity for understanding Earth’s history. We propose a quantitative framework to reconcile these competing explanations for river sinuosity. Using a coupled analysis of modeled and natural channels, we show that while a majority of observed sinuosity is consistent with randomness and limited channel migration, rivers with sinuosity ≥1.5 likely formed their geometry through sustained, ordered growth due to channel migration. This criterion frames a null hypothesis for river sinuosity that can be applied to evaluate the significance of environmental interpretations in landscapes shaped by rivers. The quantitative link between sinuosity and channel migration further informs strategies for preservation and restoration of riparian habitat and guides predictions of fluvial deposits in the rock record and in remotely sensed environments from the seafloor to planetary surfaces.

58 GEOSCIENCES↗

Nonequilibrium Physics of Multiphase Flow in Porous Media: Wettability and Disorder

Carbon capture and geologic storage, dissociation of methane hydrates in permafrost, enhanced oil recovery, and water dropout in low-temperature fuel cells, all have something in common: two or more fluids flow simultaneously through a porous medium; and the displacement of one fluid by another is often unstable (either due to gravity or viscous forces). Yet, our ability to model multiphase flow in porous media has remained a challenge. The traditional equations are unable to predict, explain, or even reproduce, the formation of the complex patterns observed in experiments. The overarching goal of the project is to develop new physical understanding of the role of wettability and disorder in multiphase flow through permeable media, and to develop new mathematical and computational models at the pore scale and at the continuum scale.

58 GEOSCIENCES↗

Tropical dry forest response to nutrient fertilization: a model validation and sensitivity analysis

Abstract. Soil nutrients, especially nitrogen (N) and phosphorus (P), regulate plant growth and hence influence carbon fluxes between the land surface and atmosphere. However, how forests adjust biomass partitioning to leaves, wood, and fine roots in response to N and/or P fertilization remains puzzling. Recent work in tropical forests suggests that trees increase fine root production under P fertilization, but it is unclear whether mechanistic models can reproduce this dynamic. In order to better understand mechanisms governing nutrient effects on plant allocation and improve models, we used the nutrient-enabled ED2 model to simulate a fertilization experiment being conducted in a secondary tropical dry forest in Costa Rica. We evaluated how different allocation parameterizations affected model performance. These parameterizations prescribed a linear relationship between relative allocation to fine roots and soil P concentrations. The slope of the linear relationship was allowed to be positive, negative, or zero. Some parameterizations realistically simulated leaf, wood, and fine root production, and these parameterizations all assumed a positive relationship between relative allocation to fine roots and soil P concentration. Model simulations of a 30-year timeframe indicated strong sensitivity to parameterization and fertilization treatment. Without P fertilization, the simulated aboveground biomass (AGB) accumulation was insensitive to the parameterization. With P fertilization, the model was highly sensitive to the parameterization and the greatest AGB accumulation occurred when relative allocation to fine roots was independent of soil P. Our study demonstrates the need for simultaneous measurements of leaf, wood, and fine root production in nutrient fertilization experiments and for longer-term experiments. Models that do not accurately represent allocation to fine roots may be highly biased in their simulations of AGB, especially on multi-decadal timescales.

Environmental Sciences & Ecology↗