Strength, Fracture Evolution, and Permeability Changes from Confined Brazilian Tests on Sandstone
Understand failure conditions of rock in the confined extension region
SEARCH · Engineering Papers
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.
Understand failure conditions of rock in the confined extension region
The incompleteness of the fossil record hinders the inference of evolutionary rates and patterns. Here, we derive relationships among true taxonomic durations, preservation probability, and observed taxonomic ranges. We use these relationships to estimate original distributions of taxonomic durations, preservation probability, and completeness (proportion of taxa preserved), given only the observed ranges. No data on occurrences within the ranges of taxa are required. When preservation is random and the original distribution of durations is exponential, the inference of durations, preservability, and completeness is exact. However, reasonable approximations are possible given non-exponential duration distributions and temporal and taxonomic variation in preservability. Thus, the approaches we describe have great potential in studies of taphonomy, evolutionary rates and patterns, and genealogy. Analyses of Upper Cambrian-Lower Ordovician trilobite species, Paleozoic crinoid genera, Jurassic bivalve species, and Cenozoic mammal species yield the following results: (1) The preservation probability inferred from stratigraphic ranges alone agrees with that inferred from the analysis of stratigraphic gaps when data on the latter are available. (2) Whereas median durations based on simple tabulations of observed ranges are biased by stratigraphic resolution, our estimates of median duration, extinction rate, and completeness are not biased.(3) The shorter geologic ranges of mammalian species relative to those of bivalves cannot be attributed to a difference in preservation potential. However, we cannot rule out the contribution of taxonomic practice to this difference. (4) In the groups studied, completeness (proportion of species [trilobites, bivalves, mammals] or genera [crinoids] preserved) ranges from 60% to 90%. The higher estimates of completeness at smaller geographic scales support previous suggestions that the incompleteness of the fossil record reflects loss of fossiliferous rock more than failure of species to enter the fossil record in the first place.
The Mars2020 rover includes a Sampling & Caching Subsystem (SCS) that is required to acquire core samples of hard rocks and to abrade rocky surfaces flat and smooth in order to facilitate surface science. Hard, high strength rocks exhibit brittle failure modes, and can be drilled and chiseled efficiently using percussive impacts. The coring drill percussion mechanism was designed to provide those impacts with variable output force and to be used with various drill bits. A unique hammering mechanism, using a base-driven hammer, powered by a conventional rotary actuator driving a Scotch yoke, was designed to robustly provide the required function in a Martian environment for a life exceeding 2 million impacts.
Abstract Micromechanical modeling of geomaterials is challenging because of the complex geometry of discontinuities and potentially large number of deformable material bodies that contact each other dynamically. In this study, we have developed a numerical approach for micromechanical analysis of deformable geomaterials with dynamic contacts. In our approach, we detect contacts among multiple blocks with arbitrary shapes, enforce different contact constraints for three different contact states of separated, bonded, and sliding, and iterate within each time step to ensure convergence of contact states. With these features, we are able to simulate the dynamic contact evolution at the microscale for realistic geomaterials having arbitrary shapes of grains and interfaces. We demonstrate the capability with several examples, including a rough fracture with different geometric surface asperity characteristics, settling of clay aggregates, compaction of a loosely packed sand, and failure of an intact marble sample. With our model, we are able to accurately analyze (1) large displacements and/or deformation, (2) the process of high stress accumulated at contact areas, (3) the failure of a mineral cemented rock samples under high stress, and (4) post-failure fragmentation. The analysis highlights the importance of accurately capturing (1) the sequential evolution of geomaterials responding to stress as motion, deformation, and high stress; (2) large geometric features outside the norms (such as large asperities and sharp corners) as such features can dominate the micromechanical behavior; and (3) different mechanical behavior between loosely packed and tightly packed granular systems.
Rocks under stress deform by creep mechanisms that include formation and slip on small-scale internal cracks. Intragranular cracks and slip along grain contacts release energy as elastic waves termed acoustic emissions (AE). AEs are thought to contain predictive information that can be used for fault failure forecasting. Here, we present a method using unsupervised classification and an attention network to forecast labquakes using AE waveform features. Our data were generated in a laboratory setting using a biaxial shearing device with granular fault gouge intended to mimic the conditions of tectonic faults. Here, we analyzed the temporal evolution of AEs generated throughout several hundred laboratory earthquake cycles. We used a Conscience Self-Organizing Map (CSOM) to perform topologically ordered vector quantization based on waveform properties. The resulting map was used to interactively cluster AEs. We examined the clusters over time to identify those with predictive ability. Finally, we used a variety of LSTM and attention-based networks to test the predictive power of the AE clusters. By tracking cumulative waveform features over the seismic cycle, the network is able to forecast the time-to-failure (TTF) of lab earthquakes. Our results show that analyzing the data to isolate predictive signals and using a more sophisticated network architecture are key to robustly forecasting labquakes. In the future, this method could be applied on tectonic faults to monitor earthquakes and augment early warning systems.
Crack initiation, growth, and coalescence in flawed rocks have been extensively studied for 2D (planar, penetrating) flaws under uniaxial/biaxial compression. However, little is known as to the mechanisms and processes of cracking from 3D flaws under true triaxial compression, where the intermediate principal stress (σ 2 ) is distinguished from the major and minor principal stresses. In this work, we systematically investigate the effects of σ 2 on the 3D cracking behavior of rock specimens with preexisting flaws, through the use of mechanistic simulations of mixed-mode fracture in rocks. We explore how two characteristics of σ 2 , namely, (i) its orientation with respect to the flaw and (ii) its magnitude, affect two aspects of the cracking behavior, namely, (i) the cracking pattern and (ii) the peak stress. Results show that the orientation of σ 2 exerts more control over the cracking pattern than the flaw inclination angle. The peak stress becomes highest when σ 1 is parallel to the flaw, whereas it becomes lowest when σ 2 is parallel to the flaw. Also, the effects of σ 2 magnitude are more significant when σ 2 becomes more oblique to the flaw plane. On the basis of our observations, we propose mechanisms underlying the cracking behavior of 3D flawed rocks under true triaxial compression.
The Source Physics Experiment (SPE) is a long-term NNSA research and development effort designed to improve nonproliferation verification and monitoring capabilities. The overarching goals of the SPE program are to improve understanding of prompt signals and physical signatures that develop from underground chemical explosions and associated modeling capabilities. Our work focuses on a primary factor controlling chemical explosion induced signals and signatures: the material properties of the rocks in which the chemical explosion takes place. This document reports on material property determinations of legacy core USGS Test Well F and outcrop analogs for the subsurface stratigraphy for the third phase of SPE in the Rock Valley (RV) area of the Nevada National Security Site (NNSS). The objective of this work is to establish a baseline set of lithologic descriptions and material properties expected prior to observatory borehole drilling in support of the SPE-RVDC (Rock Valley Direct Comparison) experiment. We determine for each rock type the compressional failure envelope, elastic properties as a function of stress (bulk modulus versus mean stress, shear modulus versus shear stress, Young’s modulus versus axial stress and Poisson’s ratio versus axial stress), indirect tensile strength, and porosity. Geologic characterization, both at the core-scale and microscale, provides context for using the data in modeling efforts and to inform interpretations for the material properties testing.
In this paper, a scale-invariant bond model for simulating two-dimensional large deformation of continua is presented. An accurate torque form caused by relative rotation between bonded particles is derived based on beam theory. Validations of the proposed bond model are carried out by large deformation simulations of a simply supported beam subjected to central loading and a cantilever beam with end loading. Comparison of simulations with analytical solutions show a very good agreement. Later, Brazilian tests are simulated to demonstrate the ability of the framework to model failure and fracture process of rock.
The objective of this work is to utilize integrated geomechanics, field vertical seismic profile (VSP) and microseismic data to characterize the complex subsurface stress conditions at the Farnsworth Unit (FWU). The model is based on a five-spot sector model extracted from a primary geomechanical model. The five-spot well injection pattern is characterized by extensive reservoir characterization data, such well logs, extracted cores and borehole geophone data, to facilitate the detailed examination of stress changes and microseismic event occurrences. The study utilizes field vertical seismic volumes acquired from the injection well 13-10A. The seismic volumes successfully provided snapshots of the behavior of the reservoir at distinct times. The use of VSP and microseismic data provided direct and indirect estimates of the dynamic stress changes occurring in the overburden, reservoir and underburden rock formations. In order to illuminate the stress regions and identify rocks that have undergone inelastic failure, microseismic event occurrences were utilized. Microseismic activity has been detected at the FWU; further study of its locations, timing, and magnitude was needed to deduce the nature of the changing stress state. The results of the study revealed that microseismic events were successfully modeled within the Morrow B formation. Moment magnitudes of seismic events were within the same magnitudes for events in the reservoir, suggesting the suitability of the model. The results of the study showed that the computed moment magnitudes for seismic events were insignificant to warrant safety concerns. The study findings showed the usefulness of coupled hydromechanical models in predicting the subsurface stress changes associated with CO 2 injection. The knowledge gained from this study will serve as a guideline for industries planning to undertake underground CO 2 storage, and characterize the subsurface stress changes.
Research on the mechanics of lunar soils and on developing probes to determine the properties of lunar surface materials is summarized. The areas of investigation include the following: soil simulation, soil property determination using an impact penetrometer, soil stabilization using urethane foam or phenolic resin, effects of rolling boulders down lunar slopes, design of borehole jack and its use in determining failure mechanisms and properties of rocks, and development of a permeability probe for measuring fluid flow through porous lunar surface materials.
A number of geologic processes, particularly seismic faulting, impact crater slumping, and long runout landslides, require the failure of geologic materials under differential stresses much smaller than expected on the basis of conventional rock mechanics. This paper proposes that the low strengths apparent in these phenomena are due to a state of 'acoustic fluidization' induced by a transient strong acoustic wave field. The strain rates possible in such a field are evaluated, and it is shown that acoustically fluidized debris behaves as a newtonian fluid with a viscosity in the range 100,000 to 10,000,000 P for plausible conditions. Energy gains and losses in the acoustic field are discussed, and the mechanism is shown to be effective if internal dissipation in the field gives a Q approximately greater than 100. Whether such values for Q are realized is not known at present. However, acoustic fluidization provides a qualitatively correct description of the failure of rock debris under low differential stresses in the processes of faulting, crater slumping, and long runout landslides. Acoustic fluidization thus deserves serious consideration as a possible explanation of these phenomena.
Deformation localization is a widely observed, but rarely quantified process in the crust. Recent observations suggest that the localization of seismicity and fracture networks can help identify the approach to catastrophic failure. Here, we quantify the localization processes of the volumetric and deviatoric strain components in twelve triaxial compression experiments imaged with X-ray tomography. We capture three-dimensional images of the rock cores during triaxial compressing toward failure, and then calculate the local strain components using digital volume correlation. The divergence and curl of the incremental displacement vector field provide the volumetric and deviatoric components of the strain field. We quantify localization using the proportion of the rock occupied by high magnitudes of the volumetric and deviatoric strains, and the Gini coefficient of these high magnitude strains, which measures the deviation from a uniform process. We find that the vast majority, but not all, of the experiments experience strain localization toward failure. The rocks typically experience their maximum degree of strain localization not immediately preceding failure, but on average at 90% of the failure stress. The volumetric strain tends to localize more than the deviatoric strain. Finally, these observations support using the localization of the volumetric strain, along with the deviatoric strain, to identify the evolution of the precursory phase preceding earthquakes.
The principles of operation of holographic interferometry are reviewed, and some considerations affecting the design of high-pressure experiments using this technique are discussed. Some applications of holographic interferometry in high-pressure experiments are described, including strain measurements, investigation of failure phenomena, and investigation of piston rotation and sample bending in the early stages of a supposed uniaxial compression test. Specific attention is given to an experiment in which holographic interferometry was used to study the deformation of a jacketed rock sample while it was stressed to failure.
Surface properties, failure processes, and thermodynamic properties of rock in simulated lunar environments
Strength of rocks in the confined tension region, where the minimum principal stress is tensile, has only infrequently been measured and is not well understood. Quasi-static confined Brazilian tests under a range of confining stresses (2.76 to 27.58 MPa) where used to determine the strength of sandstone in the confined tension region. The test results indicate that the strength in the confined tension region was a strong function of the intermediate principal stress: increasing the intermediate principal stress significantly increased the strength of the sandstone. The strength data were well fit by the Mogi–Coulomb criterion, which accounts for the intermediate principal stress. Unconfined Brazilian strength data were not well fit to the Mogi–Coulomb criterion derived from the confined Brazilian test data, consistent with a transition from tensile to shear processes dominating failure with increasing confining pressure. Observations of post-failure fracture surfaces reveal more indication of shear processes with increasing confining pressure. Numerical simulations from combined finite-discrete element method are compared to the experimental results and reflect similar conditions for failure compared to the experimental tests in the confined tension region.
In the vincinity of earthquake focal regions, conditions may not be equal. Crustal rocks stressed to more than approximately 0.6 of their failure strength exhibit material properties over and above that of linear elasticity. Interactions between the earth tide and crustal rocks that are under high tectonic stress are discussed in terms of simple phenomenological models. In particular, the difference between a nonlinear elastic model of dilatancy and a dilatancy model that exhibits hysteresis is noted. It is concluded that the small changes in stress produced by the earth tide act as a probe of the properties of crustal rocks. Observations of earth tide tilts and strains in such high stress zones may, therefore, provide keys to the constitutive properties and the tectonic stress rate tensor of these zones.
Landslides are recognized as common triggers of glacial lake outburst floods in high mountain regions, but the remoteness of dangerous glacial lakes makes monitoring and disaster management very difficult. Using a very high resolution digital elevation model produced from a series of stereo satellite images, we simulate debris flows and rock avalanche runout around two very dangerous lakes in high mountain Nepal. We constrain particularly likely sources of these slope failures to simulate 79,800 potential debris flows and rock avalanches of volumes ranging from 1,000 m3 to 1,000,000 m3. This allowed us to identify zones in and around dangerous glacial lakes that are most exposed to landslide runout hazards, along with the most frequent source areas that produced landslides that may run into either of the two glacial lakes. Although our empirical-statistical model presented limitations that numerical models attempt to overcome, our probabilistic approach combined with a very high resolution DEM may shed new light on how to identify sources and impact areas of debris flows and rock avalanches that may cascade into downstream disasters.
Stress-strain tests were performed both on shock-damaged gabbro and limestone. The effective Young's modulus decreases with increasing initial damage parameter value, and an apparent work-softening process occurs prior to failure. To further characterize shock-induced microcracks, the longitudinal elastic wave velocity behavior of shock-damaged gabbro in the direction of compression up to failure was measured using an acoustic transmission technique under uniaxial loading. A dramatic increase in velocity was observed for the static compressive stress range of 0-50 MPa. Above that stress range, the velocity behavior of lightly damaged (D(sub 0) less than 0.1) gabbro is almost equal to unshocked gabbro. The failure strength of heavily-damaged (D(sub 0) greater than 0.1) gabbro is approx. 100-150 MPa, much lower than that of lightly damaged and unshocked gabbros (approx. 230-260 MPa). Following Nur's theory, the crack shape distribution was analyzed. The shock-induced cracks in gabbro appear to be largely thin penny-shaped cracks with c/a values below 5 x 10(exp -4). Moreover, the applicability of Ashby and Sammis's theory relating failure strength and damage parameter of shock-damaged rocks was examined and was found to yield a good estimate of the relation of shock-induced deficit in elastic modulus with the deficit in compressive strength.