Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geoscience”

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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Underground hydrogen storage leakage detection and characterization based on machine learning of sparse seismic data

Underground hydrogen storage (UHS) is considered as a scalable approach for massive storage and seasonal extraction of hydrogen (H 2 ). Although conventional leakage detection and characterization methods based on time-lapse seismic imaging and inversion generally apply to H 2 leakage detection problem, a high-fidelity yet cost effective geophysics approach is still missing to reliably inform leakage location and properties based on very sparse data. In response, we develop a novel supervised machine learning method to detect and characterize H 2 leakage from UHS. The input to our neural network are sparse time-lapse seismic waveforms, while the output from the neural network includes the spatial location and physical properties of a H 2 leakage. Here, we generate high-quality time-lapse waveforms using the elastic-wave equations to train the neural network. We train and validate our machine learning model and find that it attains high accuracy in using extremely sparse time-lapse seismic data to detect and characterize H 2 leakage. Our investigation is the first systematic study that focuses on applying machine learning to subsurface H 2 leakage detection and characterization and could potentially serve as a cost-effective geophysical tool for underground hydrogen leakage detection and characterization with high fidelity.

08 HYDROGEN↗
DOE OSTI2026

Economic assessment of seismic monitoring for underground hydrogen storage

Underground hydrogen storage (UHS) plays a key role in the energy landscape. However, like other subsurface engineering technologies, UHS may cause leakage into the groundwater or atmosphere and possibly induce local seismicity. To reduce these risks, seismic monitoring could be a viable technique to track the UHS plume, detect leakages, and locate induced seismicity events. Seismic monitoring has been proposed to safely monitor UHS, but research in this area is still new and requires field studies. Lab and theoretical studies have demonstrated the validity of seismic monitoring for UHS. Therefore, it is imperative to analyze the economic feasibility of seismic monitoring for UHS. Hence, we develop a cost model and open-source Python code for seismic monitoring that considers types of seismometers, comprehensive operational scenarios, detection thresholds, and long-term leakage monitoring. A case study is further provided to validate the cost model on reservoir simulations of UHS. We find that the levelized cost for a 10-year operating UHS site will range on the order of ∼0.003 $\$$/kg. The methods developed in this study could also be applied to the monitoring of groundwater, gas, and/or wastewater injection.

08 HYDROGEN↗
DOE OSTI2023

Accuracy and precision of ESI-Orbitrap-IRMS observations of hours to tens of hours via reservoir injection

Here, orbitrap isotope ratio mass spectrometry (Orbitrap-IRMS) has recently been applied to high-precision, natural-abundance isotope ratio measurements of a diverse range of compounds, including amino acids, oxyanions, fatty acids, and metals. These measurements can characterize many isotope ratios simultaneously at high (≈1.0‰) precision. In a successful experiment, observed precision will track the shot-noise limit and be limited by experimental time. Some isotope ratios, for example those involving 17 O in organic compounds or multiply-substituted (‘clumped’) isotopologues, require experimental times of hours to tens of hours to achieve desired precision, while current sample introduction techniques focus on observations on the order of seconds to tens of minutes. In this study, we characterize Orbitrap-IRMS performance for three long duration measurements (individual acquisitions ≥1 h and as long as 24 h) using an automated reservoir injection system coupled to a Q Exactive HF Orbitrap with an electrospray ionization (ESI) source. First, we characterize long-term intra-measurement stability through a 24-h long measurement of acetone. We report the following isotope ratios and precisions (as acquisition errors, errors on the observed ratio within this measurement (σ AE )): 13 C/ 12 C (σ AE = 0.07‰), 17 O/ 16 O (σ AE = 1.1‰), 18 O/ 16 O (σ AE = 0.3‰), and 13 C 13 C/ 12 C (σ AE = 0.65‰). The σ AE of each tracks the shot noise limit throughout and is limited by the challenging conditions (high resolution and low numbers of ions per scan) required for 17 O/ 16 O measurement in the presence of 13 C via Orbitrap. Second, we characterize inter-measurement stability via a sequence of seven 75-min analyses of perchlorate. We observe the following ratios and acquisition errors: 37 Cl/ 35 Cl (σ AE = 0.09‰); 17 O/ 16 O (σ AE = 1.6‰); 18 O/ 16 O (σ AE = 0.7‰), 37 Cl 17 O/ 35 Cl 16 O (σ AE = 2.7‰), and 37 Cl 18 O/ 35 Cl 16 O (σ AE = 1.2‰). However, we find that inter-measurement drift between acquisitions limits our accuracy and precision for standardized measurements (i.e., error on reported δ values) to ≈1‰ for the 37 Cl/ 35 Cl measurement. Hence, the benefits of low σ AE may not be fully realized. Third, we demonstrate accuracy via sample/standard comparisons of a methionine sample with 13 C enrichment of ≈20‰ relative to a known standard. Using a sequence of seven 60-min analyses, we recover the following isotope ratios and standardized precisions (i.e., error on reported δ values, denoted propagated acquisition errors, σ PAE ): 33 S/ 32 S (σ PAE = 1.0‰), 34 S/ 32 S (σ PAE = 0.7‰), 15 N/ 14 N (σ PAE = 2.1‰), 2 H/ 1 H (σ PAE = 3.2‰), 13 C/ 12 C (σ PAE = 0.4‰), 18 O/ 16 O (σ PAE = 1.6‰), & 13 C 13 C/ 12 C (σ PAE = 2.8‰) with confirmation of accurate results for the known 13 C/ 12 C and 13 C 13 C/ 12 C enrichments. Together, our results demonstrate the viability of Orbitrap-IRMS for long duration measurements of diverse sample types via an automated reservoir injection system. Inter-measurement stability remains a challenge; we expect our methods to be most applicable to extended measurements of hard-to-observe properties, such as 17 O in organics and clumped isotopologues.

58 GEOSCIENCES↗
DOE OSTI2025

Non-traditional stable isotope measurements using laser desorption Orbitrap mass spectrometry: Implications for planetary missions.

Isotopic fractionation recorded in planetary materials provides insights into physical, chemical, and/or potential biological processes occurring on Solar System bodies. As we enter into the next decades of planetary exploration, the crucial information revealed by isotopic compositions of rocky and icy samples mandates that next generation spaceflight instrumentation possess the capability to measure isotope ratios in situ with sufficient precision/accuracy to distinguish between such processes. Here, in addition to identifying and fingerprinting complex organic materials with high accuracy and ultrahigh mass resolutions, laser desorption Orbitrap™ mass spectrometry (LD-O-MS) has the capability to determine the elemental and isotopic composition of solid planetary materials such as rock, regolith, organics, ice, etc. Here, we use a space-qualified LD-O-MS instrument comprising a 266 nm ultraviolet (UV) laser and Orbitrap mass analyzer ruggedized for planetary applications to investigate the stable isotopic composition of Ti and Zn metal plates. Based on the isotopic analyses, we constrain the performance of the CORALS instrument to sub per mille (‰) level accuracy and precision for Ti and at the per mille (‰) level for Zn. An LD-O-MS instrument is a versatile instrument capable of measuring isotopic composition of a variety of planetary samples and would constitute a critical instrument in the exploration of various planetary bodies including but not limited to Moon, Mars, Enceladus and other ocean or icy worlds, Ceres and other asteroids, and comets, thereby answering several high-priority questions pertaining to the formation and evolution of our Solar System.

58 GEOSCIENCES↗
DOE OSTI2020

Experimental investigation of multi-scale strain-field heterogeneity in rocks

Rocks have inherently heterogeneous micro-structures. The mechanical response of these heterogeneous micro-structures primarily governs the overall macroscopic behavior of intact rocks, as these features cause heterogeneity in the distribution of stress and strain within a rock volume. These attributes have generally been studied through numerical tools, but in this study, the stress-induced multi-scale spatial heterogeneity in the strain-field of three types of rocks was analyzed experimentally. The full-field non-contact strain measurement approach of 2D Digital Image Correlation (DIC) was employed for real-time evaluation of the progression of strain-field heterogeneity across the surface of the rock specimens when subjected to uniaxial compression. The results show that the strain-field heterogeneity is a characteristic of the inherent heterogeneity of the rock specimens, and that it amplifies with the magnitude of stress applied on the rock specimens. The changes in the strain-field heterogeneity with loading were consistent with established microcrack evolution processes in the rock specimens. The statistical fluctuation in the strain-field heterogeneity of the three rocks at varying gauge lengths (GL) was then analyzed for the determination of the Representative Volume Element (RVE) length-scale at increasing load levels. The evolution of the RVE length-scale with loading facilitates an objective way of determination of the stress at which the strain-field heterogeneity increases and begins to influence the overall global response of the material. The RVE length-scale variation with loading showed that at the RVE length-scale reaches the specimen-scale at load levels well below the peak strength.

2D-digital image correlation↗
DOE OSTI2020

A geospatial model for the analysis of time-dependent land subsidence induced by reservoir depletion

Land subsidence due to fluid depletion is an outcome of physical processes operating across a wide range of time and length scales. Although geomechanical models are crucial to simulate reservoir compaction and predict its long-term fate, their use across large regions often bears prohibitive computational costs. To overcome this obstacle, this paper proposes a simplified modelling framework consisting of (i) a near-field numerical solver simulating the coupling between fluid flow and rock deformation with reference to a simplified one-dimensional geometry and (ii) a far-field geospatial algorithm mapping ground settlements across a region through the superposition of poroelastic computations at multiple wells. The model computes the delay between depletion history and reservoir compaction in proximity of a producing well by assuming basal depletion of a fluid-saturated deformable disk, while the Geertsma solution of nucleus of strain is used to extrapolate the impact of such time-varying reservoir compaction around the well. This approach has been used to back-analyze the spatio-temporal progression of subsidence at the Groningen gas field. The results are compared against measurements collected over 50 year of production at 25 benchmark locations scattered over an area of 900 km 2 . It is shown that coupled simulations based on average values of rock compressibility and permeability lead to nonlinear trends of subsidence evolution in good agreement with field measurements, while uncoupled analyses overpredict settlements by more than 70%. Lastly, synthetic forecasts based on different rates of depletion were provided. The results suggest that slower depletion rates lead to lower subsidence at a given time, and that residual subsidence may continue to develop for several decades after interruption of production activities.

58 GEOSCIENCES↗
DOE OSTI2020

A generalized power-law criterion for rocks based on Mohr failure theory

Nonlinearity of Mohr envelopes confirmed by past laboratory test data can be of considerable importance in theoretical interpretation of rock strength behaviors and many practical applications. Current nonlinear failure criteria are introduced with parameters to empirically fit test data, and the physical and mechanical meanings of these fitting parameters are often not known. Predicted strengths by theoretical criteria (e.g., Griffith criterion) generally do not agree well with test data. Considering the nonlinear properties of a Mohr envelope, a nonlinear failure criterion in a generalized power-law form is derived in this paper based on strength magnitude coefficient (a) and failure envelope curvature coefficient (n) which are expressed in terms of uniaxial tensile strength, uniaxial compressive strength, and angle of internal friction. Test data on granite, marble, dolomite, and andesite from published literature were chosen to test the performance of the proposed criterion along with Griffith and Hoek-Brown criteria. Statistical analysis was performed to evaluate the application of different criteria in predicting the strength of rocks. In addition, some Griffith-typed criteria are mathematically shown to be the special cases of the proposed criterion. The proposed and Hoek-Brown criteria have similar power-law forms with a significant difference in the coefficient functions that govern the shear strength magnitude of rock. The different tensile and compressive failure mechanisms are also investigated from the perspective of the relationship between the uniaxial tension Mohr circle and the proposed failure envelope.

58 GEOSCIENCES↗
DOE OSTI2020

Modeling of fluid injection-induced fault reactivation using coupled fluid flow and mechanical interface model

The present study is aimed at developing a numerical model to reproduce coupled hydro-mechanical processes associated with fault reactivation by fluid injection in low permeability rock, as part of the DECOVALEX-2019 project Task B. We proposed a modeling approach for simulating the processes using the TOUGH-FLAC simulator, and modeled a fault reactivation experiment conducted at Mont Terri Rock Laboratory in Switzerland. The first step of the study involved benchmark calculations considering a simplified fault plane and geometry. Fluid flow along a fault was modeled using elements of aperture-sized thickness on the basis of Darcy's law and the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by zero-thickness interface elements in FLAC3D upon which a slip and/or separation is allowed. A methodology to connect a TOUGH2 volume element to a FLAC3D interface element was developed for handling the hydro-mechanical interactions on the fault during fluid injection. Two different fault models for describing the evolutions of hydraulic aperture by elastic fracture opening and failure-induced aperture increase were considered in the benchmark calculations. In the coupling process, the changes in geometrical features and hydrological properties induced by mechanical deformation were continuously updated. The transient responses of the fault and host rock to stepwise pressurization were examined during the simulation. The hydro-mechanical behavior, including the injection flow rate, pressure distribution around the borehole, stress conditions, and displacements in normal and shear directions were monitored in the surrounding rock and along the fault. The results of benchmark calculations suggest that the developed model reasonably represents the hydro-mechanical behavior of a fault and the surrounding rock. This modeling approach was applied to the fault reactivation experiment of the Mont Terri Rock Laboratory. In this interpretive modeling, a parametric study was conducted to examine the effects of input parameters regarding in situ stress and fault properties on the hydro-mechanical responses of the fault to water injection. Then, an optimal parameter set to reproduce the field experiment results was chosen by trial-and-error. The injection flow rate and pressure response during fault reactivation closely matched those obtained at the site, which indicates the capability of the model to appropriately capture the progressive pathway evolution during fault reactivation tests at the site. Finally, the anchor displacements were overestimated by the model, but a fair agreement was obtained in terms of the order of magnitude and the variation tendency.

58 GEOSCIENCES↗
DOE OSTI2020

Prediction of the thermal-hydraulic-mechanical response of a geological repository at large scale and sensitivity analyses

To better understand the Thermal-Hydraulic-Mechanical (THM) response of the rock caused by the high-level radioactive waste (HLW) released heat in the Callovo-Oxfordian formation (COx) in the near-field and far-field areas, a series of coupled THM modelling is performed. This study presents a case study based on the French geological repository for HLW (Cigéo project) to better assess modelling of deep geological repositories (DGR) within the DECOVALEX-2019 framework. In this study, the proposed coupled THM model, implemented in the platform of COMSOL, is validated by comparing the modelled THM results for a point heat source in an infinite rock mass with an analytical solution. It is also validated by comparing its thermal component for the Base Case with accurate results calculated using another numerical method. To shorten the calculation time, the DGR is modelled using six cuboid blocks representing 162 HLW cells and only six central HLW cells of interest modelled in detail. This simplification is validated. This model can easily be used to model not only the near-field THM response of the rock but also the far-field THM response. The influence of the boundary conditions applied on the gallery wall, HLW cell walls, and external surfaces of the model and the influence of different vertical dimensions of the geometry are studied. The sensitivity analyses of the THM parameters of COx on the THM response at different locations are performed. Finally, the influence of using a 2-dimensional (2D) model to represent a 3-dimensional (3D) repository is also investigated.

58 GEOSCIENCES↗
DOE OSTI2020

Determination of permeability for hydrocarbon release due to excavation-induced stress redistribution in rock salt

Due to a stress redistribution after the excavation of an underground tunnel for radioactive waste disposal, an Ed/DZ (excavation disturbed/damaged zone) will be generated in the near field of the opening, resulting in significant changes in the hydraulic and mechanical properties of the rock mass in the zone. Initially more or less randomly distributed hydrocarbons at grain boundaries in rock salt, which sometimes can only be observed with ultraviolet light, can then be mobilised and migrate at a potentially significant rate towards the excavation. Within the international cooperative project DECOVALEX 2019, the migration mechanism of such fluid inclusions in rock salt is being studied intensively. Specifically, a multi-scale modelling strategy has been developed. A macroscale coupled hydro-mechanical modelling of an underground excavation was performed to determine hydraulic and time-dependent deviatoric stress conditions, by taking into account the rock salt creep behaviour. Under the obtained macro-scale constraints, micro-scale modelling of a pathway dilation along halite grain boundaries was performed using different model strategies: a) coupled hydromechanical modelling with a consideration of hydraulic pressure-induced dilatant deformation, b) nonlinear dynamic model taking account of fluid migration, stress-dependent grain boundary wetting and shear-induced dilatancy of salt, and c) phase-field modelling of flow pathway propagation. The permeability increase resulting from the pathway dilation is estimated to be as high as two orders of magnitude. Based on the permeability determined, a series of pressure buildups measured from a borehole with a high hydrocarbon release rate, a total of 430 build-ups within a monitoring time of 938 days, can be simulated with a macro-scale compressible flow model accounting for different zones around the opening.

58 GEOSCIENCES↗