Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “underground”

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 1,369 records · Page 76

A revised model for microbially induced calcite precipitation: Improvements and new insights based on recent experiments

Abstract The model for microbially induced calcite precipitation (MICP) published by Ebigbo et al. (2012) has been improved based on new insights obtained from experiments and model calibration. The challenge in constructing a predictive model for permeability reduction in the underground with MICP is the quantification of the complex interaction between flow, transport, biofilm growth, and reaction kinetics. New data from Lauchnor et al. (2015) on whole‐cell ureolysis kinetics from batch experiments were incorporated into the model, which has allowed for a more precise quantification of the relevant parameters as well as a simplification of the reaction kinetics in the equations of the model. Further, the model has been calibrated objectively by inverse modeling using quasi‐1D column experiments and a radial flow experiment. From the postprocessing of the inverse modeling, a comprehensive sensitivity analysis has been performed with focus on the model input parameters that were fitted in the course of the model calibration. It reveals that calcite precipitation and concentrations of and are particularly sensitive to parameters associated with the ureolysis rate and the attachment behavior of biomass. Based on the determined sensitivities and the ranges of values for the estimated parameters in the inversion, it is possible to identify focal areas where further research can have a high impact toward improving the understanding and engineering of MICP.

Hommel, Johannes↗

Increased stream discharge after the 3 September 2016 M w 5.8 Pawnee, Oklahoma earthquake

Earthquakes influence hydrogeological processes and properties in Earth's crust, some of which affect surface waters. Here we document increased discharge in a stream after the 3 September 2016 M w 5.8 earthquake near Pawnee, Oklahoma, an event likely induced by underground wastewater disposal. Discharge increased by an order of magnitude and remained elevated until the change was obscured by rain 1 week later. Given the earthquake magnitude and distance from the stream, by comparison with previous examples of responses to earthquakes, increased discharge after this earthquake is expected. While the mechanism increasing discharge cannot be confirmed, the observations require changes in physical properties of the subsurface. Fluid injection may thus influence hydrogeological properties of shallow groundwater systems and aquifers indirectly by inducing seismicity, if the induced seismic events are large enough.

58 GEOSCIENCES↗

Validating induced seismicity forecast models—Induced Seismicity Test Bench

Abstract Induced earthquakes often accompany fluid injection, and the seismic hazard they pose threatens various underground engineering projects. Models to monitor and control induced seismic hazard with traffic light systems should be probabilistic, forward‐looking, and updated as new data arrive. In this study, we propose an Induced Seismicity Test Bench to test and rank such models; this test bench can be used for model development, model selection, and ensemble model building. We apply the test bench to data from the Basel 2006 and Soultz‐sous‐Forêts 2004 geothermal stimulation projects, and we assess forecasts from two models: Shapiro and Smoothed Seismicity (SaSS) and Hydraulics and Seismics (HySei). These models incorporate a different mix of physics‐based elements and stochastic representation of the induced sequences. Our results show that neither model is fully superior to the other. Generally, HySei forecasts the seismicity rate better after shut‐in but is only mediocre at forecasting the spatial distribution. On the other hand, SaSS forecasts the spatial distribution better and gives better seismicity rate estimates before shut‐in. The shut‐in phase is a difficult moment for both models in both reservoirs: the models tend to underpredict the seismicity rate around, and shortly after, shut‐in.

Király‐Proag, Eszter↗

Field characterization of elastic properties across a fault zone reactivated by fluid injection

In this paper, we studied the elastic properties of a fault zone intersecting the Opalinus Clay formation at 300 m depth in the Mont Terri Underground Research Laboratory (Switzerland). Four controlled water injection experiments were performed in borehole straddle intervals set at successive locations across the fault zone. A three-component displacement sensor, which allowed capturing the borehole wall movements during injection, was used to estimate the elastic properties of representative locations across the fault zone, from the host rock to the damage zone to the fault core. Young's moduli were estimated by both an analytical approach and numerical finite difference modeling. Results show a decrease in Young's modulus from the host rock to the damage zone by a factor of 5 and from the damage zone to the fault core by a factor of 2. In the host rock, our results are in reasonable agreement with laboratory data showing a strong elastic anisotropy characterized by the direction of the plane of isotropy parallel to the laminar structure of the shale formation. In the fault zone, strong rotations of the direction of anisotropy can be observed. Finally, the plane of isotropy can be oriented either parallel to bedding (when few discontinuities are present), parallel to the direction of the main fracture family intersecting the zone, and possibly oriented parallel or perpendicular to the fractures critically oriented for shear reactivation (when repeated past rupture along this plane has created a zone).

58 GEOSCIENCES↗

Aseismic Motions Drive a Sparse Seismicity During Fluid Injections Into a Fractured Zone in a Carbonate Reservoir

An increase in fluid pressure in faults can trigger seismicity and large aseismic motions. Understanding how fluid and faults interact is an essential goal for seismic hazard and reservoir monitoring, but this key relation remains unclear. We developed an in situ experiment of fluid injections at a 10 meter scale. Water was injected at high pressure in different geological structures inside a fault damaged zone, in limestone at 280 m depth in the Low Noise Underground Laboratory (France). Induced seismicity, as well as strains, pressure, and flow rate, was continuously monitored during the injections. Although nonreversible deformations related to fracture reactivations were observed for all injections, only a few tests generated seismicity. Events are characterized by a 0.5-to-4 kHz content and a small magnitude (approximately -3.5). They are located within 1.5 m accuracy between 1 and 12 m from the injections. Comparing strain measurements and seismicity shows that more than 96% of the deformation is aseismic. The seismic moment is also small compared to the one expected from the injected volume. Moreover, a dual seismic behavior is observed as (1) the spatiotemporal distribution of some cluster of events is clearly independent from the fluid diffusion (2) while a diffusion-type pattern can be observed for some others clusters. Furthermore, the seismicity might therefore appear as an indirect effect to the fluid pressure, driven by aseismic motion and related stress perturbation transferred through failure.

58 GEOSCIENCES↗

Permeability Variations Associated With Fault Reactivation in a Claystone Formation Investigated by Field Experiments and Numerical Simulations

We studied the relation between rupture and changes in permeability within a fault zone intersecting the Opalinus Clay formation at 300 m depth in the Mont Terri Underground Research Laboratory (Switzerland). A series of water injection experiments were performed in a borehole straddle interval set within the damage zone of the main fault. A three-component displacement sensor allowed an estimation of the displacement of a minor fault plane reactivated during a succession of step rate pressure tests. The experiment reveals that the fault hydromechanical (HM) behavior is different from one test to the other with varying pressure levels needed to trigger rupture and different slip behavior under similar pressure conditions. Numerical simulations were performed to better understand the reason for such different behavior and to investigate the relation between rupture nucleation, permeability change, pressure diffusion, and rupture propagation. Here, our main findings are as follows: (i) a rate frictional law and a rate-and-state permeability law can reproduce the first test, but it appears that the rate constitutive parameters must be pressure dependent to reproduce the complex HM behavior observed during the successive injection tests; (ii) almost similar ruptures can create or destroy the fluid diffusion pathways; (iii) a too high or too low diffusivity created by the main rupture prevents secondary rupture events from occurring whereas “intermediate” diffusivity favors the nucleation of a secondary rupture associated with the fluid diffusion. However, because rupture may in certain cases destroy permeability, this succession of ruptures may not necessarily create a continuous hydraulic pathway.

58 GEOSCIENCES↗

Modeling the distribution of the endangered Jemez Mountains salamander (Plethodon neomexicanus) in relation to geology, topography, and climate

The Jemez Mountains salamander (Plethodon neomexicanus; hereafter JMS) is an endangered salamander restricted to the Jemez Mountains in north-central New Mexico, United States. This strictly terrestrial and lungless species requires moist surface conditions for activities such as mating and foraging. Threats to its current habitat include fire suppression and ensuing severe fires, changes in forest composition, habitat fragmentation, and climate change. Forest composition changes resulting from reduced fire frequency and increased tree density suggest that its current aboveground habitat does not mirror its historically successful habitat regime. However, because of its limited habitat area and underground behavior, we hypothesized that geology and topography might play a significant role in the current distribution of the salamander. We modeled the distribution of the JMS using a machine learning algorithm to assess how geology, topography, and climate variables influence its distribution. The best habitat suitability model indicates that geology type and maximum winter temperature (November to March) were most important in predicting the distribution of the salamander (23.5% and 50.3% permutation importance, respectively). Minimum winter temperature was also an important variable (21.4%), suggesting this also plays a role in salamander habitat. Our habitat suitability map reveals low uncertainty in model predictions, and we found slight discrepancies between the designated critical habitat and the most suitable areas for the JMS. Because geological features are important to its distribution, we recommend that geological and topographical data are considered, both during survey design and in the description of localities of JMS records once detected.

59 BASIC BIOLOGICAL SCIENCES↗

Enhanced prediction of Cs removal by CST from Hanford tank waste with K accountability

Abstract The treatment of Hanford tank waste is one of the most challenging environmental cleanup activities to date. To expedite the processing of liquid waste stored in underground tanks in Washington State it is necessary to remove the significant dose contributor, 137 Cs. Crystalline silicotitanate ion exchanger is currently used to remove 137 Cs from the aqueous phase of Hanford tank wastes in preparation for vitrification at the Waste Treatment and Immobilization Plant (WTP). Improving the understanding of potassium impacts on ion exchange behavior of Cs will help in the operation of a critical component of one of the most complex treatment processes in the world today. Optimization of this process can result in significant cost savings and less waste production. Toward this effort, a series of batch contact tests varied in potassium concentration were conducted to look at the impact of potassium concentration on Cs distribution. Experimental distribution ratios ( K d ) were compared to the distribution ratios predicted using the ZAM model. A significant underprediction of Cs capacity in the presence of potassium was seen with the existing model. A revision of the equilibrium constants was determined and provided a statistically better fit for determining the Cs K d values in tank waste matrices.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Hydro-mechanical measurements and x-ray imaging of sheared crystalline rock fractures from EGS Collab experiments 1 & 2

Measuring hydro-mechanical properties of natural fractures is a prerequisite for optimizing hydraulic stimulation design and well placement. We completed experiments to characterize shear on natural fractures in schist, amphibolite, and rhyolite specimens drilled from EGS Collab Project’s field sites at the Sanford Underground Research Facility (SURF) in South Dakota. A triaxial direct shear method and coupled x-ray imaging were used to perform hydroshearing and mechanical shearing at the site’s in-situ stress conditions. This produced simultaneous measurements of fracture and matrix strength, permeability, stress-dependent aperture, dilation, and friction strength. Our results identified that only a subset of the natural fractures was weak enough for hydroshearing. Generally, hydroshearing increases fracture permeability by a factor of 10 or more and the enhancement is retainable over time. However, the shear slip does not always result in permeability enhancement. High content of phyllosilicates was found to associate with exceptionally weak fractures that also exhibited poor or even negative enhancement after stimulation. Combining our measurements with site data, we can predict that most observable fractures at the two EGS Collab sites do not meet the criteria for hydroshearing before tensile opening. In some cases, the visible fractures are low permeability and as strong as the adjacent rock. To induce hydroshearing before tensile opening, injection must target known weak and favorably oriented fractures with confirmed pre-existing permeability.

58 GEOSCIENCES↗

Geologic carbon sequestration injection wells in overpressured storage reservoirs: estimating area of review

The Area of Review (AoR) under the US Environmental Protection Agency's (EPA) Class VI CO 2 injection permit is defined as the region surrounding the geologic carbon sequestration (GCS) project where underground sources of drinking water (USDWs) may be endangered. Estimation of the AoR is based on the calculated reservoir pressurization due to CO 2 injection and the associated potential to lift saline water into potable groundwater aquifers through open flow paths (e.g. wells) assuming the system is hydrostatic. In cases where the storage reservoirs are not initially hydrostatic, and in particular where they are overpressured, AoR estimation methods need to be altered. Here in this paper, we present and apply an approach to evaluating potential endangerment of USDW based on comparing brine leakage through a hypothetical open flow path in a no-injection scenario and brine leakage in a CO 2 -injection scenario. We present six possible ways to normalize injection-related leakage relative to no-injection leakage. We calculate leakage using semi-analytical solutions for single-phase flow and model reservoir pressurization and flow up (single) leaky wells located progressively farther from the injection well. For an example case of relative overpressure and using an injection-rate-based approach, results show 50-60% larger open-well-leakage rates for wells located at 2 km and 10% increase for wells located at 10 km from the injection well relative to the no-injection case. If total brine leakage is considered, the results depend strongly on the assumed pre-injection to post-injection time frames and on the methods of normalization used to calculate incremental leakage.

54 ENVIRONMENTAL SCIENCES↗

A review of geochemical–mechanical impacts in geological carbon storage reservoirs

Geological carbon storage (GCS) refers to the technology of capturing man-made carbon dioxide (CO 2 ) emissions, typically from stationary power sources, and storing such emissions in deep underground reservoirs. GCS is an approach being explored globally as a defense mechanism against climate change projections, although it is not without its critics. An important focus has been recently placed on understanding the coupling between rock–fluid geochemical alterations and mechanical changes for CO 2 storage schemes in saline aquifers. This article presents a review of the current state of knowledge regarding CO 2 -induced geochemical reactions in subsurface reservoirs, and their potential impact on mechanical properties and microseismic events at CO 2 storage sites. Therefore, this review focuses, in particular, on the current state of the art in fluid–rock interactions within the GCS context. Key issues to be addressed include geochemical reactions and the alteration of transport and mechanical properties. Specific review topics include the swelling of clays, the prediction of dissolution and precipitation reaction rates, CO 2 -induced changes in porosity and permeability, constitutive models of chemo–mechanical interactions in rock, and correlations between geochemical reactions and induced seismicity. The open questions in the field are emphasized, and new research needs are highlighted.

58 GEOSCIENCES↗

A simulation study of carbon storage with active reservoir management

As part of the Integrated Midcontinent Stacked Carbon Storage Hub (IMSCS-Hub) project led by Battelle Memorial Institute, a study was conducted to determine the feasibility of storing carbon dioxide (CO 2 ) in the stacked saline rock formations of the Sleepy Hollow Field (SHF), located in Red Willow County, southern Nebraska. A series of CO 2 injection simulation scenarios, with and without active reservoir management (ARM; brine extraction), were evaluated to investigate the feasibility of storing 50+ million tonnes (Mt) of CO 2 . The results indicated CO 2 injection combined with ARM may enable permanent storage of 50+ Mt of CO 2 . The area of review (AOR), the area in which underground sources of drinking water (USDWs) might be endangered during CO 2 injection, was assessed for the simulation scenarios. In comparison to a case without ARM, brine extraction resulted in a much smaller AOR, covering an area of 42 square miles (108.8 km 2 ), roughly one-fourth the size of an AOR resulting from CO 2 injection without ARM. Our findings, presented in this paper, indicate that CO 2 injection with ARM can improve the CO 2 storage capacity of a geologic storage complex up to 100% while also reduce the rate of pressure buildup in the subsurface, resulting in a 75% reduction in AOR. This may help in lowering carbon capture and storage project costs, risks, and effort needed to meet monitoring requirements for a storage project.

54 ENVIRONMENTAL SCIENCES↗

How Rains and Floods Become Groundwater: Understanding Recharge Pathways With Stable and Cosmogenic Isotopes

Anthropogenic climate change leads to increased precipitation intensity and exacerbated droughts in California, challenging the reliability and drought resiliency of water supply. Storing floodwater underground via managed aquifer recharge can mitigate these effects through direct infiltration or streambed infiltration. Seasonally dry streams (arroyos) already play an important part in managing groundwater recharge to the Livermore basin (CA). Understanding how, when and where stormwater and arroyo water infiltrate is critical to effectively utilise this strategy. To track water from recent storms (water year 2022–2023, WY23) into the Livermore Valley Groundwater Basin, we analysed stable water isotopes (δ 18 O and δ 2 H) in combination with naturally occurring radioactive isotopic tracers, sulphur-35 ( 35 S, t ½ = 87 days) and tritium ( 3 H, t ½ = 12.3 years). By comparing measurements of δ 18 O, 35 S and 3 H in arroyos to precipitation and groundwater, we classified the relative age and identified source of recharge to 16 wells near two arroyos. Two wells contained water with recent recharge (from WY23) from local precipitation. One well had recent recharge from variable (precipitation and imported water) sources. One well contained imported water recharge. Three wells contained water from mixed recent and older (pre-WY23) waters, from local precipitation sources. Two wells contained recent recharge from local mine settling ponds. Seven wells had older recharge from local precipitation sources. This combination of isotopes allows us to delineate where local and imported water recharges in this highly managed basin and identify locations where managed aquifer recharge is contributing to rapid groundwater infiltration. Our combined interpretation of isotopic water ages and sources in the context of land use shows that local infiltration of precipitation in open spaces is an important recharge mechanism, in addition to the managed arroyo recharge. Finally, a broader familiarity with 35 S will enable more extensive research on the infiltration of urban floodwaters.

58 GEOSCIENCES↗

Implementation and verification of a user-defined element (UEL) for coupled thermal-hydraulic-mechanical-chemical (THMC) processes in saturated geological media

Efficient and accurate modeling of the coupled thermal-hydraulic-mechanical-chemical (THMC) processes in various rock formations is indispensable for designing energy geo-structures such as underground repositories for high-level nuclear wastes. Here, this work focuses on developing and verifying an implicit finite element solver for generic coupled THMC problems in geological settings. Starting from the mass, momentum, and energy balance laws, a specialized set of governing equations and a thermoporoelastic constitutive model is derived. This system is then solved by an implicit finite element (FE) scheme. Specifically, the residuals and the Jacobians are scripted in a user-defined element (UEL) subroutine which is then combined with the general-purpose FE software Abaqus Standard to solve initial-boundary value problems. Considering the complexity of the system, the UEL development follows a stepwise manner by first solving the coupled hydraulic-mechanical (HM) and thermal-hydraulic-mechanical (THM) equations before moving on to the full THMC problem. Each implementation step consists of at least one verification test by comparing computed results with closed-form analytical solutions to ensure that the various coupling effects are correctly realized. To demonstrate the robustness of the algorithm and to validate the UEL, a three-dimensional case study is performed with reference to the in-situ heating test of ATLAS at Belgium in 1980s. A hypothetical radionuclide leakage event is then simulated by activating the chemical-concentration degree of freedom and prescribing a constant high concentration at the heater's surface. The model predicts a limited contaminated regime after six years considering both diffusion and advection effects on species transport.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

An integrated EOS, pore‐crush, strength and damage model framework for near‐field ground‐shock

Abstract An integrated Equation of State (EOS) and strength/pore‐crush/damage model framework is provided for modeling near to source (near‐field) ground‐shock response, where large deformations and pressures necessitate coupling EOS with pressure‐dependent plastic yield and damage. Nonlinear pressure‐dependence of strength up to high‐pressures is combined with a Modified Cam‐Clay‐like cap‐plasticity model in a way to allow degradation of strength from pore‐crush damage, what we call the “Yp‐Cap” model. Nonlinear hardening under compaction allows modeling the crush‐out of pores in combination with a fully saturated EOS, that is, for modeling partially saturated ground‐shock response, where air‐filled voids crush. Attention is given to algorithmic clarity and efficiency of the provided model, and the model is employed in example numerical simulations, including finite element simulations of underground explosions to exemplify its robustness and utility.

54 ENVIRONMENTAL SCIENCES↗

Visualisation and outlier detection for probability density function ensembles

Abstract Exploratory data analysis (EDA) for functional data—data objects where observations are entire functions—is a difficult problem that has seen significant attention in recent literature. This surge in interest is motivated by the ubiquitous nature of functional data, which are prevalent in applications across fields such as meteorology, biology, medicine and engineering. Empirical probability density functions (PDFs) can be viewed as constrained functional data objects that must integrate to one and be nonnegative. They show up in contexts such as yearly income distributions, zooplankton size structure in oceanography and in connectivity patterns in the brain, among others. While PDF data are certainly common in modern research, little attention has been given to EDA specifically for PDFs. In this paper, we extend several methods for EDA on functional data for PDFs and compare them on simulated data that exhibit different types of variation, designed to mimic that seen in real‐world applications. We then use our new methods to perform EDA on the breakthrough curves observed in gas transport simulations for underground fracture networks.

97 MATHEMATICS AND COMPUTING↗

The CUORE bolometric detector for neutrinoless double beta decay searches

The Cryogenic Underground Observatory for Rare Events (CUORE) is the first bolometric experiment reaching the 1-ton scale. The detector consists of an array of 988 TeO2 crystals arranged in 19 towers. The construction of the experiment and the installation of the detector was completed in August 2016. In this paper, the technical challenges of the construction, the design choices and measured performance of the electronic instrumentation are presented.

Cassina, L↗