SEARCH · Engineering Papers
Results for “Geothermal energy”
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.
Electrodialysis through nafion membranes coated with polyelectrolyte multilayers yields >99% pure monovalent ions at high recoveries
Separation of monovalent and multivalent ions is important in processes ranging from water softening to recycling of Li. Previous studies show that cation-exchange membranes coated with multilayer polyelectrolyte films exhibit monovalent/divalent ion selectivities as high as 1000 in electrodialysis (ED) under ideal conditions with low ion recoveries. This work shows that even at 60-80% monovalent-ion recoveries, ED through Nafion membranes coated with protonated poly(allylamine)/poly(4-styrenesulfonate) multilayers can separate equimolar K + and Mg 2+ or Li + and Mg 2+ to give 99.5% pure monovalent ion. Moreover, the use of a low-volume receiving phase allows ED to concentrate the monovalent ion. Typical current efficiencies are around 70% for these separations. Selectivities appear to increase when the membrane-coating procedure uses a higher polyelectrolyte concentration in the deposition solution. Because the limiting current declines with ED time due to depletion of the monovalent ion in the source phase, use of a gradually decreasing applied current enhances current efficiency compared to ED at a constant high current. Further, the membranes also show Li + /Mg 2+ selectivities around 100 when Mg 2+ is in a 20-fold excess in the source phase at an ionic strength of 0.6 M, but selectivity decreases at even higher ionic strengths.
Numerical study of multi-component flow and mixing in a scaled fission product venting system
Numerical simulation of isothermal multi-component gas flow was conducted to study fluid flow and mixing in a scaled low-pressure, low-temperature test facility for the fission product venting system (FPVS) of a gas cooled fast reactor (GFR). The geometry of FPVS test facility was an open loop including gradual expansion coupling and two 90°pipe elbows. First, the validation of large eddy simulation (LES) wall-adapting local eddy-viscosity (WALE) turbulent model was performed for transitional flow in a circular pipe with gradual expansion. The reactingFOAM solver in OpenFOAM v8 was employed. By introducing appropriate turbulence disturbance at the flow inlet, such as turbulence intensity, the numerical results showed good agreement in (1) the velocity profile downstream of the pipe expansion measured by particle image velocimetry (PIV) in this study, and (2) the reattachment length reported in the literature. Using this validated model, the numerical simulation of flow and component mixing was performed for the FPVS with an inlet flowrate corresponding to a Reynolds number of 2,400 to investigate the flow behavior and component mixing. Standard deviation of the mass fraction of component, referred to as absolute mixing index (AMI), was calculated to quantify the mixing of components, showing that the component mixing length determined by AMI is related to the reattachment length downstream of the expansion. Additionally, this study was able to identify the best location in the FPVS test facility to measure the velocity and concentration profies where multicomponent flow is well mixed.
Silica–Biomacromolecule Interactions: Toward a Mechanistic Understanding of Silicification
Not Available
Potential Link Between 2020 Mentone, West Texas M5 Earthquake and Nearby Wastewater Injection: Implications for Aquifer Mechanical Properties
Abstract The M5 Mentone earthquake that occurred on March 26, 2020, was the largest event recorded over the last 2 decades in West Texas within the Delaware Basin, a U.S. major petroleum‐producing area. Also, numerous hydrofracturing and wastewater disposal wells are spread across this region. Within a 30 km distance to mainshock, eight class‐II injection wells for industrial wastewater disposal target the deep porous Ellenburger aquifer at an average rate of 1.36 × 10 6 barrel (BBL) per month during 2012–2020. Poroelastic models of fluid diffusion show these nearby injectors collectively imparted up to 80.5 kPa of Coulomb stress at the mainshock location, capable of triggering this M5 event. Assuming the Mentone event occurs when pore‐pressure increase is maximum, the time delay between peak injection and the M5 occurrence corresponds with an optimal permeability of 6.76 × 10 ‐14 m 2 for the Ellenburger aquifer layer, in agreement with independent estimates.
Frequency–Magnitude Statistics of Laboratory Foreshocks Vary With Shear Velocity, Fault Slip Rate, and Shear Stress
Understanding the nucleation phase of earthquakes is key for advancing earthquake hazard assessment and improving earthquake early warning systems. However, little progress has been made in this area due to a poor understanding of nucleation processes and incomplete seismic and fault zone measurements. The ability to integrate measured fault zone properties with seismic data could significantly improve our understanding of how earthquakes begin and whether there are systematic variations in seismic properties preceding failure. In this work, we use high-resolution laboratory measurements of fault zone properties along with acoustic emission data to document temporal variations of foreshock properties. Our data show that foreshock size increases with shear stress, loading rate, and fault slip rate. We propose that the preseismic fault slip rate and fault zone thickness (i.e., porosity) work in concert to modulate foreshock properties.
Predicting Thermal Performance of an Enhanced Geothermal System From Tracer Tests in a Data Assimilation Framework
Abstract Predicting the thermal performance of an enhanced geothermal system (EGS) requires a comprehensive characterization of the underlying fracture flow patterns from practically available data such as tracer data. However, due to the inherent complexities of subsurface fractures and the generally insufficient geological/geophysical data, interpreting tracer data for fracture flow characterization and thermal prediction remains a challenging task. The present study aims to tackle the challenge by leveraging a data assimilation method to maximize the utilization of information inherently contained in tracer data, and meanwhile maintain the flexibility to handle various uncertainties. A tracer data interpretation framework was proposed with the following three components integrated: (a) We use principal component analysis (PCA) to reduce the dimensionality of model parameter space. (b) We use ES‐MDA (ensemble smoother with multiple data assimilation) to invert for fracture aperture/flow fields and obtain posterior model ensembles for uncertainty quantification. Various data types are assimilated jointly to improve the predictive ability of the posterior ensemble. (c) The inverted fracture aperture fields are then incorporated into reservoir models to predict thermal performance. We developed a field‐scale EGS model to verify the ability of the framework to characterize highly heterogeneous fracture aperture/flow fields and predicting thermal performance. We also applied the framework to a mesoscale field experiment to demonstrate its potential application in real‐world geothermal reservoirs. The results indicate that the proposed framework can effectively retrieve fracture flow information from tracer data for thermal prediction and uncertainty quantification, and thus provide informative guidance for EGS optimization and risk management.
Non-Synchronous Rotation on Europa Driven by Ocean Currents
It has been suggested that the ice shell of Jupiter's moon Europa may drift non-synchronously due to tidal torques. Here we argue that torques applied by the underlying ocean are also important and can result in non-synchronous rotation. The resulting spin rate can be slightly slower than the synchronous angular rate that would have kept the same point of the ice shell facing Jupiter. We develop an ice shell rotation model, driven by ocean stress calculated using a high-resolution state-of-the-art ocean general circulation model, and take into account the viscoelastic deformation of the ice shell. We use the ice shell model results together with observed limits on the ice shell drift speed to constrain ice shell parameters such as effective viscosity, which is currently uncertain by at least four orders of magnitude. Our results suggest, at best, sluggish ice shell convection. Depending on the relaxation time scale of the ice shell and on the ocean currents, the ice shell may exhibit negligible drift, constant drift, or oscillatory drift superimposed on random fluctuations. The expected rotation rate exceeds ~30 m/yr; future spacecraft observations can be used to test these predictions and yield insight into the properties of the ice shell and underlying ocean.
Association Between Injection and Microseismicity in Geothermal Fields With Multiple Wells: Data-Driven Modeling of Rotokawa, New Zealand, and Húsmúli, Iceland
Understanding injection-induced microseismicity in geothermal systems can provide insight into reservoir connectedness. However, fault and reservoir complexity are difficult to represent in simple analytical models, which makes it difficult to discern clear relationships from incidental associations. Here, we have used data-driven models to study how fluid injection and microseismicity are related in the Rotokawa (New Zealand) and Hellisheiði (Iceland) geothermal fields. We tested two classes of model: (a) lagged linear regression of seismicity rate as a function of well injection rates; and (b) systematic extraction of injection time series features that are then evaluated for associations with the seismicity. These models allowed us to determine which wells had the greatest correlation with microseismicity and to explain this association in a reservoir context. Finally, exploring different data types and transformations, we were unable to establish a link between rapid changes in injection rate and seismicity spikes, as suggested by some theoretical models.
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.
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.
Role of Intersections in Fracture Connectivity
Networks of intersecting fractures often provide the flow paths through subsurface reservoirs. Assessing network connectivity is challenging because fracture intersections compose a vanishingly small fraction of the network void volume. In this paper, motivated by 3D X-ray imaging of the simplest element of fracture network, that is, two orthogonal fractures, we perform a percolation and finite-size scaling analysis to study the connectivity provided by fracture intersections. The conditions when an intersection enhances connectivity across a sample depend on spatial correlations in the fracture aperture distributions, on the stress state, and on the direction of flow. Here we consider three flow directions: (a) across intersections, (b) parallel to intersections and (c) around corners. For (a), intersections provide minimal enhancement of connectivity because they contribute little additional void area. For (b), intersections increase the probability of a connected path near threshold by enabling 3D connected pathways that are not possible in parallel fractures. Flow around corners, (c), is fundamentally the result of the intersection connecting two fractures in series and spatial correlations are broken around corners, suppressing the connectivity relative to (a). When the connected fractures are stressed equally, a joint percolation threshold emerges that continues to have scale invariance. However, when the fractures are stressed unequally, the system has mixed percolation without clearly defined percolation thresholds. In all cases, percolation probabilities are found to be scale dependent which has important consequences for the connectivity of larger fracture networks composed of the fundamental element studied here.
Crustal permeability generated through microearthquakes is constrained by seismic moment
We link changes in crustal permeability to informative features of microearthquakes (MEQs) using two field hydraulic stimulation experiments where both MEQs and permeability evolution are recorded simultaneously. The Bidirectional Long Short-Term Memory (Bi-LSTM) model effectively predicts permeability evolution and ultimate permeability increase. Our findings confirm the form of key features linking the MEQs to permeability, offering mechanistically consistent interpretations of this association. Transfer learning correctly predicts permeability evolution of one experiment from a model trained on an alternate dataset and locale, which further reinforces the innate interdependency of permeability-to-seismicity. Models representing permeability evolution on reactivated fractures in both shear and tension suggest scaling relationships in which changes in permeability (Δk) are linearly related to the seismic moment (M) of individual MEQs as Δk ∝ M. This scaling relation rationalizes our observation of the permeability-to-seismicity linkage, contributes to its predictive robustness and accentuates its potential in characterizing crustal permeability evolution using MEQs.
Geologic stress modulates fluid mixing at fracture intersections
Fracture intersections are critical links that enable flow and transport in subsurface fracture networks, and their behavior strongly influences fluid mixing in a network. Although all subsurface fractures are subjected to geological stress, we lack a fundamental understanding of how fracture intersection geometry evolves under stress and how these changes influence fluid mixing. Here, we combine 3D printing, 3D X-ray tomographic imaging, and 3D pore-scale numerical simulations to reveal stress-induced changes in intersection geometry and their impact on mixing. Mixing is found to be strongly affected by partial closure of an intersection under stress. As an intersection closes, the void area for fluid flow and diffusion decreases leading to substantial deviations between conventional mixing models and full pore-scale modeling. To address this, we propose a modified mixing model that accounts for intersection deformation, which is essential for accurate modeling of solute transport and mixing through fracture networks.
In situ transmission electron microscopy observations of CaCO 3 crystallization onto polysaccharide-coated nanoparticles
Polysaccharides and proteoglycans are widely associated with the organic matrix at sites of CaCO 3 biomineralization, and previous studies indicate that these macromolecules may confer greater roles in mineral nucleation than previously recognized. This investigation uses in situ liquid-phase transmission electron microscopy (LP-TEM) to observe CaCO 3 nucleation onto aminated silica (SiO 2 –NH 3 + ) nanoparticles treated with a layer of chitosan (near-neutral derivative of chitin) or heparin (a carboxylated and highly sulfated glycosaminoglycan). In the absence of polysaccharides, few CaCO 3 particles formed and exhibited mobility. However, the SiO 2 –NH 3 + nanoparticles were enveloped in a region of higher mass density relative to the bulk solution, suggesting the development of a local solute-rich environment that surrounds the charged NH 3 + groups. The heparin- or chitosan-coated silica particles also exhibited regions of higher mass density around the nanoparticles. In the presence of these polysaccharide coatings, we observed the nucleation of abundant CaCO 3 particles whereby the polyanionic heparin promoted more nucleation than the weakly cationic chitosan. Many crystallites appeared to form at the polysaccharide–TEM cell membrane–solution interface, further indicating interfacial and macromolecule-specific control on crystallization. The combined results demonstrate that chitosan and heparin have an appreciable effect on the timing, size, and location of CaCO 3 nucleation compared to the polysaccharide-free nanoparticles.
AnisONet: A deep neural operator-based anisotropic permeability upscaler from pore to Darcy scale
Directional permeability variations, which govern directional fluid flow in porous media with anisotropy, are important to accurately predict flow behavior, reactive transport, and fluid–solid interactions for various processes such as enhanced geothermal systems, energy storage devices, and biological systems. However, the intricate architecture of porous media makes it difficult to predict directional permeabilities. In this work, we present a novel machine learning (ML) framework, AnisONet, built upon an integration of a convolutional neural network, Swin transformer, and the deep operator network architecture, designed to predict anisotropic permeability and upscale predictions to larger spatial domains. First, AnisONet was evaluated with three classes of two-dimensional (2D) porous media, including synthetic circular and elliptical grains and natural sandstone grains from micro-computed tomography images. A lattice Boltzmann model (LBM) was used to calculate directional permeabilities at every 10° angle, producing 19 data points per image of porous media. AnisONet is then trained to predict permeability as a function of rotation angle. AnisONet showed strong predictive capability of directional permeability. Second, we tested our model for five upscaling cases with a large image size in the finite-element method (FEM) for 2D Darcy flow with various permeability tensor construction methods. Overall, upscaled permeability tensors in FEM simulations produce a reasonably good match with LBM results, highlighting the importance of selecting appropriate tensor formation strategies for accurate permeability upscaling. AnisONet, as a directional permeability estimator, could be further developed for more complex geometries, with the potential to develop a foundational ML model for various applications in porous media.
Sub-Degree-Per-Hour MEMS Gyroscope for Measurement While Drilling at 300°C
The orientation module of MWD tool offers the critical drill bit orientation information, including azimuth, inclination and toolface in order to control the path of wellbore. Although high end navigation grade gyroscope can meet the accuracy requirement for azimuth finding, the requirements of cost, size and reliability under harsh environment have largely limited the deployment of gyroscopes in MWD tools. To overcome these constraints, GE Research has developed a low cost, MEMS based Multi-Ring Gyroscope (MRG) capable of azimuth-seeking in MWD applications at 300°C. The MRG prototype has been demonstrated to achieve angular random walk (ARW) better than 0.003 deg/rt(hr) and bias instability of 0.01 deg/hr, capable of meeting azimuth finding accuracy better than 0.25 deg. It has further been tested to remain fully operational at 300 °C.
Development of a 250 Degree C 15kV Supercascode Switch Using SiC JFET Technology
Tetra Corporation is developing a pulsed power drilling system for the geothermal industry. The system requires a high temperature switching technology to replace the existing technology based on Si thyristors with a temperature limit of 150 degrees C. This work focuses on the development of such switch based on a Supercascode configuration using SiC JFET technology. The development presented contains thermal models, packaging design and processes, circuit simulations and prototypes of both the Supercascode switch and the high temperature packaged SiC JFETs.