The Solar Wind Electron (SWE) instrument for the Interstellar Mapping and Acceleration Probe Mission
Explore the source record for details and available documents.
58 GEOSCIENCES↗
SEARCH · Search NASA
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.
Characterizing Reactive Transport Behavior in a Three-Dimensional Discrete Fracture Network
While several studies have linked network and in-fracture scale properties to conservative transport behavior in subsurface fractured media, studies on reactive transport cases remain relatively underdeveloped. In this study, we explore the behavior of an irreversible kinetic reaction during the interaction of two solute plumes, one consisting of species A and the other species B. When the plumes converge, these species react kinetically to form a new species C via A+B→kC. This reactive system is studied using a three-dimensional discrete fracture network (DFN) model coupled with reactive Lagrangian particle tracking. We find that the interplay of network topology and chemical properties of the reactive solutes controls reactive transport processes. The network topology drives species A and B together, and the chemical properties dictate whether and how quickly a reaction occurs. Results demonstrate that reactions are most likely to occur in high-velocity fractures that make up the network backbone. The interplay between species’ chemical properties and transport is characterized by a non-dimensional Damköhler (Da) number. We show that the spatial distribution of reactions is sensitive to Da, which subsequently influences late-time tailing behavior in outlet breakthrough time distributions. The results of this study provide initial insights into how an irreversible reaction occurs during transport in a fracture network, using a methodology that can be applied to study reactive transport in a wide range of fractured media environments and contexts.
58 GEOSCIENCES↗
Measuring Gas Transport and Sorption in Large Intact Geologic Specimens via the Piezometric Method
Uptake of noble gases into heterogeneous geologic core samples was measured using a piezometric methodology. In addition to measuring accessible porosity—as with gas pycnometry—by monitoring the rate of pressure decay, this method can also be used to estimate the gas effective diffusivity in the sample. In contrast to previous applications of this method, where milligram quantities of fractured grains are characterized, here approximately kilogram core samples were left intact when tested. In doing so, a more representative sample of the heterogeneous field geology is provided. Additionally, alteration of the pore structure and connectivity during sample preparation is avoided. Furthermore, to scale the piezometric method from milligrams to kilograms, the system was designed to operate at medium vacuum (1 to 100 Pa) to restrict transport in pores less than approximately 60 µm to large Knudsen numbers. To test the system performance, two samples of interest were selected: a rhyolitic welded tuff from Blue Canyon Dome at the Energetic Materials Research and Testing Center and a zeolitized non-welded rhyolitic tuff from the Nevada National Security Site. Three noble gases were utilized in this test series; Argon and xenon as they are of direct interest to nuclear monitoring efforts and helium as it is a weakly adsorbing reference standard. Additionally, mercury intrusion porosimetry measurements were made on subsamples of the core to compare the observed porosity by the two methods and to discuss gas transport rates in the context of the measured pore distribution.
58 GEOSCIENCES↗
Lattice Boltzmann Modeling of Drying of Porous Media Considering Contact Angle Hysteresis
Drying of porous media is governed by a combination of evaporation and movement of the liquid phase within the porous structure. Contact angle hysteresis induced by surface roughness is shown to influence multi-phase flows, such as contact line motion of droplet, phase distribution during drainage and coffee ring formed after droplet drying in constant contact radius mode. However, the influence of contact angle hysteresis on liquid drying in porous media is still an unanswered question. Lattice Boltzmann model (LBM) is an advanced numerical approach increasingly used to study phase change problems including drying. In this paper, based on a geometric formulation scheme to prescribe contact angle, we implement a contact angle hysteresis model within the framework of a two-phase pseudopotential LBM. The capability and accuracy of prescribing and automatically measuring contact angles over a large range are tested and validated by simulating droplets sitting on flat and curved surfaces. Afterward, the proposed contact angle hysteresis model is validated by modeling droplet drying on flat and curved surfaces. Then, drying of two connected capillary tubes is studied, considering the influence of different contact angle hysteresis ranges on drying dynamics. Finally, the model is applied to study drying of a dual-porosity porous medium, where phase distribution and drying rate are compared with and without contact angle hysteresis. The proposed model is shown to be capable of dealing with different contact angle hysteresis ranges accurately and of capturing the physical mechanisms during drying in different porous media including flat and curved geometries.
36 MATERIALS SCIENCE↗
Effect of Shear Displacement and Stress Changes on Fracture Hydraulic Aperture and Flow Anisotropy
Fluid flow in fractures has the potential to drastically change the economic, environmental, and safety risks associated with a subsurface operation. This work focuses on the analysis of experimental shear fracturing and subsequent shear displacement of a shale specimen that contains multiple pre-existing natural fractures. A triaxial direct shear apparatus with concurrent µCT imaging allowed continuous monitoring of sample permeability and changes in fracture geometry during fracture creation, displacement, and changes in stress. Steady-state, low Reynolds number fluid flow simulations were used to determine the changes of each individual fracture’s hydraulic aperture due to displacement and effective stress. The development of anisotropy in the hydraulic aperture of a generated shear fracture was determined for shear displacements of 0.38, 0.74, and 2.04 mm. At high confining stress (30 MPa), shear displacement significantly reduced the hydraulic aperture of the shear fracture to below flow detectable limits and X-ray imaging resolution (23 × 23 × 23 µm 3 ). The shear generated fracture was found to be more sensitive to changes in effective stress compared to the pre-existing natural fractures. Existing models derived from fluid flow equations accounting for stochastic fracture roughness were able to partially fit the simulated flow properties of the pre-existing fractures, but no model was found to fit the flow properties of the newly formed shear fracture. These results are the first report of simultaneous changes in the geometry and flow properties of surrounding, parallel natural fractures as a shear fracture is generated with varying stress and shear displacement conditions.
58 GEOSCIENCES↗
Two-Phase Flow Mechanisms Controlling CO2 Intrusion into Shaly Caprock
Abstract Geologic carbon storage in deep saline aquifers has emerged as a promising technique to mitigate climate change. CO 2 is buoyant at the storage conditions and tends to float over the resident brine jeopardizing long-term containment goals. Therefore, the caprock sealing capacity is of great importance and requires detailed assessment. We perform supercritical CO 2 injection experiments on shaly caprock samples (intact caprock and fault zone) under representative subsurface conditions. We numerically simulate the experiments, satisfactorily reproducing the observed evolution trends. Simulation results highlight the dynamics of CO 2 flow through the specimens with implications to CO 2 leakage risk assessment in field practices. The large injection-induced overpressure drives CO 2 in free phase into the caprock specimens. However, the relative permeability increase following the drainage path is insufficient to provoke an effective advancement of the free-phase CO 2 . As a result, the bulk CO 2 front becomes almost immobile. This implies that the caprock sealing capacity is unlikely to be compromised by a rapid capillary breakthrough and the injected CO 2 does not penetrate deep into the caprock. In the long term, the intrinsically slow molecular diffusion appears to dominate the migration of CO 2 dissolved into brine. Nonetheless, the inherently tortuous nature of shaly caprock further holds back the diffusive flow, favoring safe underground storage of CO 2 over geological time scales.
58 GEOSCIENCES↗
Comparative Analysis of Imaging and Measurements of Micrometer-Scale Fracture Aperture Fields Within a Heterogeneous Rock Using PET and X-ray CT
Knowledge of the spatial distribution of fracture apertures is essential for reliable characterization of flow and transport processes in fractured systems and for better understanding of physicochemical matrix–fracture interactions. Here, we propose and test two image-based methods, thereby extending the current experimental capabilities to characterize aperture size distribution in structurally heterogeneous geologic porous media noninvasively. The first approach utilizes an inversion method based on the dataset acquired from positron emission tomography (PET) and the second approach considers an extension of the classic missing attenuation technique that relies on clinical X-ray computed tomography (X-ray CT). Independent sets of imaging experiments are conducted on a fractured basalt core with heterogeneous matrix properties and aperture distributions to compare the two methodologies. A repeat of each experiment is conducted to verify the proposed workflows. The performance of these two imaging techniques is systematically evaluated through the analysis of signal-to-noise ratio, minimum fracture size detectability, and measurement errors. While both approaches provide a reliable estimation of fracture aperture distributions, PET yields a signal-to-noise ratio that is substantially higher than the corresponding X-ray CT measurements. Furthermore, uncertainties of the aperture values for PET are considerably lower ($\bar {\sigma}_{\text {d}} = 15\%$) compared to those obtained from X-ray CT ($\bar {\sigma}_{\text {d}} = 29\%$), allowing for the detection of minimum aperture sizes of 20 $\mu$m with 70% confidence level. Finally, these approaches provide key experimental tools for better understanding dynamic hydromechanical fracture properties in geologic systems.
58 GEOSCIENCES↗
Explosive Byproduct Gas Transport Through Sorptive Geomedia
Current underground nuclear explosion (UNE) detection strategies rely heavily on atmospheric noble gas sampling of radioxenon. However, discriminating nuclear weapons testing programs from civilian sources is difficult due to highly variable atmospheric radioxenon backgrounds and processes affecting subsurface transport of parent radionuclides. Here, we aim to study the transport of gases produced by subsurface explosions as novel stable signatures for underground nuclear explosion (UNE) monitoring. These gases may be produced in large quantities with distinct molecular ratios, which will be impacted by subsurface transport processes. To demonstrate how ratios of gases produced by explosions can change during transport in geomaterials, we conducted laboratory benchtop experiments on the transport of carbon dioxide (CO 2 ) and hydrogen (H 2 ) gases through variably saturated zeolitic tuff, which is abundant at the historic US testing site. We observed that zeolitic tuff sorbs substantial quantities of CO 2 while allowing H 2 to transport more freely, leading to changes in the molecular ratios of the two gases along the transport pathway. Gas uptake in the dry zeolitic tuff core was 72.3% for CO 2 , compared with 53.4% for xenon and 7.6% for H 2 . The presence of 20% water saturation disrupted the CO 2 sorption process, though to a lesser extent than observed for noble gases, with a 36.7% drop in xenon sorption compared with a 21.9% drop for CO 2 . These results represent the first observations of zeolite sorption altering explosive gas ratios during transport through geomedia relevant to nuclear proliferation monitoring.
54 ENVIRONMENTAL SCIENCES↗
Effect of particle characteristics on the evolution of particle size, particle morphology, and fabric of sands loaded under uniaxial compression
This work presents the results and analyses of uniaxial compression experiments performed on three silica sands. The sands have comparable particle-size distributions, but their particles differ in morphology and strength. Cylindrical samples of the three sands were compressed in a loading device placed inside an X-ray microscope (XRM) and scanned at multiple stress levels during uniaxial compression. 3D tomography data of the samples obtained from the XRM at different stress levels were then analyzed to obtain the distributions of particle size, particle morphology, and interparticle contact normals within the sample. Results indicate that: (1) the compressibility of the sands loaded under uniaxial compression is closely tied to particle morphology and strength and (2) the anisotropy in the orientations of interparticle contact normals generally increases with axial stress; however, this increase is limited by the occurrence of particle crushing in the sample.
3D X-ray computed tomography↗
Thermophysical properties of the Ghareb formation relevant for nuclear waste disposal
The Ghareb formation, a shallowly buried porous chalk in Israel, is currently a candidate for nuclear waste disposal. The potential repository is somewhat unique for its host rock and emplacement in shallow (500 m) engineered large diameter boreholes. Herein, the thermal properties of the Ghareb are determined to support design and performance assessment; the relevant properties measured are thermal conductivity, specific heat capacity, thermal diffusivity, and thermal expansion coefficient along with their relationship with varying temperature. For the temperature range of 40 to 275 °C, the thermal conductivity ranges from 0.30 to 1.10 W/m·K, the thermal diffusivity ranges from 0.20 to 0.72 mm 2 /s, and volumetric heat capacity ranges from 0.86 to 2.00 MJ/m 3 ·K. Thermal strain measurements were used to estimate the linear thermal expansion coefficient to be 6·10 –4 –9·10 –2 °C -1 from 40 to 300 °C. These measured properties were used in a thermomechanical model to estimate near-field stresses an hour and 10 years after waste emplacement; the borehole was found to be stable. Thermal loading after 10 years was predicted to elevate local pore pressures by 1–1.5 MPa. The laboratory measurements coupled with analyses are the first attempts at performance assessment characterization for this first of its kind potential repository setting with this chosen host rock.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
Seismic identification of geothermal prospecting in Harrat Rahat, Northern Arabian Shield
Harrat Rahat volcanic field lies along the Makkah-Madinah-Nafud active volcanic line. It has extensive linear weak zones of wrench faults that could provide ideal conduits for upwelling magma and hydrothermal fluids as well. Here, passive seismic data recorded by seismic stations installed within and around Harrat Rahat has been processed for the detection of microseismic events. Moreover, ambient noise was also gathered and analyzed to construct a subsurface velocity model. This model indicates low velocity zone in Harrat Rahat, and two areas with coincident basement faults and possibly with enhanced subsurface temperatures were identified. It is concluded that Harrat Rahat has high geothermal potential and promising for geothermal renewable power.
58 GEOSCIENCES↗
Seismic source characterization of the Arabian Peninsula and Zagros Mountains from regional moment tensor and coda envelopes
Reliable estimates of moment magnitude and source mechanism for seismic events in the Middle East can be challenging due to a small number of openly available stations, the complex tectonic setting, and regions of high attenuation. Access to high-quality waveform data from well-calibrated regional seismic stations is fundamental in producing robust and stable estimates of earthquake source parameters, particularly when measurements of absolute ground motion amplitudes are required. Earthquake source mechanism and moment magnitude are invaluable information in the assessment of seismic hazard, plate motions, and the characterization of faults and regional stress field. The expansion of the Kingdom of Saudi Arabia’s national seismic network in the past several decades provides an opportunity to develop the capabilities of routine focal mechanism and moment magnitude estimations. In this study, we use time-domain full waveform moment tensor inversion and coda envelope–derived amplitude measurements to solve for earthquake source mechanism, moment magnitudes, and their source type. Furthermore, we compared the moment magnitudes calculated from the two methods and publicly available earthquake catalogs and discuss the implications of the obtained source parameters.
58 GEOSCIENCES↗
Geothermal potential of Harrat Rahat, Northern Arabian Shield: geological constraints
Arabian Shield offers a strong potential for geothermal energy. However, only a limited amount of exploration effort has taken place. Harrat Rahat is a recent volcanic area located along the Makkah-Madinah-Nafud volcanic line in western Saudi Arabia. This study is oriented to identify the geothermal potential of Harrat Rahat. The geologic data was collected and processed. Here, we focus on evaluating the potential sources of subsurface heat and groundwater. Subsurface heat is based on the age of erupted volcanoes and the expected volume of subsurface magma chambers. An analogy with an existing geothermal field in a similar setting as Mexico suggests that the faults and the coincident seismicity may be associated with higher permeability. It is concluded that Harrat Rahat is characterized by a good geothermal potential for future geothermal energy source production.
58 GEOSCIENCES↗
A novel framework for elastoplastic behaviour of anisotropic solids
A general framework for developing nonlinear hyperelastic/plastic constitutive laws for anisotropic solids experiencing large strains and strain rates has been developed. The proposed framework does not rely on the “a priori” known strain energy function, but instead introduces a physical decomposition of the material element into seven physically independent stress bearing mechanisms, each of which has a constitutive law in terms of internal moments described by a scalar function of a single variable. Furthermore, the model has been encoded into a combined finite-discrete element method and tested against static geomechanical test data. The numerical validation experiments show the model can reproduce plastic anisotropic behavior in both biaxial and uniaxial loading of a geomaterial.
58 GEOSCIENCES↗
HOSS: an implementation of the combined finite-discrete element method
Nearly thirty years since its inception, the combined finite-discrete element method (FDEM) has made remarkable strides in becoming a mainstream analysis tool within the field of Computational Mechanics. FDEM was developed to effectively “bridge the gap” between two disparate Computational Mechanics approaches known as the finite and discrete element meth-ods. At Los Alamos National Laboratory (LANL) researchers developed the Hybrid Optimization Software Suite (HOSS) as a hybrid multi-physics platform, based on FDEM, for the simulation of solid material behavior complemented with the latest technological enhancements for full fluid–solid interaction. Furthermore, in HOSS, several newly developed FDEM algorithms have been implemented that yield more accurate material deformation formulations, inter-particle interaction solvers, and fracture and fragmentation solutions. Additionally, an explicit computational fluid dynamics solver and a novel fluid–solid interaction algorithms have been fully integrated (as opposed to coupled) into the HOSS’ solid mechanical solver, allowing for the study of an even wider range of problems. Advancements such as this are leading HOSS to become a tool of choice for multi-physics problems. Finally, HOSS has been successfully applied by a myriad of researchers for analysis in rock mechanics, oil and gas industries, engineering application (structural, mechanical and biomedical engineering), mining, blast loading, high velocity impact, as well as seismic and acoustic analysis. This paper intends to summarize the latest development and application efforts for HOSS.
58 GEOSCIENCES↗
Integration of Large-Scale Electrical Imaging into Geological Framework Development and Refinement
Geologic framework models (GFMs) are critical to the construction of reliable simulation models of groundwater flow and contaminant transport. To support GFM development, direct information (e.g., core samples, fluid samples, hydraulic testing) tends to be sparse and separated by large distances relative to the spatial scales of aquifer heterogeneity. There are additional challenges associated with highly contaminated legacy waste sites, where drilling is particularly costly, and invasive sampling requires specialized handling and disposal of hazardous materials. At these sites in particular, non-invasive geophysical imaging can play an important role in filling spatial gaps between boreholes and reducing characterization costs by optimizing and minimizing the number of necessary boreholes. Here this paper presents a case study demonstrating the use of large-scale (> 30 km 2 ) electrical mapping to identify hydrostratigraphy and potential paleochannels at the Hanford Site, located in Washington State, USA. In two field campaigns, over 36 line-kilometers of electrical resistivity tomography (ERT) data were collected along 14 transects. ERT surveys were sited and performed to image critical aspects (e.g., paleochannels, stratigraphic contacts) of the subsurface, demonstrating a general workflow for integrating ERT with GFM development. Inconsistencies between the GFM and ERT were catalogued to provide a basis for future site characterization using complementary geophysical methods and (or) direct sampling.
58 GEOSCIENCES↗
A spatial-statistical investigation of surface expressions associated with cyclic steaming in the Midway-Sunset Oil Field, California
In the Midway Sunset Oil Field in Central California, operators inject steam into the shallow diatomite formation to enhance heavy oil recovery through imbibition, wettability alteration, and viscosity reduction, among other mechanisms. The injected steam, however, does not always remain in the reservoir or return through the wells. In two zones in the study area, the steam comes out at the surface, creating sinkholes, seeps, and steam outlets. These phenomena, called “surface expressions,” pose safety and environmental hazards. Even though these surface expressions are a widespread problem in Central California, they are not well documented and understood. Possible causes of the surface expressions include: high injection pressure, structurally controlled flow patterns, leakage of steam through old improperly abandoned wells, high injection volumes, or flow along naturally occurring faults, among other possible factors. This work examines attributes of the zones with surface expressions in order to determine factors that may contribute to their occurrence. Spatial statistical analysis using logistic regression, random forests, and classification trees is used to explore the relationship between the surface expressions and geological and production-related attributes. The results point to a significant spatial correlation between the surface expressions and two predictors: concentration of plugged wells and geologic seal thickness. The results guide follow-up studies to further investigate the role of well abandonment and seal thickness in the occurrence of surface expressions.
02 PETROLEUM↗
Borehole breakout modeling in arkose and granite rocks
The existence of a deep borehole in the Earth’s crust disturbs the local stresses and creates a stress concentration that may result in breakout and damage to the borehole. Maintaining wellbore integrity mitigates environmental impacts such as groundwater contamination, gas leakage to the atmosphere, and fluid spills and seepage at the surface. In this paper, the stability of deep boreholes (5 km) is examined by laboratory experiments and numerical models in the context of nuclear waste disposal in Israel. Two rock types in southern Israel are considered: the crystalline basement (granite) and the Zenifim Formation (arkose). A series of room-temperature triaxial rock deformation experiments were conducted at different confining pressures. This mechanical characterization was then used to parameterize the elastic properties and damage behavior of the rocks. This facilitated modeling the stability of the deep boreholes by two different formulations of damage rheology: a dynamic-oriented formulation used to model deformation immediately after the creation of the open hole and a quasi-static formulation used to model longer stress corrosion regime. The calibrated modeling results indicate greater stability with Zenifim arkose than the crystalline granite for deep borehole conditions despite the granite having a greater triaxial compressive strength. Dissipation associated with dilation and porous compaction in the arkose during deformation plays a significant stabilizing role in the borehole compared to crystalline rocks. These results suggest that common strength-based borehole stability assessment may lead to inaccurate predictions. Three-dimensional modeling of bottom-hole stress conditions and the effects of transient borehole geometry show conventional two-dimensional analysis may not be conservative when predicting borehole damage.
58 GEOSCIENCES↗