Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ROCK”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Evolution of Pore Structure and Permeability of Rocks Under Hydrothermal Conditions (Final Report)

The physical and transport properties of porous rocks can be altered by a variety of diagenetic, metamorphic, and tectonic processes, and the changes that result are of critical importance to such industrial applications as resource recovery, carbon dioxide sequestration, and waste isolation in geologic formations. The associated interrelationships between rocks, pore fluids and deformation are also key to understanding many natural processes, including dynamic metamorphism, fault mechanics, fault stability and pressure solution creep. In this project, we investigated the changes of permeability and pore geometry owing to inelastic deformation by solution-transfer, brittle fracturing and dislocation creep. In particular, we studied the coupling between pore fluids and deformation in fluid-filled quartz, calcite, mudstones and ultramafic rocks and examine the effects of coupled mechanical and chemical processes on the evolution of porosity and permeability under hydrothermal conditions. The investigations used a combination of techniques, that included triaxial laboratory experiments; uniaxial compressive loading with in situ observations of microstructure; permeability under hydrostatic compression; observations of microstructure using optical and electron microscopes, and micro-CT imaging; and numerical calculations. Laboratory experiments were designed to provide mechanical and transport data under conditions that isolate the particular mechanisms responsible for the changes. The data obtained were used to quantify changes in surface roughness, porosity, pore dimensions and their spatial fluctuations. The results of the experiments and data from image analyses were compared to the results of network, finite-difference and other numerical models to verify the validity of experimentally established relations between permeability and other rock properties. A bibliography for the period from 2015-2019 is appended below. New results obtained in the final year and the unfunded extension period are detailed below.

58 GEOSCIENCES↗
DOE OSTI2016

Discrete element modeling of indentation tests to investigate mechanisms of CO 2 -related chemomechanical rock alteration

During CO 2 injection into geological formations, petrophysical and geomechanical properties of host formations can be altered due to mineral dissolution and precipitation. Field and laboratory results have shown that sandstone and siltstone can be altered by CO 2 -water mixtures, but few quantitative studies have been performed to fully investigate underlying mechanisms. Based on the hypothesis that CO 2 -water mixtures alter the integrity of rock structure by attacking cements rather than grains, we attempt to explain the degradation of cementation due to long-term contact with CO 2 and water and mechanisms for changes in rock mechanical properties. Many sandstones, including calcite-cemented quartzitic sandstone, chlorite-cemented quartzitic sandstone, and hematite-cemented quartzitic sandstone, contain interparticle cements that are more readily affected by CO 2 -water mixtures than grains. A model that couples the discrete element method and the bonded-particle model is used to perform simulations of indentation tests on synthetic rocks with crystal and random packings. The model is verified against the analytical cavity expansion model and validated against laboratory indentation tests on Entrada sandstone with and without CO 2 alteration. Sensitivity analysis is performed for cementation microscopic parameters including stiffness, size, axial, and shear strength. The simulation results indicate that the CO 2 -related degradation of mechanical properties in bleached Entrada sandstone can be attributed to the reduction of cement size rather than cement strength. Finally, our study indicates that it is possible to describe the CO 2 -related rock alteration through particle-scale mechanisms.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Extended Barton–Bandis model for rock joints under cyclic loading: Formulation and implicit algorithm

In this paper, the Barton–Bandis model for rock joints is extended to cyclic loading conditions, without any new material parameter. Also developed herein is an implicit solution algorithm for the extended Barton–Bandis model, which can also be used for the original Barton–Bandis model for which an implicit algorithm has been unavailable. To this end, we first cast the Barton–Bandis model into an incremental elasto-plastic framework, deriving an expression for the elastic shear stiffness being consistent with the original model formulation. We then extend the model formulation to cyclic loading conditions, incorporating the dependence of shear stress and dilation on the joint position and the shearing direction. The extension is achieved by introducing a few state-dependent variables which can be calculated with the existing material parameters. For robust and accurate utilization of the model, we also develop an implicit algorithm based on return mapping, which is unconditionally stable and guarantees the satisfaction of the strength criterion. We verify that the proposed model formulation and algorithm produce virtually the same results as the original Barton–Bandis model under monotonic shearing conditions. Here, we then validate the extended Barton–Bandis model against experimental data on natural rock joints under cycling loading conditions. The present work thus enables the Barton–Bandis model, which has been exceptionally popular in research and practice, to be applicable to a wider range of problems in rock mechanics and rock engineering.

58 GEOSCIENCES↗
DOE OSTI2021

Multi-scale Coupled Processes Modeling of Fractures as Porous, Interfacial and Granular Systems from Rock Images with the Numerical Manifold Method

Abstract The greatest challenges of rigorously modeling coupled hydro-mechanical processes in fractured rocks at different scales are associated with computational geometry. In addition, selections of continuous or discontinuous models, physical laws, and coupling priorities at different scales based on different geometric features determine the applicability of a numerical model for a certain type of problem. In this study, we present our multi-scale modeling capabilities that have been developed based on the numerical manifold method for analyzing coupled hydro-mechanical processes in fractured rocks. Based on their geometric features, the fractures are modeled as continua—finite-thickness porous zones, and discontinua—discontinuous interfaces and microscale asperities and granular systems. Different governing equations, physical laws, coupling priorities, and approaches for addressing fracture intersections and shearing are then applied to describe these. We applied these models to simulate coupled processes in fractured rocks using realistic geometry obtained from rock images at different scales. We first calculated shearing of a single fracture with different models and demonstrated the impacts of asperities on shearing. We then applied the continuous and discontinuous models to simulate a network of rough fractures, demonstrating that contact dynamics contribute significantly to the geometric, multi-physical evolution of systems where rough fractures are not mineral filled. For a discrete fracture network, our coupled processes modeling demonstrates that shearing of the discrete fractures can have a major impact on stress and pore pressure distribution. Lastly, we applied the discontinuous granular model to simulate evolution of a complex granular system with a deformation band, demonstrating that the deformation band can dominate contact dynamics, the structural and the stress evolution of the granular system.

58 GEOSCIENCES↗
DOE OSTI2021

A Time-Dependent Directional Damage Theory for Brittle Rocks Considering the Kinetics of Microcrack Growth

In this work, a novel microcrack damage theory for describing the time-dependent behavior of brittle rocks is proposed. Instead of using scalar or tensorial variables to approximate the effect of microcracks, the concept of directional damage is introduced by upscaling the directional distribution density of microcracks through the noninteracting homogenization scheme. The contribution of an individual crack on global strain is first derived from the crack-opening displacements including sliding, closure, and dilation. The model is then cast in the free energy form, and the range of parameters ensuring compliance with the second law of thermodynamics is derived rigorously. Considering that the subcritical propagation of microcracks is the dominant mechanism behind the time-dependent deformation and failure of brittle rocks, the driving force for directional damage is defined as the corresponding energy release rate, and the effect of time is introduced through a damage evolution law inspired from the kinetics of subcritical crack growth. The model is applied to simulate the behavior of basalt rock under various loading conditions. In addition to predicting the stress–strain–time response, the model offers enhanced resolution in representing the anisotropic characteristics of microcrack-induced damage in brittle rocks and their time-dependent evolutions.

58 GEOSCIENCES↗
DOE OSTI2025

Geomechanical characterization of Rock Valley carbonates

In Rock Valley, Nevada, the USA, past earthquake sequences and shallow faulting in Paleozoic carbonates remain poorly understood. Carbonates are typically more ductile rocks that do not experience large stress drops or fracture coalescence, which contradicts previous observations in the region. As part of the Source Physics Experiment, an effort has been made to experimentally characterize the petrophysical and geomechanical properties of carbonates from Rock Valley. Well core and outcrop samples were used to determine the difference between shallow-buried Tertiary limestones and deeply buried Paleozoic limestones and dolostones. The carbonates possessed porosities between 3 and 9%, with V P and V S in the Tertiary carbonates around 2700 and 1800 m/s, respectively, and 6000 and 3100 m/s, respectively, in the Paleozoic carbonates. Microstructural characterization revealed that the Paleozoic carbonates contained significantly more pre-existing damage and alteration than the Tertiary carbonates, though the deformation varies from localized to diffuse. Unconfined Brazilian tests and triaxial tests with confining pressures between 0 and 50 MPa showed that, although samples all experienced failure, the dolostones typically failed at greater stresses and possessed greater E/ν ratios than the limestones. Furthermore, the Hoek–Brown failure criterion was used to construct failure envelopes with the compressive and tensile failure tests. Velocity, density, and porosity measurements were used to construct a hypothetical velocity-depth profile and compare the results with the theoretical petrophysical measurements. Using brittleness analysis, the likely failure conditions were determined to be low-porosity dolomitic rock for fault nucleation in the Paleozoic basement.

Brittleness↗
DOE OSTI2024

Validating the rhenium proxy for rock organic carbon oxidation using weathering profiles

Chemical weathering over geological timescales acts as a source or sink of atmospheric carbon dioxide (CO 2 ), while influencing long-term redox cycling and atmospheric oxygen (O 2 ) at Earth's surface. There is a growing recognition that the oxidative weathering of rock organic carbon (OC petro ) can release more CO 2 than is locally drawn down by silicate weathering, and may vary due to changes in erosion and climate. The element rhenium (Re) has emerged as a proxy to track the oxidative weathering of OC petro , yet uncertainties in its application remain namely that we lack a systematic assessment of the comparative mobility of Re and OC petro during sedimentary rock weathering. Here we measure Re and OC petro loss across gradients in rock weathering at 9 global sites, spanning a range of initial OC petro values from ~0.2 % to >10 %. We use titanium to account for volume changes during weathering and assess Re and OC petro loss alongside major elements that reflect silicate (Na, Mg), carbonate (Ca, Mg) and sulfide (S) weathering. Across the dataset, Re loss is correlated with OC petro loss but not with loss of any other major element. Further, across the weathering profiles, the average molar ratio of OC petro to Re loss was 0.84 ± 0.15, with 8 out of 9 sites having a ratio >0.74. At one site (Marcellus Shale), the average ratio was lower at 0.58 ± 0.11. The excess loss of Re matches expectations that, typically, between ~0 and 20 % of the Re liberated by sedimentary rock weathering derives from silicate or sulfide phases, while some OC petro may be physically or chemically protected from weathering. Overall, our measurements provide validation for the Re proxy of OC petro oxidation and allow future work to further improve our knowledge of regional and global-scale rates of this important source of CO 2 in the geochemical carbon cycle.

58 GEOSCIENCES↗
DOE OSTI2018

Not all moissanites are created equal: New constraints on moissanite from metamorphic rocks of Bulgaria

Here, terrestrial moissanite (SiC) is widely reported as an ultra-high pressure mineral occurring in kimberlites, diamonds and ultramafic/mafic rocks of mantle origin. However, the conditions of crystallization remain largely unknown. Moreover, dozens of SiC occurrences have been reported from continental crust sources such as granitoids, andesite–dacite volcanic rocks and their breccia, metasomatic and metamorphic rocks, and even limestones. The validity of many of these reports is still debated primarily due to possible contaminations from the widespread use of synthetic SiC abrasives in samples preparation. Indeed, reports of well-documented in-situ occurrences of moissanite in association with co-existing minerals are still scarce. The only condition of moissanite formation that is agreed upon is that extremely reducing media are required (e.g. 4.5–6 log units below the iron-wustite buffer). Here, we report the new occurrence of moissanite that was found in-situ within the garnet–staurolite–mica schists of Topolovgrad metamorphic group of Triassic age in Southern Bulgaria. The 10–300 μm moissanite crystals are situated within 0.1–1.2 mm isolated clusters, filled with amorphous carbon and nanocrystalline graphite. Most of moissanite crystals are 15R (rhombohedral) and 6H (hexagonal) polytypes, and one prismatic crystal, found within them, exhibits unusual concentric polytypical zoning with core (15R), intermediate zone (6H) and rim (3C-cubic). Experimental data show that this type of polytypical zonation is likely due to a decrease in temperature (or/and pressure?) and changes in Si/C ratio. Indeed, amphibolite facies metamorphism (500–580°C – garnet–staurolite zone) followed by a subsequent cooling during the retrograde stage of green schist facies metamorphism (~400–500°C) could have provided a change in temperature. The SiC containing clusters exhibit evidence that they are pre-metamorphic, and we hypothesize that their protolith was a “black shale” material likely rich in carbon, hydrocarbon and terrigenous silica. The latter served as a source of isolated chemically-reduced media, which is required for SiC formation. Other concepts to explain moissanite occurrences in metasedimentary rocks are also discussed. Importantly, our findings show that the formation conditions of moissanite are likely more variable than previously recognized.

15R-SiC↗
DOE OSTI2018

Catalytic generation of methane at 60–100 °C and 0.1–300 MPa from source rocks containing kerogen Types I, II, and III

Low temperature (60 and 100 °C) and long-term (6 months to 5 years) heating of pre-evacuated and sterilized shales and coals containing kerogen Types I (Mahogany Shale), II (Mowry Shale and New Albany Shale), and III (Springfield Coal and Wilcox Lignite) with low initial maturities (vitrinite reflectance Ro 0.39–0.62%) demonstrates that catalytically generated hydrocarbons may explain the occurrence of some non-biogenic natural gas accumulations where insufficient thermal maturity contradicts the conventional thermal cracking paradigm. Extrapolation of the observed rate of catalytic methanogenesis in the laboratory suggests that significant amounts of sedimentary organic carbon can be converted to relatively dry natural gas over tens of thousands of years in sedimentary basins at temperatures as low as 60 °C. Here, our laboratory experiments utilized source rock (shale and coal) chips sealed in gold and Pyrex® glass tubes in the presence of hydrogen-isotopically contrasting waters. Parallel heating experiments applied hydrostatic pressures from 0.1 to 300 MPa. Control experiments constrained the influence of pre-existing and residual methane in closed pores of rock chips that was unrelated to newly generated methane. This study’s experimental methane yields at 60 and 100 °C are 5–11 orders of magnitude higher than the theoretically predicted yields from kinetic models of thermogenic methane generation, which strongly suggests a contribution of catalytic methanogenesis. Higher temperature, longer heating time, and lower hydrostatic pressure enhanced catalytic methanogenesis. No clear relationships were observed between kerogen type or total organic carbon content and methane yields via catalysis. Catalytic methanogenesis was strongest in Mowry Shale where methane yields at 60 °C amounted to ~2.5 μmol per gram of organic carbon after one year of hydrous heating at ambient pressure. In stark contrast to the earlier findings of hydrogen isotopic exchange between water and thermogenic methane in hydrous pyrolysis experiments above 300 °C, the hydrogen isotopic composition of added water exerted limited influence on the δ2H value of methane generated catalytically at low temperatures. We hypothesize that the catalytic sites responsible for methanogenesis are located in hydrophobic microenvironments with limited access to water. The δ13CCH4 values of methane generated catalytically at 60–100 °C range from ~–57.6 to –41.4‰ and are thus similar to typical thermogenic methane (δ13CCH4 >–50‰) and microbially generated methane (<–55‰). Future studies need to evaluate the possibility that clumped isotope characteristics of catalytically generated methane can diagnose the low-temperature regime of catalytic methanogenesis. Furthermore, testing of freshly cored anoxic rocks is needed to determine whether the use of archived, oxygen-exposed rocks in geochemical maturation/catalysis studies introduces artifacts in experimental hydrocarbon yields.

58 GEOSCIENCES↗
DOE OSTI2021

Mineralogy-dependent sulfide oxidation via polysulfide and thiosulfate pathways during weathering of mixed-sulfide bearing mine waste rock

Tracking the S oxidation pathway in sulfide-bearing mine waste-rock piles is complicated by variations in water content, O 2 and Fe 3+ concentrations, microbial diversity, mineralogy, the occurrence of a variety of S species associated with incomplete oxidation, and nonlinear coupling between physicochemical processes. Synchrotron-based S K-edge X-ray absorption near edge structure (XANES) spectroscopy facilitates the identification and quantification of S species (i.e., S 2– , $S_{2}O_{2}^{2–}$, S 8 , $S_{2}O_{3}^{2–}$, and $S_{2}O_{4}^{2–}$) in a weathered mine waste-rock pile containing pyrite, sphalerite, galena, pyrrhotite, and chalcopyrite. Mineralogy-dependent polysulfide and thiosulfate pathways both affect the S oxidation within the waste-rock pile. Currently, the polysulfide pathway, with S 8 as the intermediate S species, generates high concentrations of dissolved metals (from the rapid oxidation of monosulfides including sphalerite, galena, pyrrhotite, and chalcopyrite). As these monosulfides are depleted, the thiosulfate pathway of pyrite oxidation, including intermediate oxidation products (i.e., S 8 and $S_{2}O_{3}^{2–}$), will become the dominant oxidation pathway, resulting in the potential for generation of additional acidic drainage. This information highlights the importance of identifying the sulfide-mineral oxidation reaction pathways when attempting to estimate acid generation potential, and when developing remediation strategies for the storage and management of sulfide-bearing waste rock.

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↗
DOE OSTI2019

Approaches to multi-scale analyses of mechanically and thermally-driven migration of fluid inclusions in salt rocks

Fluid inclusions are found within mineral crystals or along grain boundaries in many sedimentary rocks, notably in evaporite formations, and can migrate along a thermal or hydro-mechanical gradient. Shale and salt rocks have been considered potential host rocks for radioactive waste disposal, due to their low permeability. Previously stagnant inclusions may become mobilised by a perturbation of the in situ state by a geotechnical installation or the emplacement of heat-generating waste. Here, the migration of fluid inclusions can thus have important impacts on the long-term performance of a geologic repository for high-level radioactive waste disposal. As a part of the international research project DECOVALEX-2019, two aspects of fluid inclusion migration in rock salt are currently investigated under different boundary conditions: a) altered hydro-mechanical conditions as a consequence of tunnel excavation or borehole drilling and b) coupled thermo-hydro-mechanical-chemical conditions during the heating period of the post-closure phase of a repository.

58 GEOSCIENCES↗
DOE OSTI2022

Characterization of Acoustic Emissions From Analogue Rocks Using Sparse Regression‐DMDc

Abstract Moisture loss in rock is known to generate acoustic emissions (AE). Phenomena that result in AE during drying are related to the movement of fluids through the pores and induced‐cracks that arise from differential mineral shrinkage, especially in clay‐bearing rock. AE from the movement of fluids occurs from the reconfiguration of fluid interfaces during drying, while AE from mineral shrinkage involves the debonding within or between minerals. Here, analogue rock samples were used to examine the differences in the AE signatures when one or both AE source‐types are present. An unsupervised sparse regression model, Dynamic Mode Decomposition with control, that extends Dynamic Mode Decomposition is used to characterize the AE signals recorded during the drying of porous analogue rock samples fabricated with ordinary Portland cement, with and without clay. This method can effectively and accurately reconstruct acoustic signals emitted from samples that only experience moisture loss without cracking. However, the method struggles to reconstruct signals from samples with intricate crack networks that formed during drying because AE generating mechanisms can emit contemporaneously, and the resulting waves propagate through drying‐induced cracks that can lead to multiple internal reflections. Thus, the differential reconstruction accuracy of time series generated by different underlying physical processes provides a robust filter for reducing large data catalogs. In general, both dynamics and sparse initiating events are learned directly from data and this method exposes a data hierarchy based on the complexity of the intrinsic dynamics.

58 GEOSCIENCES↗
DOE OSTI2022

Structural Characterization of a Toppling Rock Slab From Array-Based Ambient Vibration Measurements and Numerical Modal Analysis

Accurate assessments of the internal structure and boundary conditions of unstable rock slopes are imperative for evaluating landslide hazard scenarios. However, instability characterization at depth remains challenging and is often limited by costly or invasive subsurface investigations. Here, we develop a new approach coupling array-based ambient vibration modal analysis and numerical modeling to improve structural characterization of rock slope instabilities at depth. We used ambient noise cross-correlation on 4 hr of seismic data recorded by an array of 30 nodal geophones at a 500-m-long toppling rock slab in Utah, USA to identify modal frequencies between 0.8 and 3.5 Hz and derive modal displacements. We show that transverse and longitudinal bending modes span the length of the instability, indicating an interconnected slab. Further, statistical comparison of field results with outputs from >1,000 finite element models with varying boundary conditions showed that the instability depth varies between 40–70 and 10–20 m in the middle and lateral regions, respectively. Our approach yields new information on the structural conditions of rock cliff and column instabilities at depth, which is not easily obtained by other means but is imperative for change detection monitoring and improved hazard assessments.

58 GEOSCIENCES↗
DOE OSTI2026

Stimulated geologic hydrogen: from mechanistic control to engineered rock transformation

Geologic hydrogen (GeoH 2 ) generated from subsurface iron-rich rock–water reaction (i.e., serpentinization) is emerging as a promising candidate for the next primary energy source. Yet, accelerating GeoH 2 production from geological to human timescales via enhanced serpentinization remains a formidable scientific and technical challenge. Here, in this Review, we decipher the mechanistic control and explore strategies for accelerating in situ, engineered iron-rich rock transformation into carbon-free GeoH 2 by orders of magnitude. Serpentinization rate is hindered by low porosity and permeability of source rocks, suboptimal temperatures, unfavorable water chemistry, inefficient Fe 2+ -to-Fe 3+ conversion, thermodynamic constraints, and low reactive surface area. While closed-system experiments provide valuable mechanistic insights, open-system conditions with fluid circulation are more crucial for economically viable GeoH 2 production. We assess stimulation techniques from enhanced hydrocarbon and geothermal recovery as tools to be adopted or adapted for increasing reactive surface areas for stimulated GeoH2 production. We estimate that 7.40 × 10 5 to 1.73 × 10 6 million metric tons (Mt) of hydrogen could be engineered over 20 to 50 years from about 10% iron-rich rocks within 10 km depth of continental crust. Enabling GeoH 2 as a viable energy source requires not only advancing scientific frontiers but also forming a global GeoH 2 research network and innovation ecosystem to address the critical scientific, technical, societal, economic, and policy challenges.

08 HYDROGEN↗
DOE OSTI2013

Ion beam analyses of radionuclide migration in heterogeneous rocks

The migration of radionuclides (RN) in the environment is a topic of general interest, for its implications on public health, and it is an issue for the long-term safety studies of deep geological repositories (DGR) for high-level radioactive waste. The role played by colloids on RN migration is also of great concern. Diffusion and sorption are fundamental mechanisms controlling RN migration in rocks and many experimental approaches are applied to determine transport parameters for low sorbing RN in homogeneous rocks. However, it is difficult to obtain relevant data for high sorbing RN or colloids, for which diffusion lengths are extremely short, or within heterogeneous rocks, where transport might be different in different minerals. The ion beam techniques Rutherford Backscattering Spectrometry (RBS) and micro-Particle Induced X-Ray Emission ({mu}PIXE), rarely applied in the field, were selected for their micro-analytical potential to study RN diffusion and surface retention within heterogeneous rocks. Main achievements obtained during last 12 years are highlighted.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Challenges and opportunities for fractured rock imaging using 3D cross-borehole electrical resistivity

There is an increasing need to characterize discrete fractures away from boreholes to better define fracture distributions and monitor solute transport. We performed a 3D evaluation of static and time-lapse cross-borehole electrical resistivity tomography (ERT) data sets from a limestone quarry in which flow and transport are controlled by a bedding-plane feature. Ten boreholes were discretized using an unstructured tetrahedral mesh, and 2D panel measurements were inverted for a 3D distribution of conductivity. We evaluated the benefits of 3D versus 2.5D inversion of ERT data in fractured rock while including the use of borehole regularization disconnects (BRDs) and borehole conductivity constraints. High-conductivity halos (inversion artifacts) surrounding boreholes were removed in static images when BRDs and borehole conductivity constraints were implemented. Furthermore, applying these constraints focused transient changes in conductivity resulting from solute transport on the bedding plane, providing a more physically reasonable model for conductivity changes associated with solute transport at this fractured rock site. Assuming bedding-plane continuity between fractures identified in borehole televiewer data, we discretized a planar region between six boreholes and applied a fracture regularization disconnect (FRD). Although the FRD appropriately focused conductivity changes on the bedding plane, the conductivity distribution within the discretized fracture was nonunique and dependent on the starting homogeneous model conductivity. Synthetic studies performed to better explain field observations showed that inaccurate electrode locations in boreholes resulted in low-conductivity halos surrounding borehole locations. These synthetic studies also showed that the recovery of the true conductivity within an FRD depended on the conductivity contrast between the host rock and fractures. Our findings revealed that the potential exists to improve imaging of fractured rock through 3D inversion and accurate modeling of boreholes. However, deregularization of localized features can result in significant electrical conductivity artifacts, especially when representing features with a high degree of spatial uncertainty.

resistivity, fractured rock, inversion, imaging, g↗
DOE OSTI2016

Application of Monte Carlo Methods to Perform Uncertainty and Sensitivity Analysis on Inverse Water-Rock Reactions with NETPATH

Methods were developed to quantify uncertainty and sensitivity for NETPATH inverse water-rock reaction models and to calculate dissolved inorganic carbon, carbon-14 groundwater travel times. The NETPATH models calculate upgradient groundwater mixing fractions that produce the downgradient target water chemistry along with amounts of mineral phases that are either precipitated or dissolved. Carbon-14 groundwater travel times are calculated based on the upgradient source-water fractions, carbonate mineral phase changes, and isotopic fractionation. Custom scripts and statistical code were developed for this study to facilitate modifying input parameters, running the NETPATH simulations, extracting relevant output, postprocessing the results, and producing graphs and summaries. The scripts read userspecified values for each constituent’s coefficient of variation, distribution, sensitivity parameter, maximum dissolution or precipitation amounts, and number of Monte Carlo simulations. Monte Carlo methods for analysis of parametric uncertainty assign a distribution to each uncertain variable, sample from those distributions, and evaluate the ensemble output. The uncertainty in input affected the variability of outputs, namely source-water mixing, phase dissolution and precipitation amounts, and carbon-14 travel time. Although NETPATH may provide models that satisfy the constraints, it is up to the geochemist to determine whether the results are geochemically reasonable. Two example water-rock reaction models from previous geochemical reports were considered in this study. Sensitivity analysis was also conducted to evaluate the change in output caused by a small change in input, one constituent at a time. Results were standardized to allow for sensitivity comparisons across all inputs, which results in a representative value for each scenario. The approach yielded insight into the uncertainty in water-rock reactions and travel times. For example, there was little variation in source-water fraction between the deterministic and Monte Carlo approaches, and therefore, little variation in travel times between approaches. Sensitivity analysis proved very useful for identifying the most important input constraints (dissolved-ion concentrations), which can reveal the variables that have the most influence on source-water fractions and carbon-14 travel times. Once these variables are determined, more focused effort can be applied to determining the proper distribution for each constraint. Second, Monte Carlo results for water-rock reaction modeling showed discrete and nonunique results. The NETPATH models provide the solutions that satisfy the constraints of upgradient and downgradient water chemistry. There can exist multiple, discrete solutions for any scenario and these discrete solutions cause grouping of results. As a result, the variability in output may not easily be represented by a single distribution or a mean and variance and care should be taken in the interpretation and reporting of results.

phase dissolution and precipitation amounts↗