Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “subsurface system”

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 181 records · Page 10

Estimation of hydraulic conductivity in a watershed using sparse multi-source data via Gaussian process regression and Bayesian experimental design

Enhanced water management systems depend on accurate estimation of subsurface hydraulic properties. However, geologic formations can vary significantly, so information from a single source (e.g., widely spaced boreholes) is insufficient in characterizing subsurface aquifer properties. Therefore, multiple sources of information are needed to complement the hydrogeology understanding of a region. Here, this study presents a numerical framework in which information from different measurement sources is combined to characterize the 3D random field in a multi-fidelity prediction model. Coupled with the model, a Bayesian experimental design was used to determine the best future sampling locations. The Upper Sangamon watershed in east-central Illinois was selected as the case study site, where the multi-fidelity Gaussian process model was used to estimate the hydraulic conductivity in the region of interest. Multi-source observation data were obtained from electrical resistivity and borehole pumping tests. The accuracy of the model prediction is dependent on the locations and the distribution of both high- and low-fidelity data. Furthermore, the multi-fidelity model was compared with the single-fidelity model. The uncertainties and confidence in the measurements and parameter estimates were quantified and used to design future cycles of data collection to further improve the confidence intervals.

54 ENVIRONMENTAL SCIENCES↗

DASEventNet: AI‐Based Microseismic Detection on Distributed Acoustic Sensing Data From the Utah FORGE Well 16A (78)‐32 Hydraulic Stimulation

Abstract Distributed acoustic sensing (DAS) has emerged as a promising seismic technology for monitoring microearthquakes (MEQs) with high spatial resolution. Efficient algorithms are needed for processing large DAS data volumes. This study introduces a deep learning (DL) model based on a Residual Convolutional Neural Network (ResNet) for detecting MEQs using DAS data, named as DASEventNet. The test data were collected from the Utah FORGE 16A (78)‐32 hydraulic stimulation experiments conducted in April 2022. The DASEventNet model achieves a remarkable accuracy of 100% when discriminating MEQs from noise in the raw test set of 260 examples. Surprisingly, the model identified weak MEQ signatures that have been manually categorized as noise. The decision‐making process with the model is decoded by the classic activation map, which illuminates learning features of the DASEventNet model. These features provide clear illustrations of weak MEQs and varied noise types. Finally, we apply the trained model to the entire period (∼7 days) of continuous DAS recordings and find that it discovers >5,700 new MEQs, previously unregistered in the public Silixa DAS catalog. The DASEventNet model significantly outperforms the traditional seismic method Short‐Term Average/Long‐Term Average (STA/LTA), which detected only 1,307 MEQs. The DASEventNet detection threshold is M w −1.80 compared to the minimum magnitude of M w −1.14 detected by STA/LTA. The spatiotemporal distribution of the newly identified MEQs defines an extensive stimulation zone and more accurately characterizes fracture geometry. Our results highlight the potential of DL for long‐term, real‐time microseismic monitoring that can improve enhanced geothermal systems and other activities that include subsurface hydraulic fracturing.

15 GEOTHERMAL ENERGY↗

Introduction to this special section: Geothermal energy

Geothermal energy is a global renewable resource that has the potential to provide a significant portion of baseload energy in many regions. In the United States, it has the potential to provide 8.5% of the electric generation capacity by the middle of the century. In general, geothermal systems require heat, permeability, and water to be viable for energy generation. However, with current technologies, only heat is strictly necessary in a native system. Engineered geothermal systems (EGS) introduce water into the subsurface at elevated pressures and reduced temperatures and enhance permeability through hydraulic and/or shear fracturing. Additionally, although moderate- to high-temperature resources currently dominate geothermal energy production, low-temperature resources have been utilized for direct-use cases. When well balanced and maintained, geothermal resources can produce significant amounts of heat and achieve long-term sustainability on the order of an estimated tens to hundreds of years.

15 GEOTHERMAL ENERGY↗

High-temperature and high-pressure NMR investigations of low viscous fluids confined in mesoporous systems

Abstract In this contribution, the relaxation and diffusional behaviors of low viscous fluids, water and methanol confined into mesoporous silica and controlled size pore glass were investigated. The engineered porous systems are relevant to geologically important subsurface energy materials. The engineered porous proxies were characterized by Brunauer–Emmett–Teller (BET) surface analyzer, nuclear magnetic resonance (NMR) spectroscopy, and electron microscopy (EM) to determine surface area, pore-wall protonation and morphology of these materials, respectively. The confined behavior of the low viscous fluids was studied by varying pore diameter, fluid-to-solid ratio, temperature, and pressure, and then compared to bulk liquid state. Both relaxation and diffusion behaviors for the confined fluids showed increasing deviation from pure bulk fluids as the fluid-to-solid ratio was decreased, and surface-to-volume ratio ( S / V ) was varied. Variable pressure deuteron NMR relaxation of confined D 2 O and confined methanol, deuterated at the hydroxyl or methyl positions, were performed to exploit the sensitivity of the deuteron quadrupole moment to molecular rotation. The methanol results demonstrated greater pressure dependence than those for water only in bulk. The deviations from bulk liquid behavior arise from different reasons such as confinement and the interactions between confined fluid and the nano-pore wall. The results of the present report give insight into the behavior of low viscosity fluid in nano-confined geometries under different state conditions.

Chemistry↗

Real-Time Monitoring of Fracture Dynamics with a Contrast Agent-Assisted Electromagnetic Method

In collaboration with the Advanced Energy Consortium, our team has previously demonstrated that the placement of electrically active proppants (EAPs) in a hydraulic fracture surveyed by electromagnetic (EM) methods can enhance the imaging of the stimulated reservoir volumes during hydraulic fracturing. That work culminated in constructing a well-characterized EAP-filled fracture anomaly at the Devine field pilot site (DFPS). In subsequent laboratory studies, we observed that the electrical conductivity of our EAP correlates with changes in pressure, salinity, and flow. Thus, we postulated that the EAP could be used as an in-situ sensor for the remote monitoring of these changes in previously EAP-filled fractures. This paper presents our latest field data from the DFPS to demonstrate such correlations at an intermediate pilot scale. We conducted surface-based EM surveys during freshwater (200 ppm) and saltwater (2,500 ppm) slug injections while running surfaced-based EM surveys. Simultaneously, we measured the following: 1) bottomhole pressure and salinity in five monitoring wells; 2) injection rate using high-precision data loggers; 3) distributed acoustic sensors in four monitoring wells; and 4) tiltmeter data on the survey area. We demonstrated that injections into an EAP-filled fracture could be successfully coupled with real-time electric field measurements on the surface, leading to remote monitoring of dynamic changes within the EAP-filled fracture. Furthermore, by comparing the electrical field traces with the bottomhole pressure, flow rate, and salinity, we concluded that the observed electric field in our study is influenced by fracture dilation and flow rate. Salinity effect was observed when saltwater was injected. EM simulations solely based on assumptions of fracture conductivity changes during injection did not reproduce all of the measured electric field magnitudes. Preliminary estimates showed that including streaming potential in our geophysical model may be needed to reduce the simulation mismatch. The methods developed and demonstrated during this study will lead to a better understanding of the extent of fracture networks, formation stress states, fluid leakoff and invasion, characterizations of engineered fracture systems, and other applications where monitoring subsurface flow tracking is deemed important.

02 PETROLEUM↗

Rock Valley Accelerated Weight Drop Seismic Data Processing and Picking of P-wave and S-wave Arrival Times

Rock Valley, in the southern end of the Nevada National Security Site, hosts a fault system that was responsible for a shallow (< 3 km below surface ) magnitude 3.7 earthquake in May 1993. In order to better understand this system, seismic properties of the shallow subsurface need to be better constrained. In April and May of 2021, accelerated weight drop (AWD) active-source seismic data were recorded in order to measure P- and S-wave travel-times for the area. This report describes the processing and phase picking of the recorded seismic waveforms. In total, we picked 7,982 P-wave arrivals at offsets up to ~2500 m, and 4,369 S-wave arrivals at offsets up to ~2200 m. These travel-time picks can be inverted for shallow P-wave and S-wave velocity structure in future studies.

58 GEOSCIENCES↗

DISSPATCH H2: Assessing Subsurface Storage Potential for Hydrogen

This project evaluates the technical feasibility of long-term underground hydrogen storage in the Anadarko Basin, integrating subsurface science, reservoir engineering, and system-level analysis to build a cohesive understanding of how hydrogen may be stored and recovered in porous media. Although salt caverns currently dominate commercial-scale hydrogen storage, their limited geographic availability necessitates alternative options that can support a more resilient energy system. Depleted oil and gas reservoirs offer a widely distributed and potentially high-capacity solution, but their suitability depends on formation-specific geologic, petrophysical, microbial, and operational factors. This study provides the first coordinated assessment of these factors within the Simpson Group of Oklahoma.

03 NATURAL GAS↗

Reservoir Characterization and Techno-Economic Analysis of Enhanced Geothermal and Closed-Loop Systems in the Wattenberg Field of the Denver-Julesburg Basin, Colorado

This report provides details of reservoir characterization and techno-economic analysis for enhanced geothermal systems and closed-loop geothermal systems in the Denver-Julesburg Basin Wattenberg area. The analysis contributes to the Geothermal Limitless Approach to Drilling Efficiencies (GLADE) project, a multi-institutional research initiative focused on advancing geothermal energy development by reducing drilling costs and improving the rate of penetration. The GLADE project brings together national laboratories, academic institutions, and industry partners, including the U.S. Department of Energy's Geothermal Technologies Office, Occidental Petroleum, the National Laboratory of the Rockies (NLR), Los Alamos National Laboratory, Colorado School of Mines, Louisiana State University, Texas A&M University, and several drilling technology companies. As part of this effort, NLR developed and applied slender-body theory modeling tools and the GEOPHIRES simulator to support techno-economic analysis. The slender-body theory wellbore simulator evaluated thermal performance in U-loop, Eavor-type multilateral, and enhanced geothermal system well configurations. Results indicated that an optimized single U-loop design delivers an outlet temperature of 168.8 degrees Celsius,stabilizing at a final outlet temperature of 135.3 degrees Celsius. Heat production also improved significantly, with an average thermal output of 16.9 MWth and a final thermal output of 11.0 MWth. These results provided critical baseline data for selecting optimized systems for techno-economic analysis. Building on these outputs, NLR applied the GEOPHIRES techno-economic simulator to estimate the levelized cost of electricity for both the Eavor multilateral and enhanced geothermal systems. The sensitivity analysis showed that drilling costs exert a particularly strong impact on the levelized cost of electricity of multilateral systems. Therefore, advances in drilling technology and improved drilling rates could yield significant gains in the economic performance of these systems. In this context, the GLADE project directly supports the U.S. Department of Energy's mission to reduce the cost of geothermal energy and accelerate the deployment of next-generation geothermal technologies.

15 GEOTHERMAL ENERGY↗

Impact of organic acids and sulfate on the biogeochemical properties of soil from urban subsurface environments

Urban subsurface environments are often different from undisturbed subsurface environments due to the impacts of human activities. For example, deterioration of underground infrastructure can introduce elevated levels of Ca, Fe, and heavy metals into subsurface soils and groundwater. Likewise, leakage from sewer systems can lead to contamination by organic C, N, S, and P. However, the impact of these organic and inorganic compounds on biogeochemical processes including microbial redox reactions, mineral transformations, and microbial community transitions in urban subsurface environments is poorly understood. Here we conducted a microcosm experiment with soil samples from an urban construction site to investigate the possible biotic and abiotic processes impacted when sulfate and acetate or lactate were introduced into an urban subsurface environment. In the top-layer soil (0-0.3 m) microcosms, which were highly alkaline (pH > 10), the major impact was on abiotic processes such as secondary mineral precipitation. In the mid-layer (2-3 m) soil microcosms, the rate of Fe(III)-reduction and the amount of Fe(II) produced were greatly impacted by the specific organic acid added, and sulfate-reduction was not observed until after Fe(III)-reduction was complete. Near the end of the incubation, some genera related to syntrophic acetate oxidation and methanogenesis were observed in the lactateamended microcosms. In the bottom-layer (7-8 m) soil microcosms, the rate of Fe(III)-reduction and the amount of Fe(II) produced were affected by the concentration of amended sulfate. Sulfate-reduction was concurrent with Fe(III)-reduction, suggesting that Fe(II) production was likely due to abiotic reduction of Fe(III) by sulfide produced by microbial sulfate reduction. The slightly acidic initial pH (~5.8) of the mid-soil system was a major factor controlling sequential microbial Fe(III) and sulfate reduction versus parallel Fe(III) and sulfate reduction in the bottom soil system, which had a neutral initial pH (~7.2). Finally, 16S rRNA gene-based community analysis revealed a variety of indigenous microbial groups including alkaliphiles, dissimilatory iron and sulfate reducers, syntrophes, and methanogens tightly coupled with, and impacted by, these complex abiotic and biogeochemical processes occurring in urban subsurface environments.

54 ENVIRONMENTAL SCIENCES↗

Physical Controls on Irrigation Return Flow Contributions to Stream Flow in Irrigated Alluvial Valleys

Irrigation can be a significant source of groundwater recharge in many agricultural regions, particularly in arid and semi-arid climates. Once infiltrated, irrigation recharge can travel via subsurface flowpaths that return to the river system in a lagged manner, supplementing natural streamflow weeks, months, or even years from when the irrigation was applied. In regions that experience low flows during summer and early fall, return flows can be a significant source of supplementary streamflow. Many water planning and operations models either ignore return flows or roughly approximate them with analytical solutions. Thus, return flows represent an important but often overlooked component of the hydrological exchange and overall water balance in agricultural regions. This study uses groundwater models to explore a wide range of factors that control irrigation return flow timing in irrigated alluvial valleys. A sensitivity analysis approach is used to assess how factors such as the extent of irrigated land adjacent to a stream, irrigation recharge rate, aquifer hydraulic conductivity, aquifer thickness, water table configuration, and seasonal fluctuations in stream stage control the timing of subsurface return flows. Modeling is conducted using MODFLOW models representing an irrigated alluvial valley adjacent to a stream. While a simplification of the full complexity in real systems, the models are a significant advancement from the analytical solution and provide new insight into the timescales of return flows over a broad range of possible conditions. To contextualize our modeling results, they are compared to an analytical solution commonly used for approximating return flows to evaluate its performance. Our findings show what factors and conditions influence return flow timing and control whether they contribute to streamflow over short term (months) or longer term (seasonal) time scales.

Ferencz, Stephen B.↗

Selecting Appropriate Model Complexity: An Example of Tracer Inversion for Thermal Prediction in Enhanced Geothermal Systems

Abstract A major challenge in the inversion of subsurface parameters is the ill‐posedness issue caused by the inherent subsurface complexities and the generally spatially sparse data. Appropriate simplifications of inversion models are thus necessary to make the inversion process tractable and meanwhile preserve the predictive ability of the inversion results. In this study, we investigate the effect of model complexity on fracture aperture inversion and thermal performance prediction in a field‐scale EGS model. Principal component analysis was used to map the aperture field to a low‐dimensional latent space. The complexity of the inversion model was quantitatively represented by the percentage of total variance in the original aperture fields preserved by the latent space. Tracer, pressure and flow rate data were used to invert for fracture aperture through an ensemble‐based inversion method, and the inferred aperture field was used to predict thermal performance. With an over‐simplified aperture model, ensemble collapse occurred. The inverted aperture models failed to resolve necessary flow and transport features, leading to a biased thermal performance prediction. A complex aperture model involved excessive features and was prone to overinterpreting the inversion data. Both the tracer/pressure/flow rate data reproduction and thermal prediction showed significant uncertainties, making it difficult to properly estimate long‐term thermal performance. Fortunately, our results indicate that there exists an appropriate model complexity which can simultaneously match inversion data and predict thermal performance with an acceptable uncertainty. The quality of the fit of tracer data appears to be a useful indicator of such an appropriate model complexity.

15 GEOTHERMAL ENERGY↗

Deep-branching acetogens in serpentinized subsurface fluids of Oman

Little is known of acetogens in contemporary serpentinizing systems, despite widely supported theories that serpentinite-hosted environments supported the first life on Earth via acetogenesis. To address this knowledge gap, genome-resolved metagenomics was applied to subsurface fracture water communities from an area of active serpentinization in the Samail Ophiolite, Sultanate of Oman. Two deeply branching putative bacterial acetogen types were identified in the communities belonging to the Acetothermia (hereafter, types I and II) that exhibited distinct distributions among waters with lower and higher water–rock reaction (i.e., serpentinization influence), respectively. Metabolic reconstructions revealed contrasting core metabolic pathways of type I and II Acetothermia, including in acetogenic pathway components (e.g., bacterial- vs. archaeal-like carbon monoxide dehydrogenases [CODH], respectively), hydrogen use to drive acetogenesis, and chemiosmotic potential generation via respiratory (type I) or canonical acetogen ferredoxin-based complexes (type II). Notably, type II Acetothermia metabolic pathways allow for use of serpentinization-derived substrates and implicate them as key primary producers in contemporary hyperalkaline serpentinite environments. Phylogenomic analyses indicate that 1) archaeal-like CODH of the type II genomes and those of other serpentinite-associated Bacteria derive from a deeply rooted horizontal transfer or origin among archaeal methanogens and 2) Acetothermia are among the earliest evolving bacterial lineages. The discovery of dominant and early-branching acetogens in subsurface waters of the largest near-surface serpentinite formation provides insight into the physiological traits that likely facilitated rock-supported life to flourish on a primitive Earth and possibly on other rocky planets undergoing serpentinization.

Science & Technology - Other Topics↗

Towards real-time monitoring: data assimilated time-lapse full waveform inversion for seismic velocity and uncertainty estimation

SUMMARY Rapid development of time-lapse seismic monitoring instrumentations has made it possible to collect dense time-lapse data for tomographically retrieving time-lapse (even continuous) images of subsurface changes. While traditional time-lapse full waveform inversion (TLFWI) algorithms are designed for sparse time-lapse surveys, they lack of effective temporal constraint on time-lapse data, and, more importantly, lack of the uncertainty estimation of the TLFWI results that is critical for further interpretation. Here, we propose a new data assimilation TLFWI method, using hierarchical matrix powered extended Kalman filter (HiEKF) to quantify the image uncertainty. Compared to existing Kalman filter algorithms, HiEKF allows to store and update a data-sparse representation of the cross-covariance matrices and propagate model errors without expensive operations involving covariance matrices. Hence, HiEKF is computationally efficient and applicable to 3-D TLFWI problems. Then, we reformulate TLFWI in the framework of HiEKF (termed hereafter as TLFWI-HiEKF) to predict time-lapse images of subsurface spatiotemporal velocity changes and simultaneously quantify the uncertainty of the inverted velocity changes over time. We demonstrate the validity and applicability of TLFWI–HiEKF with two realistic CO2 monitoring models derived from Frio-II and Cranfield CO2 injection sites, respectively. In both 2-D and 3-D examples, the inverted high-resolution time-lapse velocity results clearly reveal a continuous velocity reduction due to the injection of CO2. Moreover, the accuracy of the model is increasing over time by assimilating more time-lapse data while the standard deviation is decreasing over lapsed time. We expect TLFWI-HiEKF to be equipped with real-time seismic monitoring systems for continuously imaging the distribution of subsurface gas and fluids in the future large-scale CO2 sequestration experiments and reservoir management.

58 GEOSCIENCES↗

Geophysical Monitoring Shows that Spatial Heterogeneity in Thermohydrological Dynamics Reshapes a Transitional Permafrost System

Climate change is causing rapid changes of Arctic ecosystems. Yet, data needed to unravel complex subsurface processes are very rare. Using geophysical and in-situ sensing, this study closes an observational gap associated with thermohydrological dynamics in discontinuous permafrost systems. It highlights the impact of vegetation and snow thickness distribution on subsurface thermohydrological properties and processes. Large snow accumulation near tall shrubs insulates the ground and allows for rapid and downward heat flow. Thinner snowpack above graminoid results in surficial freezing and prevents water from infiltrating into the subsurface. Analyzing short term disturbances, we found that lateral flow could be a driving factor in talik formation. Inter-annual measurements show that deep permafrost temperatures increased by about 0.2°C over two years. The results, which suggest that snow-vegetation-subsurface processes are tightly coupled, will be useful for improving predictions of Arctic feedback to climate change, including how subsurface thermohydrology influences CO 2 and CH 4 fluxes.

58 GEOSCIENCES↗