Science.gov

Sample records for rock mechanics modelling

  1. Physical vs. Mathematical Models in Rock Mechanics

    NASA Astrophysics Data System (ADS)

    Morozov, I. B.; Deng, W.

    2013-12-01

    One of the less noted challenges in understanding the mechanical behavior of rocks at both in situ and lab conditions is the character of theoretical approaches being used. Currently, the emphasis is made on spatial averaging theories (homogenization and numerical models of microstructure), empirical models for temporal behavior (material memory, compliance functions and complex moduli), and mathematical transforms (Laplace and Fourier) used to infer the Q-factors and 'relaxation mechanisms'. In geophysical applications, we have to rely on such approaches for very broad spatial and temporal scales which are not available in experiments. However, the above models often make insufficient use of physics and utilize, for example, the simplified 'correspondence principle' instead of the laws of viscosity and friction. As a result, the commonly-used time- and frequency dependent (visco)elastic moduli represent apparent properties related to the measurement procedures and not necessarily to material properties. Predictions made from such models may therefore be inaccurate or incorrect when extrapolated beyond the lab scales. To overcome the above challenge, we need to utilize the methods of micro- and macroscopic mechanics and thermodynamics known in theoretical physics. This description is rigorous and accurate, uses only partial differential equations, and allows straightforward numerical implementations. One important observation from the physical approach is that the analysis should always be done for the specific geometry and parameters of the experiment. Here, we illustrate these methods on axial deformations of a cylindrical rock sample in the lab. A uniform, isotropic elastic rock with a thermoelastic effect is considered in four types of experiments: 1) axial extension with free transverse boundary, 2) pure axial extension with constrained transverse boundary, 3) pure bulk expansion, and 4) axial loading harmonically varying with time. In each of these cases, an

  2. A rock-/ice mechanical model for the destabilisation of permafrost rocks

    NASA Astrophysics Data System (ADS)

    Krautblatter, Michael; Funk, Daniel

    2010-05-01

    The destabilisation of permafrost rocks is commonly attributed to changes in ice-mechanical properties (Davies et al. 2001). The effect of low temperatures on intact rock strength and its mechanical relevance for shear strength and brittle fracture propagation has not been considered yet. But this effect is significant since compressive and tensile strength are reduced by up to 50% and more when rock thaws (Mellor, 1973). Here we show, that the reduction of the shear resistance of rock-rock contacts in joints plays a key role for the onset of larger instabilities in thawing permafrost rocks. Based on a Mohr-Coulomb assumption, we defined a failure criterion of an ice-filled rock cleft, with cohesive rock bridges, contact of rough fracture surfaces, ductile creep of ice and with a representation of rock-ice "failure" mechanisms along the surface and inside the ice body. The synoptic models are based on the principle of superposition, i.e. that shear stress "absorbed" by one component reduces the amount of shear stress applied to the other components. Failure along existing sliding planes can be explained by the impact of temperature on shear stress uptake by creep deformation of ice, the propensity of failure along rock-ice fractures and reduced total friction along rough rock-rock contacts. This model may account for the rapid response of rockslides to warming (reaction time). In the long term, brittle fracture propagation is initialised. Warming reduces the shear stress uptake by total friction and decreases the critical fracture toughness along rock bridges. The latter model accounts for slow subcritical destabilisation of whole rock slopes over decades to millennia, subsequent to the warming impulse (relaxation time). To test the importance of reduced friction, we conducted shearing tests on homogeneous fine-grained limestone specimen taken from a permafrost site (Zugspitze, Germany). In a temperature-controlled shearing box, we repeatedly tested mechanical

  3. Rock mechanics models evaluation report. [Contains glossary

    SciTech Connect

    Not Available

    1987-08-01

    This report documents the evaluation of the thermal and thermomechanical models and codes for repository subsurface design and for design constraint analysis. The evaluation was based on a survey of the thermal and thermomechanical codes and models that are applicable to subsurface design, followed by a Kepner-Tregoe (KT) structured decision analysis of the codes and models. The primary recommendations of the analysis are that the DOT code be used for two-dimensional thermal analysis and that the STEALTH and HEATING 5/6 codes be used for three-dimensional and complicated two-dimensional thermal analysis. STEALTH and SPECTROM 32 are recommended for thermomechanical analyses. The other evaluated codes should be considered for use in certain applications. A separate review of salt creep models indicate that the commonly used exponential time law model is appropriate for use in repository design studies. 38 refs., 1 fig., 7 tabs.

  4. Strain localisation in mechanically Layered Rocks, insights from numerical modelling

    NASA Astrophysics Data System (ADS)

    Le Pourhiet, L.; Huet, B.; Agard, P.; Labrousse, L.; Jolivet, L.; Yao, K.

    2012-09-01

    Small scale deformation in stratified rocks displays a large diversity of micro-structures, from the microscopic scale to the scale of orogens. We have designed a series of fully dynamic numerical simulations aimed at assessing which parameters control this structural diversity and which underlying mechanisms lead to strain localisation. The influence of stratification orientation on the occurrence and mode of strain localisation is tested by varying the initial dip of inherited layering versus the large scale imposed simple shear. The detailed study of the models indicates that (1) the results are length-scale independent, (2) the new shear zones are always compatible with the kinematics imposed at the boundary (3) micro-structures formed encompass the full diversity of micro-structures observed in the field and chiefly depend on the direction of the initial anisotropy versus shear direction, (4) depending on the orientation of the anisotropy, the layers may deform along subtractive or additive shear bands, (5) the deformation in anisotropic media results in non-lithostatic pressure values that are on the order of the deviatoric stress in the strong layers and (6) the introduction of brittle rheology is necessary to form localised shear bands in the ductile regime.

  5. Towards a mechanical failure model for degrading permafrost rock slopes representing changes in rock toughness and infill

    NASA Astrophysics Data System (ADS)

    Mamot, Philipp; Krautblatter, Michael; Scandroglio, Riccardo

    2016-04-01

    The climate-induced degradation of permafrost in mountain areas can reduce the stability of rock slopes. An increasing number of rockfalls and rockslides originate from permafrost-affected rock faces. Discontinuity patterns and their geometrical and mechanical properties play a decisive role in controlling rock slope stability. Under thawing conditions the shear resistance of rock reduces due to lower friction along rock-rock contacts, decreasing fracture toughness of rock-ice contacts, diminishing fracture toughness of cohesive rock bridges and altered creep or fracture of the ice itself. Compressive strength is reduced by 20 to 50 % and tensile strength decreases by 15 to 70 % when intact saturated rock thaws (KRAUTBLATTER ET AL. 2013). Elevated water pressures in fractures can lead to reduced effective normal stresses and thus to lower shear strengths of fractures. However, the impact of degrading permafrost on the mechanical properties of intact or fractured rock still remains poorly understood. In this study, we develop a new approach for modeling the influence of degrading permafrost on the stability of high mountain rock slopes. Hereby, we focus on the effect of rock- and ice-mechanical changes along striking discontinuities onto the whole rock slope. We aim at contributing to a better rock-ice mechanical process understanding of degrading permafrost rocks. For parametrisation and subsequent calibration of our model, we chose a test site (2885 m a.s.l.) close by the Zugspitze summit in Germany. It reveals i) a potential rockslide at the south face involving 10E4m³ of rock and ii) permafrost occurrence due to ice-filled caves and fractures. Here we combine kinematic, geotechnical and thermal monitoring in the field with rock-mechanical laboratory tests and a 2D numerical failure modeling. Up to date, the following results underline the potential effects of thawing rock and fracture infill on the stability of steep rock slopes in theory and praxis: i. ERT and

  6. Hydro-mechanically coupled modelling of deep-seated rock slides in the surroundings of reservoirs

    NASA Astrophysics Data System (ADS)

    Lechner, Heidrun; Preh, Alexander; Zangerl, Christian

    2016-04-01

    In order to enhance the understanding of the behaviour of deep-seated rock slides in the surroundings of large dam reservoirs, this study concentrates on failure mechanisms, deformation processes and the ability of self-stabilisation of rock slides influenced by reservoirs. Particular focus is put on internal rock mass deformations, progressive topographical slope changes due to reservoir impoundment and shear displacements along the basal shear zone in relation to its shear strength properties. In this study, a two-dimensional numerical rock slide model is designed by means of the Universal Distinct Element Code UDEC and investigated concerning different groundwater flow scenarios. These include: (i) a completely drained rock slide model, (ii) a model with fully saturated rock mass below an inclined groundwater table and (iii) a saturated groundwater model with a reservoir at the slope toe. Slope displacements initiate when the shear strength properties of the basal shear zone are at or below the critical parameters for the limit-equilibrium state and continue until a numerical equilibrium is reached due to deformation- and displacement-based geometrical changes. The study focuses on the influence of a reservoir at the toe of a rock slide and tries to evaluate the degree of displacement which is needed for a re-stabilisation in relation to the geometrical characteristics of the rock slide. Besides, challenges and limitations of applied distinct element methods to simulate large strain and displacements of deep-seated rock slides are discussed. The ongoing study will help to understand the deformation behaviour of deep-seated pre-existing rock slides in fractured rock mass during initial impounding and will be part of a hazard assessment for large reservoirs.

  7. A Real Two-Phase Mechanical Model for Rock-Ice Avalanches

    NASA Astrophysics Data System (ADS)

    Pudasaini, S. P.; Krautblatter, M.

    2012-04-01

    Rock-ice avalanches in high mountain permafrost environments are a hazardous and poorly understood process. Their hazard potential derives from the large volume, high velocities, the potential entrainment of large amounts of rock-debris, ice, snow and water during the flow, high impact pressures, and unpredictable flow paths and deposition patterns. In contrast to the usual single-phase model of rock avalanches, the solid phase (ice) in rock-ice avalanches can transform to fluid (water or slurry) during the course of the debris-avalanche and fundamentally alter the multiple mechanical processes. We postulate that a real two-phase debris flow model could much better address the dynamic interaction of solid (rock and ice) and fluid (water, snow, slurry and fine particles) rather than a simple single-phase Voellmy- or Coulomb-type model. For this, we enhance the general two-phase debris flow model proposed by Pudasaini (2011) by additionally introducing two new mechanical aspects typical for the rock-ice avalanches: (a) the dynamic strength weakening including the internal fluidization and basal lubrication, as well as (b) the internal mass and momentum exchanges between the phases. In these models, the effective basal and internal friction angles are variable and are described in terms of evolving effective solid volume fraction (rock and ice), friction factors, volume fraction of the ice, true friction coefficients and the lubrication and fluidization factors. These factors are functions of several physical parameters and mechanical and dynamical variables, including the volume fractions of the solid, shear-rate and the normal stresses. Rock-ice avalanches are a unique scenario in geophysical mass flows, where phase exchange and material strength weakening occurs and can dominate the flow dynamics. Here, we present an innovative approach to model and simulate these two special aspects. Additionally, in the model, the inertial terms include the hydraulic pressure

  8. The three-dimension model for the rock-breaking mechanism of disc cutter and analysis of rock-breaking forces

    NASA Astrophysics Data System (ADS)

    Zhang, Zhao-Huang; Sun, Fei

    2012-06-01

    To study the rock deformation with three-dimensional model under rolling forces of disc cutter, by carrying out the circular-grooving test with disc cutter rolling around on the rock, the rock mechanical behavior under rolling disc cutter is studied, the mechanical model of disc cutter rolling around the groove is established, and the theory of single-point and double-angle variables is proposed. Based on this theory, the physics equations and geometric equations of rock mechanical behavior under disc cutters of tunnel boring machine (TBM) are studied, and then the balance equations of interactive forces between disc cutter and rock are established. Accordingly, formulas about normal force, rolling force and side force of a disc cutter are derived, and their validity is studied by tests. Therefore, a new method and theory is proposed to study rock-breaking mechanism of disc cutters.

  9. Mechanical defradation of Emplacement Drifts at Yucca Mountain- A Modeling Case Study. Part I: Nonlithophysal Rock

    SciTech Connect

    M. Lin; D. Kicker; B. Damjanac; M. Board; M. Karakouzian

    2006-07-05

    This paper outlines rock mechanics investigations associated with mechanical degradation of planned emplacement drifts at Yucca Mountain, which is the designated site for the proposed U.S. high-level nuclear waste repository. The factors leading to drift degradation include stresses from the overburden, stresses induced by the heat released from the emplaced waste, stresses due to seismically related ground motions, and time-dependent strength degradation. The welded tuff emplacement horizon consists of two groups of rock with distinct engineering properties: nonlithophysal units and lithophysal units, based on the relative proportion of lithophysal cavities. The term 'lithophysal' refers to hollow, bubble like cavities in volcanic rock that are surrounded by a porous rim formed by fine-grained alkali feldspar, quartz, and other minerals. Lithophysae are typically a few centimeters to a few decimeters in diameter. Part I of the paper concentrates on the generally hard, strong, and fractured nonlithophysal rock. The degradation behavior of the tunnels in the nonlithophysal rock is controlled by the occurrence of keyblocks. A statistically equivalent fracture model was generated based on extensive underground fracture mapping data from the Exploratory Studies Facility at Yucca Mountain. Three-dimensional distinct block analyses, generated with the fracture patterns randomly selected from the fracture model, were developed with the consideration of in situ, thermal, and seismic loads. In this study, field data, laboratory data, and numerical analyses are well integrated to provide a solution for the unique problem of modeling drift degradation.

  10. Effect of Particle Shape on Mechanical Behaviors of Rocks: A Numerical Study Using Clumped Particle Model

    PubMed Central

    Rong, Guan; Liu, Guang; Zhou, Chuang-bing

    2013-01-01

    Since rocks are aggregates of mineral particles, the effect of mineral microstructure on macroscopic mechanical behaviors of rocks is inneglectable. Rock samples of four different particle shapes are established in this study based on clumped particle model, and a sphericity index is used to quantify particle shape. Model parameters for simulation in PFC are obtained by triaxial compression test of quartz sandstone, and simulation of triaxial compression test is then conducted on four rock samples with different particle shapes. It is seen from the results that stress thresholds of rock samples such as crack initiation stress, crack damage stress, and peak stress decrease with the increasing of the sphericity index. The increase of sphericity leads to a drop of elastic modulus and a rise in Poisson ratio, while the decreasing sphericity usually results in the increase of cohesion and internal friction angle. Based on volume change of rock samples during simulation of triaxial compression test, variation of dilation angle with plastic strain is also studied. PMID:23997677

  11. Modeling the Progressive Failure of Jointed Rock Slope Using Fracture Mechanics and the Strength Reduction Method

    NASA Astrophysics Data System (ADS)

    Zhang, Ke; Cao, Ping; Meng, Jingjing; Li, Kaihui; Fan, Wenchen

    2015-03-01

    The fracturing process during the progressive failure of a jointed rock slope is numerically investigated by using fracture mechanics and the strength reduction method (SRM). A displacement discontinuity method containing frictional elements is developed for the calculation of the stress intensity factor (SIF). The failure initiation of the jointed rock slope is analyzed by evaluating the SIF. A new joint model is proposed by combining solid elements with interface elements in the commercial software FLAC3D. These represent the discontinuous planes in a rock mass on which sliding or separation can occur. The progressive failure process is simulated by reducing the shear strength of the rock mass, which includes the process of stress concentration, crack initiation, crack propagation, slip weakening, and coalescence of failure surfaces. The factor of safety (FS) and location of the critical failure surface are determined by the SRM. The influence of the joint inclination is investigated using the FS and the SIF. Laboratory experiments on specimens containing an inclined flaw under compression-shear stress are also conducted to investigate the effect of the angle between the shear direction and the flaw inclination, which provides an experimental explanation for the shear behavior of jointed rock. The results show that the joint inclination dominates the failure behavior of jointed rock slope, and two failure patterns have been classified.

  12. Experimental investigations for the modeling of chemo-mechanical processes of anhydritic rock

    NASA Astrophysics Data System (ADS)

    Huber, Tara; Pimentel, Erich; Anagnostou, Georg

    2015-04-01

    When anhydritic rock comes into contact with water, the anhydrite dissolves and gypsum precipitates as a result of an oversaturation of the sulphate and calcium ions in the water. This anhydrite to gypsum transformation (AGT) leads to an increase in the solid volume by roughly 61% and possibly also of the pore volume, thus resulting in macroscopic swelling. In tunneling, swelling rock can cause massive damage, since it can exhibit high pressures on the lining or result in large deformations of the lining. Even though this phenomenon has already been observed more than a century ago, AGT in sulphatic rock still raises many open questions. One question we focus on is the mechanical (i.e. stress, strain) behavior of anhydritic rock during swelling. The coupled chemo-mechanical processes, i.e. the development of stresses and strains during AGT, has not been fully understood so far. A chemo-mechanical model must be considered which take anhydrite dissolution, gypsum precipitation as well as stresses and strains into account. In this contribution we present some fundamental research which was done in order to establish a relationship between AGT and the resulting or applied strains and stresses. The research contains experiments on samples consisting of anhydrite and kaolin under oedometric conditions. In order to reduce uncertainties due to swelling of clay as well as inhomogeneous compositions and structures of the natural rock samples, the experiments at the present stage of this research are performed on artificially created, reproducible samples. The samples contain 40% industrial anhydritic powder and 60% Polwhite E Chinaclay (of which the main component is kaolinite). The powders are mixed and compacted in a steel ring under high axial pressure, thus creating intact discs with a dry density of roughly 1.9 g/cm3. In a first series of Oedometer tests the swelling strain under various constant axial stress is measured until the maximal strain is reached. At the

  13. E. coli RS2GFP Retention Mechanisms in Laboratory-Scale Fractured Rocks: A Statistical Model

    NASA Astrophysics Data System (ADS)

    Rodrigues, S. N.; Qu, J.; Dickson, S. E.

    2011-12-01

    during the visualization experiments will help to draw conclusions regarding the retention mechanisms at play in the fractures. This paper will present a statistical model for the retention of E. coli RS2GFP in fractures as a function of aperture field characteristics and flow rate. The images of the E. coli RS2GFP travelling through the casts, in combination with direct aperture field measurements, will provide insight into the different mechanisms instrumental to particulate retention within the fractures. The results of these experiments will partially bridge the knowledge gap in fractured aquifers by furthering the understanding of transport and retention mechanisms of microorganisms within fractured rocks.

  14. A Dynamic Damage Mechanics Source Model for Explosions in Crystalline Rock

    NASA Astrophysics Data System (ADS)

    Mihaly, J. M.; Bhat, H. S.; Sammis, C. G.; Rosakis, A.

    2011-12-01

    The micromechanical damage mechanics formulated by Ashby and Sammis [PAGEOPH, 1990] and generalized by Deshpande and Evans [J. Mech. Phys. Solids, 2008] has been extended to allow for a more generalized stress state and to incorporate an experimentally motivated crack growth (damage evolution) law that is valid over a wide range of loading rates. This law is sensitive to both the crack tip stress field and its time derivative, and thus produces strain-rate sensitivity in the constitutive response. The model is experimentally verified by predicting the failure strength of Dionysus-Pentelicon marble over strain rates ranging from to . This rate-dependent damage mechanics has been implemented in the ABAQUS dynamic finite element code and used to explore the effects of burn rate (loading rate) and lithostatic stress on the spatial extent of fracture damage and S waves generated by explosions in crystalline rock. Slower rise times and longer pressure pulses produce more damage and stronger S waves.

  15. A numerical model of hydro-thermo-mechanical coupling in a fractured rock mass

    SciTech Connect

    Bower, K.M.

    1996-06-01

    Coupled hydro-thermo-mechanical codes with the ability to model fractured materials are used for predicting groundwater flow behavior in fractured aquifers containing thermal sources. The potential applications of such a code include the analysis of groundwater behavior within a geothermal reservoir. The capability of modeling hydro-thermo systems with a dual porosity, fracture flow model has been previously developed in the finite element code, FEHM. FEHM has been modified to include stress coupling with the dual porosity feature. FEHM has been further developed to implicitly couple the dependence of fracture hydraulic conductivity on effective stress within two dimensional, saturated aquifers containing fracture systems. The cubic law for flow between parallel plates was used to model fracture permeability. The Bartin-Bandis relationship was used to determine the fracture aperture within the cubic law. The code used a Newton Raphson iteration to implicitly solve for six unknowns at each node. Results from a model of heat flow from a reservoir to the moving fluid in a single fracture compared well with analytic results. Results of a model showing the increase in fracture flow due to a single fracture opening under fluid pressure compared well with analytic results. A hot dry rock, geothermal reservoir was modeled with realistic time steps indicating that the modified FEHM code does successfully model coupled flow problems with no convergence problems.

  16. Observations, models, and mechanisms of failure of surface rocks surrounding planetary surface loads

    NASA Technical Reports Server (NTRS)

    Schultz, R. A.; Zuber, M. T.

    1994-01-01

    Geophysical models of flexural stresses in an elastic lithosphere due to an axisymmetric surface load typically predict a transition with increased distance from the center of the load of radial thrust faults to strike-slip faults to concentric normal faults. These model predictions are in conflict with the absence of annular zones of strike-slip faults around prominent loads such as lunar maria, Martian volcanoes, and the Martian Tharsis rise. We suggest that this paradox arises from difficulties in relating failure criteria for brittle rocks to the stress models. Indications that model stresses are inappropriate for use in fault-type prediction include (1) tensile principal stresses larger than realistic values of rock tensile strength, and/or (2) stress differences significantly larger than those allowed by rock-strength criteria. Predictions of surface faulting that are consistent with observations can be obtained instead by using tensile and shear failure criteria, along with calculated stress differences and trajectories, with model stress states not greatly in excess of the maximum allowed by rock fracture criteria.

  17. THERMO-HYDRO-MECHANICAL MODELING OF WORKING FLUID INJECTION AND THERMAL ENERGY EXTRACTION IN EGS FRACTURES AND ROCK MATRIX

    SciTech Connect

    Robert Podgorney; Chuan Lu; Hai Huang

    2012-01-01

    Development of enhanced geothermal systems (EGS) will require creation of a reservoir of sufficient volume to enable commercial-scale heat transfer from the reservoir rocks to the working fluid. A key assumption associated with reservoir creation/stimulation is that sufficient rock volumes can be hydraulically fractured via both tensile and shear failure, and more importantly by reactivation of naturally existing fractures (by shearing), to create the reservoir. The advancement of EGS greatly depends on our understanding of the dynamics of the intimately coupled rock-fracture-fluid-heat system and our ability to reliably predict how reservoirs behave under stimulation and production. Reliable performance predictions of EGS reservoirs require accurate and robust modeling for strongly coupled thermal-hydrological-mechanical (THM) processes. Conventionally, these types of problems have been solved using operator-splitting methods, usually by coupling a subsurface flow and heat transport simulators with a solid mechanics simulator via input files. An alternative approach is to solve the system of nonlinear partial differential equations that govern multiphase fluid flow, heat transport, and rock mechanics simultaneously, using a fully coupled, fully implicit solution procedure, in which all solution variables (pressure, enthalpy, and rock displacement fields) are solved simultaneously. This paper describes numerical simulations used to investigate the poro- and thermal- elastic effects of working fluid injection and thermal energy extraction on the properties of the fractures and rock matrix of a hypothetical EGS reservoir, using a novel simulation software FALCON (Podgorney et al., 2011), a finite element based simulator solving fully coupled multiphase fluid flow, heat transport, rock deformation, and fracturing using a global implicit approach. Investigations are also conducted on how these poro- and thermal-elastic effects are related to fracture permeability

  18. Strain localisation in mechanically layered rocks beneath detachment zones: insights from numerical modelling

    NASA Astrophysics Data System (ADS)

    Le Pourhiet, L.; Huet, B.; Labrousse, L.; Yao, K.; Agard, P.; Jolivet, L.

    2013-04-01

    We have designed a series of fully dynamic numerical simulations aimed at assessing how the orientation of mechanical layering in rocks controls the orientation of shear bands and the depth of penetration of strain in the footwall of detachment zones. Two parametric studies are presented. In the first one, the influence of stratification orientation on the occurrence and mode of strain localisation is tested by varying initial dip of inherited layering in the footwall with regard to the orientation of simple shear applied at the rigid boundary simulating a rigid hanging wall, all scaling and rheological parameter kept constant. It appears that when Mohr-Coulomb plasticity is being used, shear bands are found to localise only when the layering is being stretched. This corresponds to early deformational stages for inital layering dipping in the same direction as the shear is applied, and to later stages for intial layering dipping towards the opposite direction of shear. In all the cases, localisation of the strain after only γ=1 requires plastic yielding to be activated in the strong layer. The second parametric study shows that results are length-scale independent and that orientation of shear bands is not sensitive to the viscosity contrast or the strain rate. However, decreasing or increasing strain rate is shown to reduce the capacity of the shear zone to localise strain. In the later case, the strain pattern resembles a mylonitic band but the rheology is shown to be effectively linear. Based on the results, a conceptual model for strain localisation under detachment faults is presented. In the early stages, strain localisation occurs at slow rates by viscous shear instabilities but as the layered media is exhumed, the temperature drops and the strong layers start yielding plastically, forming shear bands and localising strain at the top of the shear zone. Once strain localisation has occured, the deformation in the shear band becomes extremely penetrative but

  19. Rock Cavern Stability Analysis Under Different Hydro-Geological Conditions Using the Coupled Hydro-Mechanical Model

    NASA Astrophysics Data System (ADS)

    Chen, H. M.; Zhao, Z. Y.; Choo, L. Q.; Sun, J. P.

    2016-02-01

    Rock cavern stability has a close relationship with the uncertain geological parameters, such as the in situ stress, the joint configurations, and the joint mechanical properties. Therefore, the stability of the rock cavern should be studied with variable geological conditions. In this paper, the coupled hydro-mechanical model, which is under the framework of the discontinuous deformation analysis, is developed to study the underground cavern stability when considering the hydraulic pressure after excavation. Variable geological conditions are taken into account to study their impacts on the seepage rate and the cavern stability, including the in situ stress ratio, joint spacing, and joint dip angle. In addition, the two cases with static hydraulic pressure and without hydraulic pressure are also considered for the comparison. The numerical simulations demonstrate that the coupled approach can capture the cavern behavior better than the other two approaches without the coupling effects.

  20. U.S. National Committee for Rock Mechanics; and Conceptual model of fluid infiltration in fractured media. Project summary, July 28, 1997--July 27, 1998

    SciTech Connect

    1998-09-01

    The title describes the two tasks summarized in this report. The remainder of the report contains information on meetings held or to be held on the subjects. The US National Committee for Rock Mechanics (USNC/RM) provides for US participation in international activities in rock mechanics, principally through adherence to the International Society for Rock Mechanics (ISRM). It also keeps the US rock mechanics community informed about new programs directed toward major areas of national concern in which rock mechanics problems represent critical or limiting factors, such as energy resources, excavation, underground storage and waste disposal, and reactor siting. The committee also guides or produces advisory studies and reports on problem areas in rock mechanics. A new panel under the auspices of the US National Committee for Rock Mechanics has been appointed to conduct a study on Conceptual Models of Fluid Infiltration in Fractured Media. The study has health and environmental applications related to the underground flow of pollutants through fractured rock in and around mines and waste repositories. Support of the study has been received from the US Nuclear Regulatory Commission and the Department of Energy`s Yucca Mountain Project Office. The new study builds on the success of a recent USNC/RM report entitled Rock Fractures and Fluid Flow: Contemporary Understanding and Applications (National Academy Press, 1996, 551 pp.). A summary of the new study is provided.

  1. Reaction-transport-mechanical (RTM) simulator Sym.CS: Putting together water-rock interaction, multi-phase and heat flow, composite petrophysics model, and fracture mechanics

    NASA Astrophysics Data System (ADS)

    Paolini, C.; Park, A. J.; Mellors, R. J.; Castillo, J.

    2009-12-01

    A typical CO2 sequestration scenario involves the use of multiple simulators for addressing multiphase fluid and heat flow, water-rock interaction and mass-transfer, rock mechanics, and other chemical and physical processes. The benefit of such workflow is that each model can be constrained rigorously; however, the drawback is final modeling results may achieve only a limited extent of the theoretically possible capabilities of each model. Furthermore, such an approach in modeling carbon sequestration cannot capture the nonlinearity of the various chemical and physical processes. Hence, the models can only provide guidelines for carbon sequestration processes with large margins of error. As an alternative, a simulator is being constructed by a multi-disciplinary team with the aim of implementing a large array of fundamental phenomenologies, including, but not limited to: water-rock interaction using elemental mass-balance and explicit mass-transfer and reaction coupling methods; multi-phase and heat flow, including super-critical CO2 and oil; fracture mechanics with anisotropic permeabilities; rheological rock mechanics based on incremental stress theory; and a composite petrophysics model capable of describing changing rock composition and properties. The modules representing the processes will be solved using a layered iteration method, with the goal of capturing the nonlinear feedback among all of the processes. The simulator will be constructed using proven optimization and modular, object-oriented, and service-oriented programming methods. Finally, a novel AJAX (asynchronous JavaScript and XML) user interface is being tested to host the simulator that will allow usage through an Internet browser. Currently, the water-rock interaction, composite petrophysics, and multi-phase fluid and heat flow modules are available for integration. Results of the water-rock interaction and petrophysics coupling has been used to model interaction between a CO2-charged water and

  2. Mechanism of Rock Burst Occurrence in Specially Thick Coal Seam with Rock Parting

    NASA Astrophysics Data System (ADS)

    Wang, Jian-chao; Jiang, Fu-xing; Meng, Xiang-jun; Wang, Xu-you; Zhu, Si-tao; Feng, Yu

    2016-05-01

    Specially thick coal seam with complex construction, such as rock parting and alternative soft and hard coal, is called specially thick coal seam with rock parting (STCSRP), which easily leads to rock burst during mining. Based on the stress distribution of rock parting zone, this study investigated the mechanism, engineering discriminant conditions, prevention methods, and risk evaluation method of rock burst occurrence in STCSRP through setting up a mechanical model. The main conclusions of this study are as follows. (1) When the mining face moves closer to the rock parting zone, the original non-uniform stress of the rock parting zone and the advancing stress of the mining face are combined to intensify gradually the shearing action of coal near the mining face. When the shearing action reaches a certain degree, rock burst easily occurs near the mining face. (2) Rock burst occurrence in STCSRP is positively associated with mining depth, advancing stress concentration factor of the mining face, thickness of rock parting, bursting liability of coal, thickness ratio of rock parting to coal seam, and difference of elastic modulus between rock parting and coal, whereas negatively associated with shear strength. (3) Technologies of large-diameter drilling, coal seam water injection, and deep hole blasting can reduce advancing stress concentration factor, thickness of rock parting, and difference of elastic modulus between rock parting and coal to lower the risk of rock burst in STCSRP. (4) The research result was applied to evaluate and control the risk of rock burst occurrence in STCSRP.

  3. Rock mechanics for hard rock nuclear waste repositories

    SciTech Connect

    Heuze, F.E.

    1981-09-01

    The mined geologic burial of high level nuclear waste is now the favored option for disposal. The US National Waste Terminal Storage Program designed to achieve this disposal includes an extensive rock mechanics component related to the design of the wastes repositories. The plan currently considers five candidate rock types. This paper deals with the three hard rocks among them: basalt, granite, and tuff. Their behavior is governed by geological discontinuities. Salt and shale, which exhibit behavior closer to that of a continuum, are not considered here. This paper discusses both the generic rock mechanics R and D, which are required for repository design, as well as examples of projects related to hard rock waste storage. The examples include programs in basalt (Hanford/Washington), in granitic rocks (Climax/Nevada Test Site, Idaho Springs/Colorado, Pinawa/Canada, Oracle/Arizona, and Stripa/Sweden), and in tuff (Nevada Test Site).

  4. Rock Properties Model

    SciTech Connect

    C. Lum

    2004-09-16

    The purpose of this model report is to document the Rock Properties Model version 3.1 with regard to input data, model methods, assumptions, uncertainties and limitations of model results, and qualification status of the model. The report also documents the differences between the current and previous versions and validation of the model. The rock properties model provides mean matrix and lithophysae porosity, and the cross-correlated mean bulk density as direct input to the ''Saturated Zone Flow and Transport Model Abstraction'', MDL-NBS-HS-000021, REV 02 (BSC 2004 [DIRS 170042]). The constraints, caveats, and limitations associated with this model are discussed in Section 6.6 and 8.2. Model validation accomplished by corroboration with data not cited as direct input is discussed in Section 7. The revision of this model report was performed as part of activities being conducted under the ''Technical Work Plan for: The Integrated Site Model, Revision 05'' (BSC 2004 [DIRS 169635]). The purpose of this revision is to bring the report up to current procedural requirements and address the Regulatory Integration Team evaluation comments. The work plan describes the scope, objectives, tasks, methodology, and procedures for this process.

  5. Training and Research on Probabilistic Hydro-Thermo-Mechanical Modeling of Carbon Dioxide Geological Sequestration in Fractured Porous Rocks

    SciTech Connect

    Gutierrez, Marte

    2013-05-31

    Colorado School of Mines conducted research and training in the development and validation of an advanced CO{sub 2} GS (Geological Sequestration) probabilistic simulation and risk assessment model. CO{sub 2} GS simulation and risk assessment is used to develop advanced numerical simulation models of the subsurface to forecast CO2 behavior and transport; optimize site operational practices; ensure site safety; and refine site monitoring, verification, and accounting efforts. As simulation models are refined with new data, the uncertainty surrounding the identified risks decrease, thereby providing more accurate risk assessment. The models considered the full coupling of multiple physical processes (geomechanical and fluid flow) and describe the effects of stochastic hydro-mechanical (H-M) parameters on the modeling of CO{sub 2} flow and transport in fractured porous rocks. Graduate students were involved in the development and validation of the model that can be used to predict the fate, movement, and storage of CO{sub 2} in subsurface formations, and to evaluate the risk of potential leakage to the atmosphere and underground aquifers. The main major contributions from the project include the development of: 1) an improved procedure to rigorously couple the simulations of hydro-thermomechanical (H-M) processes involved in CO{sub 2} GS; 2) models for the hydro-mechanical behavior of fractured porous rocks with random fracture patterns; and 3) probabilistic methods to account for the effects of stochastic fluid flow and geomechanical properties on flow, transport, storage and leakage associated with CO{sub 2} GS. The research project provided the means to educate and train graduate students in the science and technology of CO{sub 2} GS, with a focus on geologic storage. Specifically, the training included the investigation of an advanced CO{sub 2} GS simulation and risk assessment model that can be used to predict the fate, movement, and storage of CO{sub 2} in

  6. Fundamental Study on Applicability of Powder-Based 3D Printer for Physical Modeling in Rock Mechanics

    NASA Astrophysics Data System (ADS)

    Fereshtenejad, Sayedalireza; Song, Jae-Joon

    2016-06-01

    Applications of 3D printing technology become more widespread in many research fields because of its rapid development and valuable capabilities. In rock mechanics and mining engineering, this technology has the potential to become a useful tool that might help implement a number of research studies previously considered impractical. Most commercial 3D printers cannot print prototypes with mechanical properties that match precisely those of natural rock samples. Therefore, some additional enhancements are required for 3D printers to be effectively utilized for rock mechanics applications. In this study, we printed and studied specimens using a powder-based commercial ZPrinter® 450 with ZP® 150 powder and Zb® 63 binder used as raw materials. The specimens printed by this 3D printer exhibited relatively low strength and ductile behavior, implying that it needs further improvements. Hence, we focused on several ways to determine the best combination of printing options and post-processing including the effects of the printing direction, printing layer thickness, binder saturation level, and heating process on the uniaxial compressive strength (UCS) and stress-strain behavior of the printed samples. The suggested procedures have demonstrated their effectiveness by obtaining the printed samples that behave similarly to the natural rocks with low UCS. Although our optimization methods were particularly successful, further improvements are required to expand 3D printer application in the area of rock mechanics.

  7. Rock mechanics contributions from defense programs

    SciTech Connect

    Heuze, F.E.

    1992-02-01

    An attempt is made at illustrating the many contributions to rock mechanics from US defense programs, over the past 30-plus years. Large advances have been achieved in the technology-base area covering instrumentation, material properties, physical modeling, constitutive relations and numerical simulations. In the applications field, much progress has been made in understanding and being able to predict rock mass behavior related to underground explosions, cratering, projectile penetration, and defense nuclear waste storage. All these activities stand on their own merit as benefits to national security. But their impact is even broader, because they have found widespread applications in the non-defense sector; to name a few: the prediction of the response of underground structures to major earthquakes, the physics of the earth`s interior at great depths, instrumentation for monitoring mine blasting, thermo-mechanical instrumentation useful for civilian nuclear waste repositories, dynamic properties of earthquake faults, and transient large-strain numerical modeling of geological processes, such as diapirism. There is not pretense that this summary is exhaustive. It is meant to highlight success stories representative of DOE and DOD geotechnical activities, and to point to remaining challenges.

  8. Rock Mechanical Properties from Logs Petrophysics : Concepts and Results

    NASA Astrophysics Data System (ADS)

    Gaillot, Philippe; Crawford, Brian; Alramahi, Bashar; Karner, Steve

    2010-05-01

    The objective of the "geomechanics from logs" (GML) research project is to develop model-driven predictive software for determining rock mechanical properties (specifically rock strength, compressibility and fracability) from other, more easily measured, rock properties (e.g. lithology, porosity, clay volume, velocity) routinely derived from nuclear, resistivity and acoustic logging tools. To this end, geomechanics from logs seeks to increase fundamental understanding of the primary geologic controls on rock mechanical properties and to translate this new insight into novel predictive tools. In detail, GML predictors rely on (i) the generation of relational rock mechanical properties databases incorporating QC'd core-based laboratory measurements (both in-house and high-precision published data); (ii) the use of established rock physics models (e.g. friable sand, contact cement models) to investigate theoretical relationships between geologic processes, reservoir environment, rock microstructure and elastic, bulk and transport petrophysical attributes/properties; (iii) the subdivision of database rocks into generic lithotypes (e.g. sand, shaly sand, sandy shale, shale) with common petrophysical attributes/properties; (iv) the use of multivariate statistics to generate lithotype-dependent empirical predictive relationships between mechanical properties and log-derived petrophysical attributes/properties; (v) the estimation of uncertainties associated with predictive function parameters; (vi) the application and validation of mechanical properties predictive tools to well-documented case studies (e.g. sand strength for perforation stability, rock compressibility for reservoir simulation) to test overall performance and quantify uncertainty in predictions. This paper presents the results of various rock strength, rock compressibility and rock fracability case studies conducted in wells of different stratigraphic age and depositional environment. Overall, GML (i

  9. RATDAMPER - A Numerical Model for Coupling Mechanical and Hydrological Properties within the Disturbed Rock Zone at the Waste Isolation Pilot Plant

    SciTech Connect

    RATH,JONATHAN S.; PFEIFLE,T.W.; HUNSCHE,U.

    2000-11-27

    A numerical model for predicting damage and permeability in the disturbed rock zone (DRZ) has been developed. The semi-empirical model predicts damage based on a function of stress tensor invariant. For a wide class of problems hydrologic/mechanical coupling is necessary for proper analysis. The RATDAMPER model incorporates dilatant volumetric strain and permeability. The RATDAMPER model has been implemented in a weakly coupled code, which combines a finite element structural code and a finite difference multi-phase fluid flow code. Using the development of inelastic volumetric strain, a value of permeability can be assigned. This flexibility allows empirical permeability functional relationships to be evaluated.

  10. Integrated Experimental and Modeling Studies of Mineral Carbonation as a Mechanism for Permanent Carbon Sequestration in Mafic/Ultramafic Rocks

    SciTech Connect

    Wang, Zhengrong; Qiu, Lin; Zhang, Shuang; Bolton, Edward; Bercovici, David; Ague, Jay; Karato, Shun-Ichiro; Oristaglio, Michael; Zhu, Wen-Iu; Lisabeth, Harry; Johnson, Kevin

    2014-09-30

    A program of laboratory experiments, modeling and fieldwork was carried out at Yale University, University of Maryland, and University of Hawai‘i, under a DOE Award (DE-FE0004375) to study mineral carbonation as a practical method of geologic carbon sequestration. Mineral carbonation, also called carbon mineralization, is the conversion of (fluid) carbon dioxide into (solid) carbonate minerals in rocks, by way of naturally occurring chemical reactions. Mafic and ultramafic rocks, such as volcanic basalt, are natural candidates for carbonation, because the magnesium and iron silicate minerals in these rocks react with brines of dissolved carbon dioxide to form carbonate minerals. By trapping carbon dioxide (CO2) underground as a constituent of solid rock, carbonation of natural basalt formations would be a secure method of sequestering CO2 captured at power plants in efforts to mitigate climate change. Geochemical laboratory experiments at Yale, carried out in a batch reactor at 200°C and 150 bar (15 MPa), studied carbonation of the olivine mineral forsterite (Mg2SiO4) reacting with CO2 brines in the form of sodium bicarbonate (NaHCO3) solutions. The main carbonation product in these reactions is the carbonate mineral magnesite (MgCO3). A series of 32 runs varied the reaction time, the reactive surface area of olivine grains and powders, the concentration of the reacting fluid, and the starting ratio of fluid to olivine mass. These experiments were the first to study the rate of olivine carbonation under passive conditions approaching equilibrium. The results show that, in a simple batch reaction, olivine carbonation is fastest during the first 24 hours and then slows significantly and even reverses. A natural measure of the extent of carbonation is a quantity called the carbonation fraction, which compares the amount of carbon removed from solution, during a run, to the maximum amount

  11. Generalized Models for Rock Joint Surface Shapes

    PubMed Central

    Du, Shigui; Hu, Yunjin; Hu, Xiaofei

    2014-01-01

    Generalized models of joint surface shapes are the foundation for mechanism studies on the mechanical effects of rock joint surface shapes. Based on extensive field investigations of rock joint surface shapes, generalized models for three level shapes named macroscopic outline, surface undulating shape, and microcosmic roughness were established through statistical analyses of 20,078 rock joint surface profiles. The relative amplitude of profile curves was used as a borderline for the division of different level shapes. The study results show that the macroscopic outline has three basic features such as planar, arc-shaped, and stepped; the surface undulating shape has three basic features such as planar, undulating, and stepped; and the microcosmic roughness has two basic features such as smooth and rough. PMID:25152901

  12. Processes, mechanisms, parameters, and modeling approaches for partially saturated flow in soil and rock media; Yucca Mountain Site Characterization Project

    SciTech Connect

    Wang, J.S.Y.; Narasimhan, T.N.

    1993-06-01

    This report discusses conceptual models and mathematical equations, analyzes distributions and correlations among hydrological parameters of soils and tuff, introduces new path integration approaches, and outlines scaling procedures to model potential-driven fluid flow in heterogeneous media. To properly model the transition from fracture-dominated flow under saturated conditions to matrix-dominated flow under partially saturated conditions, characteristic curves and permeability functions for fractures and matrix need to be improved and validated. Couplings from two-phase flow, heat transfer, solute transport, and rock deformation to liquid flow are also important. For stochastic modeling of alternating units of welded and nonwelded tuff or formations bounded by fault zones, correlations and constraints on average values of saturated permeability and air entry scaling factor between different units need to be imposed to avoid unlikely combinations of parameters and predictions. Large-scale simulations require efficient and verifiable numerical algorithms. New path integration approaches based on postulates of minimum work and mass conservation to solve flow geometry and potential distribution simultaneously are introduced. This verifiable integral approach, together with fractal scaling procedures to generate statistical realizations with parameter distribution, correlation, and scaling taken into account, can be used to quantify uncertainties and generate the cumulative distribution function for groundwater travel times.

  13. Rock.XML - Towards a library of rock physics models

    NASA Astrophysics Data System (ADS)

    Jensen, Erling Hugo; Hauge, Ragnar; Ulvmoen, Marit; Johansen, Tor Arne; Drottning, Åsmund

    2016-08-01

    Rock physics modelling provides tools for correlating physical properties of rocks and their constituents to the geophysical observations we measure on a larger scale. Many different theoretical and empirical models exist, to cover the range of different types of rocks. However, upon reviewing these, we see that they are all built around a few main concepts. Based on this observation, we propose a format for digitally storing the specifications for rock physics models which we have named Rock.XML. It does not only contain data about the various constituents, but also the theories and how they are used to combine these building blocks to make a representative model for a particular rock. The format is based on the Extensible Markup Language XML, making it flexible enough to handle complex models as well as scalable towards extending it with new theories and models. This technology has great advantages as far as documenting and exchanging models in an unambiguous way between people and between software. Rock.XML can become a platform for creating a library of rock physics models; making them more accessible to everyone.

  14. Modeling Transport of Viruses in Fractured Rock

    NASA Astrophysics Data System (ADS)

    Sleep, B. E.; Mondal, P. K.

    2011-12-01

    Fractured rock aquifers are frequently used for water supply for human consumption. In many instances the fractured rock aquifers are vulnerable to contamination by pathogens, including viruses, due to co-location of on-site septic systems, wastewater discharges, biosolids and agricultural activities. Approximately half of the illnesses associated with groundwater consumption in the Unites States have been attributed to viral contamination. A number of these cases have been related to transport of viruses from septic systems to drinking water wells. Despite the potential for rapid transport of viruses through rock fractures to drinking water wells, the understanding of virus transport in fractured rock is limited. In particular, the impacts of virus size, fracture aperture variability and roughness, matrix porosity, groundwater velocity, and geochemical conditions have not been well studied. In this study, a multidimensional model for virus transport in variable aperture fractures is presented. The model is applied to laboratory experiments on transport of virus-sized latex microspheres (0.02 and 0.2 microns) and bacteriophages (MS2 and PR772) in artificially fractured dolomite rocks. In these experiments significant impacts of particle size, fracture characteristics, groundwater velocity, and geochemistry were observed. Given the variability in aperture distribution and associated spatial variation in groundwater flow field, one-dimensional models were not suitable for a comprehensive evaluation of the mechanisms governing the microsphere and bacteriophage transport. Various relationships for virus retention (attachment and detachment) are evaluated to provide insight into the governing processes in virus transport in fractured rock. In addition, the role of virus size, fracture aperture variability, fracture roughness, fracture surface charge, matrix porosity, groundwater velocity, and ionic strength in virus transport are evaluated. Scale-up to the field is

  15. Mechanical changes in thawing permafrost rocks and their influence on rock stability at the Zugspitze summit, Germany - a research concept

    NASA Astrophysics Data System (ADS)

    Mamot, Philipp; Scandroglio, Riccardo; Krautblatter, Michael

    2015-04-01

    strength tests as well as P-wave velocity measurements of dozens of frozen and unfrozen Zugspitze limestone samples. Our future tasks are as follows: i) To assess the spatial permafrost distribution in the slope we plan to conduct further laboratory-calibrated ERT and SRT. A dense rock temperature measuring network as well as nearby weather stations will supply input data for a simple thermal model of the rock slope. ii) To assess the spatial and temporal pattern of rock instability at the test site we will continue measuring discontinuity movements. iii) Undertaking rock-mechanical laboratory tests on Zugspitze limestone and to focus on temperature related changes of friction along rock discontinuities without ice infill and fracture toughness KIIc of intact rock bridges. These tests will be carried out with a direct shear box in unfrozen and frozen state. The measurement of P-wave velocity of the same rock samples will help to upscale rock toughness values to the rock slope at the study site. We aim at developing and calibrating a discontinuum mechanical model of stability changes in thawing permafrost rocks. Krautblatter, M., Funk, D., Günzel, F. K. (2013): Why permafrost rocks become unstable: a rock-ice-mechanical model in time and space. Earth Surface Processes and Landforms 38, 876-887.

  16. Mechanical and hydraulic properties of rocks related to induced seismicity

    USGS Publications Warehouse

    Witherspoon, P.A.; Gale, J.E.

    1977-01-01

    Witherspoon, P.A. and Gale, J.E., 1977. Mechanical and hydraulic properties of rocks related to induced seismicity. Eng. Geol., 11(1): 23-55. The mechanical and hydraulic properties of fractured rocks are considered with regard to the role they play in induced seismicity. In many cases, the mechanical properties of fractures determine the stability of a rock mass. The problems of sampling and testing these rock discontinuities and interpreting their non-linear behavior are reviewed. Stick slip has been proposed as the failure mechanism in earthquake events. Because of the complex interactions that are inherent in the mechanical behavior of fractured rocks, there seems to be no simple way to combine the deformation characteristics of several sets of fractures when there are significant perturbations of existing conditions. Thus, the more important fractures must be treated as individual components in the rock mass. In considering the hydraulic properties, it has been customary to treat a fracture as a parallel-plate conduit and a number of mathematical models of fracture systems have adopted this approach. Non-steady flow in fractured systems has usually been based on a two-porosity model, which assumes the primary (intergranular) porosity contributes only to storage and the secondary (fracture) porosity contributes only to the overall conductivity. Using such a model, it has been found that the time required to achieve quasi-steady state flow in a fractured reservoir is one or two orders of magnitude greater than it is in a homogeneous system. In essentially all of this work, the assumption has generally been made that the fractures are rigid. However, it is clear from a review of the mechanical and hydraulic properties that not only are fractures easily deformed but they constitute the main flow paths in many rock masses. This means that one must consider the interaction of mechanical and hydraulic effects. A considerable amount of laboratory and field data is now

  17. Some observations on the mechanism of aircraft wing rock

    NASA Technical Reports Server (NTRS)

    Hwang, C.; Pi, W. S.

    1978-01-01

    A pressure scale model of Northrop F-5A was tested in NASA Ames Research Center Eleven-Foot Transonic Tunnel to simulate the wing rock oscillations in a transonic maneuver. For this purpose, a flexible model support device was designed and fabricated which allowed the model to oscillate in roll at the scaled wing rock frequency. Two tunnel entries were performed to acquire the pressure (steady state and fluctuating) and response data when the model was held fixed and when it was excited by flow to oscillate in roll. Based on these data, a limit cycle mechanism was identified which supplied energy to the aircraft model and caused the Dutch roll type oscillations, commonly called wing rock. The major origin of the fluctuating pressures which contributed to the limit cycle was traced to the wing surface leading edge stall and the subsequent lift recovery. For typical wing rock oscillations, the energy balance between the pressure work input and the energy consumed by the model aerodynamic and mechanical damping was formulated and numerical data presented.

  18. Some observations on the mechanism of aircraft wing rock

    NASA Technical Reports Server (NTRS)

    Hwang, C.; Pi, W. S.

    1979-01-01

    A scale model of the Northrop F-5A was tested in NASA Ames Research Center Eleven-Foot Transonic Tunnel to simulate the wing rock oscillations in a transonic maneuver. For this purpose, a flexible model support device was designed and fabricated, which allowed the model to oscillate in roll at the scaled wing rock frequency. Two tunnel entries were performed to acquire the pressure (steady state and fluctuating) and response data when the model was held fixed and when it was excited by flow to oscillate in roll. Based on these data, a limit cycle mechanism was identified, which supplied energy to the aircraft model and caused the Dutch roll type oscillations, commonly called wing rock. The major origin of the fluctuating pressures that contributed to the limit cycle was traced to the wing surface leading edge stall and the subsequent lift recovery. For typical wing rock oscillations, the energy balance between the pressure work input and the energy consumed by the model's aerodynamic and mechanical damping was formulated and numerical data presented.

  19. Reconstruction of Sedimentary Rock Based on MechanicalProperties

    SciTech Connect

    Jin, Guodong; Patzek, Tad W.; Silin, Dmitry B.

    2004-05-04

    We describe a general, physics-based approach to numericalreconstruction of the geometrical structure and mechanical properties ofnatural sedimentary rock in 3D. Our procedure consists of three mainsteps: sedimentation, compaction, and diagenesis, followed by theverification of rock mechanical properties. The dynamic geologicprocesses of grain sedimentation and compaction are simulated by solvinga dimensionless form of Newton's equations of motion for an ensemble ofgrains. The diagenetic rock transformation is modeled using a cementationalgorithm, which accounts for the effect of rock grain size on therelative rate of cement overgrowth. Our emphasis is on unconsolidatedsand and sandstone. The main input parameters are the grain sizedistribution, the final rock porosity, the type and amount of cement andclay minerals, and grain mechanical properties: the inter-grain frictioncoefficient, the cement strength, and the grain stiffness moduli. We usea simulated 2D Fontainebleau sandstone to obtain the grain mechanicalproperties. This Fontainebleau sandstone is also used to study theinitiation, growth, and coalescence of micro-cracks under increasingvertical stress. The box fractal dimension of the micro-crackdistribution, and its variation with the applied stress areestimated.

  20. Laboratory rock mechanics testing manual. Public draft

    SciTech Connect

    Shuri, F S; Cooper, J D; Hamill, M L

    1981-10-01

    Standardized laboratory rock mechanics testing procedures have been prepared for use in the National Terminal Waste Storage Program. The procedures emphasize equipment performance specifications, documentation and reporting, and Quality Assurance acceptance criteria. Sufficient theoretical background is included to allow the user to perform the necessary data reduction. These procedures incorporate existing standards when possible, otherwise they represent the current state-of-the-art. Maximum flexibility in equipment design has been incorporated to allow use of this manual by existing groups and to encourage future improvements.

  1. Simulation of Electrical Transport in Rocks under Mechanical Action

    NASA Astrophysics Data System (ADS)

    Salgueiro da Silva, M. A.; Seixas, T. M.

    2015-12-01

    Rock's electrical properties can be changed by mechanical action, especially when deformation is accompanied by micro-fracturing processes. Knowing how electrical charge is generated in inelastically deformed rocks, the nature and properties of the generated charge carriers, and their spatial distribution and propagation is crucial to gain insight into the origin of seismo-electromagnetic signals. In this work, we describe briefly a model for the numerical simulation of electrical transport in rocks under mechanical action, assuming that high and low mobility charge carriers of opposite signs can be simultaneously generated by micro-fracturing processes and recombine, diffuse and drift across the sample rock. The electrical behavior can then be described using an adaptation of the formalism applied to semiconductors. We provide simulation results on a one-dimensional lattice using finite-difference discretization. Our results show that a large mobility contrast among charge carriers allows charge separation inside the deformation region, which leads to the formation of charged layers of alternate signs. Inside these layers, rapid electric field variations are observed which can lead to the emission of electromagnetic radiation. With proper positioning of current electrodes inside the deformation region, it is possible to collect electrical current even without any applied voltage. We discuss our results in the light of available experimental results on the generation of electrical and electromagnetic signals in deformed rocks.

  2. Investigations of Near-Field Thermal-Hydrologic-Mechanical-Chemical Models for Radioactive Waste Disposal in Clay/Shale Rock

    SciTech Connect

    Liu, H.H.; Li, L.; Zheng, L.; Houseworth, J.E.; Rutqvist, J.

    2011-06-20

    Clay/shale has been considered as potential host rock for geological disposal of high-level radioactive waste throughout the world, because of its low permeability, low diffusion coefficient, high retention capacity for radionuclides, and capability to self-seal fractures. For example, Callovo-Oxfordian argillites at the Bure site, France (Fouche et al., 2004), Toarcian argillites at the Tournemire site, France (Patriarche et al., 2004), Opalinus Clay at the Mont Terri site, Switzerland (Meier et al., 2000), and Boom clay at the Mol site, Belgium (Barnichon and Volckaert, 2003) have all been under intensive scientific investigation (at both field and laboratory scales) for understanding a variety of rock properties and their relationships to flow and transport processes associated with geological disposal of radioactive waste. Figure 1-1 presents the distribution of clay/shale formations within the USA.

  3. Lithophysal Rock Mass Mechanical Properties of the Repository Host Horizon

    SciTech Connect

    D. Rigby

    2004-11-10

    The purpose of this calculation is to develop estimates of key mechanical properties for the lithophysal rock masses of the Topopah Spring Tuff (Tpt) within the repository host horizon, including their uncertainties and spatial variability. The mechanical properties to be characterized include an elastic parameter, Young's modulus, and a strength parameter, uniaxial compressive strength. Since lithophysal porosity is used as a surrogate property to develop the distributions of the mechanical properties, an estimate of the distribution of lithophysal porosity is also developed. The resulting characterizations of rock parameters are important for supporting the subsurface design, developing the preclosure safety analysis, and assessing the postclosure performance of the repository (e.g., drift degradation and modeling of rockfall impacts on engineered barrier system components).

  4. High-pressure mechanical instability in rocks

    USGS Publications Warehouse

    Byerlee, J.D.; Brace, W.F.

    1969-01-01

    At a confining pressure of a few kilobars, deformation of many sedimentary rocks, altered mafic rocks, porous volcanic rocks, and sand is ductile, in that instabilities leading to audible elastic shocks are absent. At pressures of 7 to 10 kilobars, however, unstable faulting and stick-slip in certain of these rocks was observed. This high pressure-low temperature instability might be responsible for earthquakes in deeply buried sedimentary or volcanic sequences.

  5. Relevance of rock mechanics to gun propellants

    NASA Astrophysics Data System (ADS)

    Costantino, M.

    1983-10-01

    Over the past 30 years, earth scientists made a considerable effort to understand the mechanical response of geologic materials to large static and dynamic stresses. These materials range from loose soils to hard rock, with wide ranges of porosity, composition, particle size, and physical properties. The materials are described with fairly sophisticated equations of state that treat elastic and large strain plastic loading of porous, inhomogeneous, anisotropic bodies. Rate and temperature effects also are included, but are not as well understood. Well-developed laboratory methodology exists for studying compaction of porous bodies under hydrostatic loading, failure under shear loading, uniaxial strain loading, fluid flow through connected porosity, pressure dependence of sound speeds and moduli, and fracture of composite assemblies. Elaborate computer codes in one, two, and three dimensions are available to analyze these experiments.

  6. ROCK PROPERTIES MODEL ANALYSIS MODEL REPORT

    SciTech Connect

    Clinton Lum

    2002-02-04

    The purpose of this Analysis and Model Report (AMR) is to document Rock Properties Model (RPM) 3.1 with regard to input data, model methods, assumptions, uncertainties and limitations of model results, and qualification status of the model. The report also documents the differences between the current and previous versions and validation of the model. The rock properties models are intended principally for use as input to numerical physical-process modeling, such as of ground-water flow and/or radionuclide transport. The constraints, caveats, and limitations associated with this model are discussed in the appropriate text sections that follow. This work was conducted in accordance with the following planning documents: WA-0344, ''3-D Rock Properties Modeling for FY 1998'' (SNL 1997, WA-0358), ''3-D Rock Properties Modeling for FY 1999'' (SNL 1999), and the technical development plan, Rock Properties Model Version 3.1, (CRWMS M&O 1999c). The Interim Change Notice (ICNs), ICN 02 and ICN 03, of this AMR were prepared as part of activities being conducted under the Technical Work Plan, TWP-NBS-GS-000003, ''Technical Work Plan for the Integrated Site Model, Process Model Report, Revision 01'' (CRWMS M&O 2000b). The purpose of ICN 03 is to record changes in data input status due to data qualification and verification activities. These work plans describe the scope, objectives, tasks, methodology, and implementing procedures for model construction. The constraints, caveats, and limitations associated with this model are discussed in the appropriate text sections that follow. The work scope for this activity consists of the following: (1) Conversion of the input data (laboratory measured porosity data, x-ray diffraction mineralogy, petrophysical calculations of bound water, and petrophysical calculations of porosity) for each borehole into stratigraphic coordinates; (2) Re-sampling and merging of data sets; (3) Development of geostatistical simulations of porosity; (4

  7. MECHANICAL DEGRADATION OF EMPLACEMENT DRIFTS AT YUCCA MOUNTAIN - A CASE STUDY IN ROCK MECHANICS, PART 1: NONLITHOPHYSAL ROCK, PART 2: LITHOPHYSAL ROCK

    SciTech Connect

    M. Lin, D. Kicker, B. Damjanac, M. Board, and M. Karakouzian

    2006-02-27

    This paper outlines rock mechanics investigations associated with mechanical degradation of planned emplacement drifts at Yucca Mountain, which is the designated site for a US high-level nuclear waste repository. The factors leading to drift degradation include stresses from the overburden, stresses induced by the heat released from the emplaced waste, stresses due to seismically related ground motions, and time-dependent strength degradation. The welded tuff emplacement horizon consists of two groups of rock with distinct engineering properties: nonlithophysal units and lithophysal units, based on the relative proportion of lithophysal cavities. Part I of the paper concentrates on the generally hard, strong, and fractured nonlithophysal rock. The degradation behavior of the tunnels in the nonlithophysal rock is controlled by the occurrence of keyblocks. A statistically equivalent fracture model was generated based on extensive underground fracture mapping data from the Exploratory Studies Facility at Yucca Mountain. Three-dimensional distinct block analyses, generated with the fracture patterns randomly selected from the fracture model, were developed with the consideration of in situ, thermal, seismic loads. In this study, field data, laboratory data, and numerical analyses are well integrated to provide a solution for the unique problem of modeling drift degradation throughout the regulatory period for repository performance.

  8. Influence of sedimentary environments on mechanical properties of clastic rocks

    NASA Astrophysics Data System (ADS)

    Meng, Zhaoping; Zhang, Jincai; Peng, Suping

    2006-10-01

    The sedimentary environments are the intrinsic factor controlling the mechanical properties of clastic rocks. Examining the relationship between rock sedimentary environments and rock mechanical properties gives a better understanding of rock deformation and failure mechanisms. In this study, more than 55 samples in coal measures were taken from seven different lithologic formations in eastern China. Using the optical microscope the sedimentary characteristics, such as components of clastic rocks and sizes of clastic grains were quantitatively tested and analyzed. The corresponding mechanical parameters were tested using the servo-controlled testing system. Different lithologic attributes in the sedimentary rocks sampled different stress-strain behaviors and failure characteristics under different confining pressures, mainly due to different compositions and textures. Results demonstrate that clastic rocks have the linear best-fit for Mohr-Coulomb failure criterion. The elastic moduli in clastic rocks are highly dependent upon confining pressures, unlike hard rocks. The envelope lines of the mechanical properties versus the contents of quartz, detritus of the grain diameter of more than 0.03 mm, and grain size in clastic rocks are given. The compressive strength or elastic modulus and the grain diameter have a non-monotonic relation and demonstrate the “grain-diameter softening” effect.

  9. Carbonate rock depositional models: A microfacies approach

    SciTech Connect

    Carozzi, A.V.

    1988-01-01

    Carbonate rocks contain more than 50% by weight carbonate minerals such as calcite, dolomite, and siderite. Understanding how these rocks form can lead to more efficient methods of petroleum exploration. Micofacies analysis techniques can be used as a method of predicting models of sedimentation for carbonate rocks. Micofacies in carbonate rocks can be seen clearly only in thin sections under a microscope. This section analysis of carbonate rocks is a tool that can be used to understand depositional environments, diagenetic evolution of carbonate rocks, and the formation of porosity and permeability in carbonate rocks. The use of micofacies analysis techniques is applied to understanding the origin and formation of carbonate ramps, carbonate platforms, and carbonate slopes and basins. This book will be of interest to students and professionals concerned with the disciplines of sedimentary petrology, sedimentology, petroleum geology, and palentology.

  10. Mechanical and acoustic properties of weakly cemented granular rocks

    SciTech Connect

    Nakagawa, S.; Myer, L.R.

    2001-05-09

    This paper presents the results of laboratory measurements on the mechanical and acoustic properties of weakly cemented granular rock. Artificial rock samples were fabricated by cementing sand and glass beads with sodium silicate binder. During uniaxial compression tests, the rock samples showed stress-strain behavior which was more similar to that of soils than competent rocks, exhibiting large permanent deformations with frictional slip. The mechanical behavior of the samples approached that of competent rocks as the amount of binder was increased. For very weak samples, acoustic waves propagating in these rocks showed very low velocities of less than 1000 m/sec for compressional waves. A borehole made within this weakly cemented rock exhibited a unique mode of failure that is called ''anti-KI mode fracture'' in this paper. The effect of cementation, grain type, and boundary conditions on this mode of failure was also examined experimentally.

  11. A Refined Model for Solid Particle Rock Erosion

    NASA Astrophysics Data System (ADS)

    Momber, A. W.

    2016-02-01

    A procedure for the estimation of distribution parameters of a Weibull distribution model K 1 = f( K Ic 12/4 / σ C 23/4 ) for solid particle erosion, as recently suggested in Rock Mech Rock Eng, doi: 10.1007/s00603-014-0658-x, 2014, is derived. The procedure is based on examinations of elastic-plastically responding rocks (rhyolite, granite) and plastically responding rocks (limestone, schist). The types of response are quantified through SEM inspections of eroded surfaces. Quantitative numbers for the distribution parameter K 1 are calculated for 30 rock materials, which cover a wide range of mechanical properties. The ranking according to the parameter K 1 is related to qualitative rock classification schemes. A modified proposal for the erosion of schist due to solid particle impingement at normal incidence is introduced.

  12. An Elastoplastic Model for Partially Saturated Collapsible Rocks

    NASA Astrophysics Data System (ADS)

    Ma, Jianjun

    2016-02-01

    A unified elastoplastic model for describing the stress-strain behavior of partially saturated collapsible rocks is proposed. The elastic-plastic response due to loading and unloading is captured using bounding surface plasticity. The coupling effect of hydraulic and mechanical responses is addressed by applying the effective stress concept. Special attention is paid to the rock-fluid characteristic curve (RFCC), effective stress parameter, and suction hardening. A wide range of saturation degree is considered. The characteristics of mechanical behavior in partially saturated collapsible rocks are captured for all cases considered.

  13. Laboratory modeling of seismoelectric effects in rock

    NASA Astrophysics Data System (ADS)

    Besedina, Alina; Kocharyan, Gevorg

    2010-05-01

    It is well known that deformation of rock by seismic waves is accompanied by a complex of various electromagnetic effects of different physical natures. These phenomena are widely used in exploration geophysics and well tests, besides that these effects are of great interest in forecasting catastrophic events, such as rock bursts in mining operations. For adequate interpretation of experimental results it is necessary to understand the physical nature of the seismoelectric effect. Despite of a considerable amount of performed investigations, no general model of the phenomenon has been developed yet. The known sources of electric signals in rock are electrokinetic phenomena, piezoelectric phenomena, triboelectricity, contact electrification, induction phenomena and the effect of charged edge dislocation oscillations. One of the urgent questions is studying the relationships between form and amplitude of the seismic pulse and the electric signal. In this work an experimental investigation is presented of the process of electric signal origination in hard rock, which does not contain fluid in an explicit form. The constructed laboratory set-up allows to make experiments with compressional waves of a wide range of amplitudes and frequencies. It also allows to simulate both continuous media, and fractured rock. Marble, granite and a model material made of hyposulphite mixed with granite crumb were used in this research. Longitudinal waves of different intensities were initiated in the model by impacts of balls of different masses. The constructed one-dimensional model - the rod - provides conditions for formation and propagation of a plane wave as well as a noticeable delay of the arrival of the tension wave reflected from the free end of the rod. This permitted to sort out clearly the electric signals accompanying propagation of a longitudinal compressional wave and to find out the degree of correlation between the parameters of electric and mechanical signals. It is

  14. Interaction of thermal and mechanical processes in steep permafrost rock walls: A conceptual approach

    NASA Astrophysics Data System (ADS)

    Draebing, D.; Krautblatter, M.; Dikau, R.

    2014-12-01

    Degradation of permafrost rock wall decreases stability and can initiate rock slope instability of all magnitudes. Rock instability is controlled by the balance of shear forces and shear resistances. The sensitivity of slope stability to warming results from a complex interplay of shear forces and resistances. Conductive, convective and advective heat transport processes act to warm, degrade and thaw permafrost in rock walls. On a seasonal scale, snow cover changes are a poorly understood key control of the timing and extent of thawing and permafrost degradation. We identified two potential critical time windows where shear forces might exceed shear resistances of the rock. In early summer combined hydrostatic and cryostatic pressure can cause a peak in shear force exceeding high frozen shear resistance and in autumn fast increasing shear forces can exceed slower increasing shear resistance. On a multiannual system scale, shear resistances change from predominantly rock-mechanically to ice-mechanically controlled. Progressive rock bridge failure results in an increase of sensitivity to warming. Climate change alters snow cover and duration and, hereby, thermal and mechanical processes in the rock wall. Amplified thawing of permafrost will result in higher rock slope instability and rock fall activity. We present a holistic conceptual approach connecting thermal and mechanical processes, validate parts of the model with geophysical and kinematic data and develop future scenarios to enhance understanding on system scale.

  15. Prediction of Fracture Behavior in Rock and Rock-like Materials Using Discrete Element Models

    NASA Astrophysics Data System (ADS)

    Katsaga, T.; Young, P.

    2009-05-01

    The study of fracture initiation and propagation in heterogeneous materials such as rock and rock-like materials are of principal interest in the field of rock mechanics and rock engineering. It is crucial to study and investigate failure prediction and safety measures in civil and mining structures. Our work offers a practical approach to predict fracture behaviour using discrete element models. In this approach, the microstructures of materials are presented through the combination of clusters of bonded particles with different inter-cluster particle and bond properties, and intra-cluster bond properties. The geometry of clusters is transferred from information available from thin sections, computed tomography (CT) images and other visual presentation of the modeled material using customized AutoCAD built-in dialog- based Visual Basic Application. Exact microstructures of the tested sample, including fractures, faults, inclusions and void spaces can be duplicated in the discrete element models. Although the microstructural fabrics of rocks and rock-like structures may have different scale, fracture formation and propagation through these materials are alike and will follow similar mechanics. Synthetic material provides an excellent condition for validating the modelling approaches, as fracture behaviours are known with the well-defined composite's properties. Calibration of the macro-properties of matrix material and inclusions (aggregates), were followed with the overall mechanical material responses calibration by adjusting the interfacial properties. The discrete element model predicted similar fracture propagation features and path as that of the real sample material. The path of the fractures and matrix-inclusion interaction was compared using computed tomography images. Initiation and fracture formation in the model and real material were compared using Acoustic Emission data. Analysing the temporal and spatial evolution of AE events, collected during the

  16. Modelling Fracture Propagation in Anisotropic Rock Mass

    NASA Astrophysics Data System (ADS)

    Shen, Baotang; Siren, Topias; Rinne, Mikael

    2015-05-01

    Anisotropic rock mass is often encountered in rock engineering, and cannot be simplified as an isotropic problem in numerical models. A good understanding of rock fracturing processes and the ability to predict fracture initiation and propagation in anisotropic rock masses are required for many rock engineering problems. This paper describes the development of the anisotropic function in FRACOD—a specialized fracture propagation modelling software—and its recent applications to rock engineering issues. Rock anisotropy includes strength anisotropy and modulus anisotropy. The level of complexity in developing the anisotropic function for strength anisotropy and modulus anisotropy in FRACOD is significantly different. The strength anisotropy function alone does not require any alteration in the way that FRACOD calculates rock stress and displacement, and therefore is relatively straightforward. The modulus anisotropy function, on the other hand, requires modification of the fundamental equations of stress and displacement in FRACOD, a boundary element code, and hence is more complex and difficult. In actual rock engineering, the strength anisotropy is often considered to be more pronounced and important than the modulus anisotropy, and dominates the stability and failure pattern of the rock mass. The modulus anisotropy will not be considered in this study. This paper discusses work related to the development of the strength anisotropy in FRACOD. The anisotropy function has been tested using numerical examples. The predicted failure surfaces are mostly along the weakest planes. Predictive modelling of the Posiva's Olkiluoto Spalling Experiment was made. The model suggests that spalling is very sensitive to the direction of anisotropy. Recent observations from the in situ experiment showed that shear fractures rather than tensile fractures occur in the holes. According to the simulation, the maximum tensile stress is well below the tensile strength, but the maximum

  17. Analyzing failure modes of rock mass based on statistical mechanics of rock mass

    NASA Astrophysics Data System (ADS)

    Bao, H.; Wu, F.

    2015-12-01

    Joints influence mechanical properties of rock mass. Based on the strength criterion of statistical mechanics of rock mass, we analyzed the four different failure modes of rock mass with a group of joints by combining with Mohr-Coulomb strength criterion. We also deduced an expression of the critical confining pressure for explaining the transformation from structure control to stress control of rock mass strength. On this basis, rock mass with a group of joints were divided into four classes according to the relations between rock mass and joints parameters. Then, the possible failure modes and their corresponding conditions were discussed. At last, the strength characteristics of diorite with a group of joints were analyzed. The results showed that the diorite belonged to class rock mass and performed significant anisotropy in compressive strength. At the condition of vertical pressure, the rock start failing after joints at the critical confining pressure of 9.12MPa. However, with the confining pressure increasing, the anisotropy of strength became weak, and the strength of diorite would convert from structure control to stress control under some particular loading directions.

  18. Probing Mechanical Properties of Rock with InSAR

    NASA Astrophysics Data System (ADS)

    Jónsson, S.

    2012-04-01

    Interferometric Synthetic Aperture Radar (InSAR) observations from satellites have revolutionized our crustal deformation measurement capabilities with its high spatial resolution, global coverage, and low cost. The high spatial resolution (typically 5-20 m) allows us to map many small-scale surface deformation phenomena in great detail. These include surface faulting, fissuring, fault creep, and other strain localization phenomena. Another advantage of the small-scale deformation mapping is that it can provide information about mechanical properties of near-surface rocks. Several studies have already been published on using InSAR to probe material properties of rock. Strain localizations at fault zones have been observed in co-seismic deformation fields near to large earthquakes and interpreted as expressions of weak fault zone materials that are a factor of two more compliant than the surrounding unbroken rock [Fialko et al., 2002]. Peltzer et al. [1999] argued that asymmetries in coseismic deformation patterns observed by InSAR showed evidence for non-linear elasticity, i.e. that the elastic moduli of shallow crustal material are different for compression and extension, due to small-scale cracks in the medium. This interpretation was later disputed by Funning et al. [2007], who provided an alternative explanation for observed deformation pattern based on along-strike variations in fault geometry and slip. In addition, observations and modeling of poro-elastic rebound after earthquakes have provided information about the difference in undrained and drained Poisson's ratio values of the near-surface rocks [Peltzer et al., 1996; Jónsson et al., 2003]. More recently we have used InSAR observations to put bounds on the tensional bulk strength of surface rocks. A dyke intrusion that took place in western Saudi Arabia in 2009 caused many moderate-sized earthquakes and extensive surface faulting. InSAR data of the area show that large-scale (40 km x 40 km) east

  19. Modeling of nuclear waste disposal by rock melting

    SciTech Connect

    Heuze, F.E.

    1982-04-01

    Today, the favored option for disposal of high-level nuclear wastes is their burial in mined caverns. As an alternative, the concept of deep disposal by rock melting (DRM) also has received some attention. DRM entails the injection of waste, in a cavity or borehole, 2 to 3 kilometers down in the earth crust. Granitic rocks are the prime candidate medium. The high thermal loading initially will melt the rock surrounding the waste. Following resolidification, a rock/waste matrix is formed, which should provide isolation for many years. The complex thermal, mechanical, and hydraulic aspects of DRM can be studied best by means of numerical models. The models must accommodate the coupling of the physical processes involved, and the temperature dependency of the granite properties, some of which are subject to abrupt discontinuities, during ..cap alpha..-..beta.. phase transition and melting. This paper outlines a strategy for such complex modeling.

  20. Rock mechanics issues in completion and stimulation operations

    SciTech Connect

    Warpinski, N.R.

    1992-02-01

    Rock mechanisms parameters such as the in situ stresses, elastic properties, failure characteristics, and poro-elastic response are important to most completion and stimulation operations. Perforating, hydraulic fracturing, wellbore stability, and sand production are examples of technology that are largely controlled by the rock mechanics of the process. While much research has been performed in these areas, there has been insufficient application that research by industry. In addition, there are new research needs that must be addressed for technology advancement.

  1. Rock mechanics issues in completion and stimulation operations

    SciTech Connect

    Warpinski, N.R.

    1992-01-01

    Rock mechanisms parameters such as the in situ stresses, elastic properties, failure characteristics, and poro-elastic response are important to most completion and stimulation operations. Perforating, hydraulic fracturing, wellbore stability, and sand production are examples of technology that are largely controlled by the rock mechanics of the process. While much research has been performed in these areas, there has been insufficient application that research by industry. In addition, there are new research needs that must be addressed for technology advancement.

  2. Significance of grain sliding mechanisms for ductile deformation of rocks

    NASA Astrophysics Data System (ADS)

    Dimanov, A.; Bourcier, M.; Gaye, A.; Héripré, E.; Bornert, M.; Raphanel, J.; Ludwig, W.

    2013-12-01

    Ductile shear zones at depth present polyphase and heterogeneous rocks and multi-scale strain localization patterns. Most strain concentrates in ultramylonitic layers, which exhibit microstructural signatures of several concomitant deformation mechanisms. The latter are either active in volume (dislocation creep), or in the vicinity and along interfaces (grain sliding and solution mass transfer). Because their chronology of appearance and interactions are unclear, inference of the overall rheology seems illusory. We have therefore characterized over a decade the rheology of synthetic lower crustal materials with different compositions and fluid contents, and for various microstructures. Non-Newtonian flow clearly related to dominant dislocation creep. Conversely, Newtonian behavior involved grain sliding mechanisms, but crystal plasticity could be identified as well. In order to clarify the respective roles of these mechanisms we underwent a multi-scale investigation of the ductile deformation of rock analog synthetic halite with controlled microstructures. The mechanical tests were combined with in-situ optical microscopy, scanning electron microscopy and X ray computed tomography, allowing for digital image correlation (DIC) techniques and retrieval of full strain field. Crystal plasticity dominated, as evidenced by physical slip lines and DIC computed slip bands. Crystal orientation mapping allowed to identify strongly active easy glide {110} <110> systems. But, all other slip systems were observed as well, and especially near interfaces, where their activity is necessary to accommodate for the plastic strain incompatibilities between neighboring grains. We also evidenced grain boundary sliding (GBS), which clearly occurred as a secondary, but necessary, accommodation mechanism. The DIC technique allowed the quantification of the relative contribution of each mechanism. The amount of GBS clearly increased with decreasing grain size. Finite element (FE) modeling

  3. Exact effective-stress rules in rock mechanics

    SciTech Connect

    Berryman, J.G. )

    1992-09-15

    The standard paradigm for analysis of rock deformation arises from postulating the existence of an equivalent homogeneous porous rock.'' However, data on the pore-pressure dependence of fluid permeability for some rocks cannot be explained using any equivalent homogeneous porous medium. In contrast, a positive result shows that deformation measurements on both high-porosity sandstones and low-porosity granites can be explained adequately in terms of an equivalent two-constituent model of porous rocks, for which exact results have recently been discovered.

  4. Conference addresses thermo-hydro-mechanical coupling in fractured rock

    NASA Astrophysics Data System (ADS)

    Kümpel, Hans-Joachim

    Various environmental problems and the use of certain energy resources are closely related to fluid flow in and the mechanical behavior of porous or fractured rock. Subjects of obvious socioeconomic relevance are the supply and protection of groundwater, the production of hydrocarbon reservoirs, land subsidence in coastal areas, exploitation of geothermal energy and the long-term disposal of critical wastes. Efficient management of such issues is often hampered by the fact that rocks and rock formations are inherently complex. Any rock sample is an aggregate of the myriad mineral particles forming its matrix and fluid molecules residing in voids. Any two rock samples differ in many aspects, including geochemical constituents, size and shape of grains, structure of pore space, and fracture networks.

  5. Chemically- and mechanically-mediated influences on the transport and mechanical characteristics of rock fractures

    SciTech Connect

    Min, K.-B.; Rutqvist, J.; Elsworth, D.

    2009-02-01

    A model is presented to represent changes in the mechanical and transport characteristics of fractured rock that result from coupled mechanical and chemical effects. The specific influence is the elevation of dissolution rates on contacting asperities, which results in a stress- and temperature-dependent permanent closure. A model representing this pressure-dissolution-like behavior is adapted to define the threshold and resulting response in terms of fundamental thermodynamic properties of a contacting fracture. These relations are incorporated in a stress-stiffening model of fracture closure to define the stress- and temperature-dependency of aperture loss and behavior during stress and temperature cycling. These models compare well with laboratory and field experiments, representing both decoupled isobaric and isothermal responses. The model was applied to explore the impact of these responses on heated structures in rock. The result showed a reduction in ultimate induced stresses over the case where chemical effects were not incorporated, with permanent reduction in final stresses after cooling to ambient conditions. Similarly, permeabilities may be lower than they were in the case where chemical effects were not considered, with a net reduction apparent even after cooling to ambient temperature. These heretofore-neglected effects may have a correspondingly significant impact on the performance of heated structures in rock, such as repositories for the containment of radioactive wastes.

  6. Rock mechanics: Issues and research needs in the disposal of wastes in hydraulic fractures

    NASA Astrophysics Data System (ADS)

    Doe, T.; McClain, W. C.

    1984-07-01

    Proposed rock mechanics studies are outlined which are designed to answer basic questions concerning hydraulic fracturing for waste disposal. These questions are: (1) how can containment be assured for Oak Ridge or other sites; and (2) what is the capacity of a site. The suggested rock mechanics program consists of four major tasks: (1) numerical modeling; (2) laboratory testing; (3) field testing; and (4) monitoring. These tasks are described.

  7. Rock mechanics issues and research needs in the disposal of wastes in hydraulic fractures

    SciTech Connect

    Doe, T.W.; McClain, W.C.

    1984-07-01

    The proposed rock mechanics studies outlined in this document are designed to answer the basic questions concerning hydraulic fracturing for waste disposal. These questions are: (1) how can containment be assured for Oak Ridge or other sites; and (2) what is the capacity of a site. The suggested rock mechanics program consists of four major tasks: (1) numerical modeling, (2) laboratory testing, (3) field testing, and (4) monitoring. These tasks are described.

  8. Mechanical Properties of Shock-Damaged Rocks

    NASA Technical Reports Server (NTRS)

    He, Hongliang; Ahrens, T. J.

    1994-01-01

    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.

  9. Rock mechanics investigations, design and construction of the Ridracoli dam

    NASA Astrophysics Data System (ADS)

    Oberti, G.; Bavestrello, F.; Rossi, P. P.; Flamigni, F.

    1986-07-01

    The Ridracoli arch-gravity concrete dam is the major work in a multipurpose project whose main scope is the water supply to 37 communities in the Forlì and Ravenna Provinces. The particular geological and structural characteristics of the foundation rock mass, consisting of a rhythmical alternation of sandstone, silstone and marl, required a wide program of in situ and laboratory investigations in order to obtain a detailed physical and mechanical characterization of the foundation. The design criteria based on the use of a physical and mathematical model are illustrated as well as the limit equilibrium analysis of the stability conditions of the abutments. Particular attention has been devoted to the problem of the excavation stability; excavation methods and stabilizing works are illustrated in detail. The scepage problems are also presented with the description of grouting and drainage works.

  10. New Mechanisms of rock-bit wear in geothermal wells

    SciTech Connect

    Macini, Paolo

    1996-01-24

    This paper presents recent results of an investigation on failure mode and wear of rock-bits used to drill geothermal wells located in the area of Larderello (Italy). A new wear mechanism, conceived from drilling records and dull bit evaluation analysis, has been identified and a particular configuration of rock-bit has been developed and tested in order to reduce drilling costs. The role of high Bottom Hole Temperature (BHT) on rock-bit performances seems not yet very well understood: so far, only drillability and formation abrasiveness are generally considered to account for poor drilling performances. In this paper, the detrimental effects of high BHT on sealing and reservoir system of Friction Bearing Rock-bits (FBR) have been investigated, and a new bearing wear pattern for FBR's run in high BHT holes has been identified and further verified via laboratory inspections on dull bits. A novel interpretation of flat worn cutting structure has been derived from the above wear pattern, suggesting the design of a particular bit configuration. Test bits, designed in the light of the above criteria, have been prepared and field tested successfully. The paper reports the results of these tests, which yielded a new rock-bit application, today considered as a standad practice in Italian geothermal fields. This application suggests that the correct evaluation of rock-bit wear can help to improve the overall drilling performances and to minimize drilling problems through a better interpretation of the relationships amongst rock-bits, formation properties and downhole temperature.

  11. Stochastic multiscale model for carbonate rocks.

    PubMed

    Biswal, B; Oren, P-E; Held, R J; Bakke, S; Hilfer, R

    2007-06-01

    A multiscale model for the diagenesis of carbonate rocks is proposed. It captures important pore scale characteristics of carbonate rocks: wide range of length scales in the pore diameters; large variability in the permeability; and strong dependence of the geometrical and transport parameters on the resolution. A pore scale microstructure of an oolithic dolostone with generic diagenetic features is successfully generated. The continuum representation of a reconstructed cubic sample of side length 2mm contains roughly 42 x 10{6} crystallites and pore diameters varying over many decades. Petrophysical parameters are computed on discretized samples of sizes up to 1000{3}. The model can be easily adapted to represent the multiscale microstructure of a wide variety of carbonate rocks. PMID:17677251

  12. Experimental Studies on the Mechanical Behaviour of Rock Joints with Various Openings

    NASA Astrophysics Data System (ADS)

    Li, Y.; Oh, J.; Mitra, R.; Hebblewhite, B.

    2016-03-01

    The mechanical behaviour of rough joints is markedly affected by the degree of joint opening. A systematic experimental study was conducted to investigate the effect of the initial opening on both normal and shear deformations of rock joints. Two types of joints with triangular asperities were produced in the laboratory and subjected to compression tests and direct shear tests with different initial opening values. The results showed that opened rock joints allow much greater normal closure and result in much lower normal stiffness. A semi-logarithmic law incorporating the degree of interlocking is proposed to describe the normal deformation of opened rock joints. The proposed equation agrees well with the experimental results. Additionally, the results of direct shear tests demonstrated that shear strength and dilation are reduced because of reduced involvement of and increased damage to asperities in the process of shearing. The results indicate that constitutive models of rock joints that consider the true asperity contact area can be used to predict shear resistance along opened rock joints. Because rock masses are loosened and rock joints become open after excavation, the model suggested in this study can be incorporated into numerical procedures such as finite-element or discrete-element methods. Use of the model could then increase the accuracy and reliability of stability predictions for rock masses under excavation.

  13. Thermo-hydro-mechanical processes in fractured rock formations during a glacial advance

    NASA Astrophysics Data System (ADS)

    Selvadurai, A. P. S.; Suvorov, A. P.; Selvadurai, P. A.

    2015-07-01

    The paper examines the coupled thermo-hydro-mechanical (THM) processes that develop in a fractured rock region within a fluid-saturated rock mass due to loads imposed by an advancing glacier. This scenario needs to be examined in order to assess the suitability of potential sites for the location of deep geologic repositories for the storage of high-level nuclear waste. The THM processes are examined using a computational multiphysics approach that takes into account thermo-poroelasticity of the intact geological formation and the presence of a system of sessile but hydraulically interacting fractures (fracture zones). The modelling considers coupled thermo-hydro-mechanical effects in both the intact rock and the fracture zones due to contact normal stresses and fluid pressure at the base of the advancing glacier. Computational modelling provides an assessment of the role of fractures in modifying the pore pressure generation within the entire rock mass.

  14. Thermo-hydro-mechanical processes in fractured rock formations during glacial advance

    NASA Astrophysics Data System (ADS)

    Selvadurai, A. P. S.; Suvorov, A. P.; Selvadurai, P. A.

    2014-11-01

    The paper examines the coupled thermo-hydro-mechanical (THM) processes that develop in a fractured rock region within a fluid-saturated rock mass due to loads imposed by an advancing glacier. This scenario needs to be examined in order to assess the suitability of potential sites for the location of deep geologic repositories for the storage of high-level nuclear waste. The THM processes are examined using a computational multiphysics approach that takes into account thermo-poroelasticity of the intact geological formation and the presence of a system of sessile but hydraulically interacting fractures (fracture zones). The modeling considers coupled thermo-hydro-mechanical effects in both the intact rock and the fracture zones due to contact normal stresses and fluid pressure at the base of the advancing glacier. Computational modelling provides an assessment of the role of fractures that can modify the pore pressure generation within the entire rock mass.

  15. Interactive evolution concept for analyzing a rock salt cavern under cyclic thermo-mechanical loading

    NASA Astrophysics Data System (ADS)

    König, Diethard; Mahmoudi, Elham; Khaledi, Kavan; von Blumenthal, Achim; Schanz, Tom

    2016-04-01

    The excess electricity produced by renewable energy sources available during off-peak periods of consumption can be used e.g. to produce and compress hydrogen or to compress air. Afterwards the pressurized gas is stored in the rock salt cavities. During this process, thermo-mechanical cyclic loading is applied to the rock salt surrounding the cavern. Compared to the operation of conventional storage caverns in rock salt the frequencies of filling and discharging cycles and therefore the thermo-mechanical loading cycles are much higher, e.g. daily or weekly compared to seasonally or yearly. The stress strain behavior of rock salt as well as the deformation behavior and the stability of caverns in rock salt under such loading conditions are unknown. To overcome this, existing experimental studies have to be supplemented by exploring the behavior of rock salt under combined thermo-mechanical cyclic loading. Existing constitutive relations have to be extended to cover degradation of rock salt under thermo-mechanical cyclic loading. At least the complex system of a cavern in rock salt under these loading conditions has to be analyzed by numerical modeling taking into account the uncertainties due to limited access in large depth to investigate material composition and properties. An interactive evolution concept is presented to link the different components of such a study - experimental modeling, constitutive modeling and numerical modeling. A triaxial experimental setup is designed to characterize the cyclic thermo-mechanical behavior of rock salt. The imposed boundary conditions in the experimental setup are assumed to be similar to the stress state obtained from a full-scale numerical simulation. The computational model relies primarily on the governing constitutive model for predicting the behavior of rock salt cavity. Hence, a sophisticated elasto-viscoplastic creep constitutive model is developed to take into account the dilatancy and damage progress, as well as

  16. Decompaction mechanism of deep crystalline rocks under stress relief

    NASA Astrophysics Data System (ADS)

    Gorbatsevich, F. F.

    2003-07-01

    Within a geological massif in a stable geodynamical situation contacts on the grain boundaries in polycrystalline rocks at great depths are continuous and firm. The stress release of those rocks during drilling and excavation to the surface is accompanied by their disintegration (decompaction). The reason for the decompaction is generation of microcracks during stress release due to the difference between the elastic moduli of crystalline grains at their contacts. The mechanism of decompaction may occur not only in polymineral but in polycrystalline rocks as well. The method of decompaction evaluation of deep crystalline rocks under stress relief is presented. According to the calculations the initial manifestation of the decompaction effect in biotite gneisses will occur when they are extracted from the deep range of 0.8-1 km. The first microcracks arise on the grain borders between quartz-biotite and oligoclase-biotite. It is shown that the uplift of gneiss-granite varities of the rocks cut by the Kola superdeep borehole from depths exceeding 13-15 km will be possible in a form of separate mineral grains. Practical importance of the presented method is in an opportunity to evaluate the level of excavated decompaction. The method allow estimating the depth, from which the rock will be extracted only in a sludge form.

  17. Modeling capillary barriers in unsaturated fractured rock

    NASA Astrophysics Data System (ADS)

    Wu, Yu-Shu; Zhang, W.; Pan, Lehua; Hinds, Jennifer; Bodvarsson, G. S.

    2002-11-01

    This work presents a series of numerical modeling studies that investigate the hydrogeologic conditions required to form capillary barriers and the effect that capillary barriers have on fluid flow and tracer transport processes in the unsaturated fractured rock of Yucca Mountain, Nevada, a potential site for storing high-level radioactive waste. The modeling approach is based on a dual-continuum formulation of coupled multiphase fluid and tracer transport through fractured porous rock. The numerical modeling results showed that effective capillary barriers can develop where both matrix and fracture capillary gradients tend to move water upward. Under the current hydrogeologic conceptualization of Yucca Mountain, strong capillary barrier effects exist for diverting a significant amount of moisture flow through the relatively shallow Paintbrush nonwelded unit, with major faults observed at the site serving as major downward pathways for laterally diverted percolation fluxes. In addition, we used observed field liquid saturation and goechemical isotopic data to check model results and found consistent agreement.

  18. A constitutive model for predicting rock fragmentation by blasting

    NASA Astrophysics Data System (ADS)

    Liu, Liqing; Katsabanis, P. D.

    1996-05-01

    This paper describes the development of a constitutive model for predicting rock damage and fragment size distribution due to explosive loading. The model is based on continuum mechanics and statistical fracture mechanics, assuming the rock medium is an isotropic, continuous and homogeneous material with pre-existing microcracks. In the model, damage to the rock medium is defined as the probability of fracture at a given crack density which is obtained by integrating a crack density function over time. The material constants used in the crack density function are determined according to their physical meaning. The minimum damage value at which the fragments may be formed is set by assuming that there is at least one crack per unit volume. Fragment size distribution is achieved considering the equilibrium between kinetic energy and surface energy. The simulation results are in good accordance with the theory of explosive energy partitioning in a rock medium. As a result, the damage zone induced by the shock wave and stress waves, once established, remains stable. The model has been calibrated by field crater blasting and small scale bench blasting tests.

  19. Computational method for thermoviscoelasticity with application to rock mechanics

    NASA Astrophysics Data System (ADS)

    Lee, S. C.

    1984-01-01

    Large scale numerical computations associated with rock mechanics problems have required efficient and economical models for predicting temperature, stress, failure, and deformed structural configuration under various loafing conditions. To meet this requirement, the complex dependence of the properties of geological materials on the time and temperature is modified to yield a reduced time scale as a function of time and temperature under the thermorheologically simple material (TSM) postulate. The thermorheologically linear concept is adopted in the finite element formulation by uncoupling thermal and mechanical responses. The thermal responses, based on transient heat conduction or convective diffusion, are formulated by using the two point recurrence scheme and the upwinding scheme, respectively. An incremental solution procedure with the implicit time stepping scheme is proposed for the solution of the thermoviscoelastic response. The proposed thermoviscoelastic solution algorithm is based on the uniaxial creep experimental data and the corresponding temperature shift functions, and is intended to minimize computational efforts by allowing large time step size with stable solutions. A thermoelastic fracture formulation is also presented by introducing the degenerate quadratic isoparametric singular element for the thermally induced line crack problems.

  20. Statistical models of lunar rocks and regolith

    NASA Technical Reports Server (NTRS)

    Marcus, A. H.

    1973-01-01

    The mathematical, statistical, and computational approaches used in the investigation of the interrelationship of lunar fragmental material, regolith, lunar rocks, and lunar craters are described. The first two phases of the work explored the sensitivity of the production model of fragmental material to mathematical assumptions, and then completed earlier studies on the survival of lunar surface rocks with respect to competing processes. The third phase combined earlier work into a detailed statistical analysis and probabilistic model of regolith formation by lithologically distinct layers, interpreted as modified crater ejecta blankets. The fourth phase of the work dealt with problems encountered in combining the results of the entire project into a comprehensive, multipurpose computer simulation model for the craters and regolith. Highlights of each phase of research are given.

  1. A mechanism for high wall-rock velocities in rockbursts

    USGS Publications Warehouse

    McGarr, A.

    1997-01-01

    Considerable evidence has been reported for wall-rock velocities during rockbursts in deep gold mines that are substantially greater than ground velocities associated with the primary seismic events. Whereas varied evidence suggests that slip across a fault at the source of an event generates nearby particle velocities of, at most, several m/s, numerous observations, in nearby damaged tunnels, for instance, imply wall-rock velocities of the order of 10 m/s and greater. The common observation of slab buckling or breakouts in the sidewalls of damaged excavations suggests that slab flexure may be the mechanism for causing high rock ejection velocities. Following its formation, a sidewall slab buckles, causing the flexure to increase until the stress generated by flexure reaches the limit 5 that can be supported by the sidewall rock. I assume here that S is the uniaxial compressive strength. Once the flexural stress exceeds S, presumably due to the additional load imposed by a nearby seismic event, the slab fractures and unflexes violently. The peak wall-rock velocity v thereby generated is given by v=(3 + 1-??2/2)1 2 S/?????E for rock of density ??, Young's modulus E, and Poisson's ratio ??. Typical values of these rock properties for the deep gold mines of South Africa yield v= 26 m/s and for especially strong quartzites encountered in these same mines, v> 50m/s. Even though this slab buckling process leads to remarkably high ejection velocities and violent damage in excavations, the energy released during this failure is only a tiny fraction of that released in the primary seismic event, typically of magnitude 2 or greater.

  2. DEM analyses of shear behaviour of rock joints by a novel bond contact model

    NASA Astrophysics Data System (ADS)

    Jiang, M. J.; Liu, J.; Sun, C.; Chen, H.

    2015-09-01

    The failure of rock joints is one of the potential causes for the local and general rock instability, which may trigger devastating geohazards such as landslide. In this paper, the Distinct Element Method (DEM) featured by a novel bond contact model was utilized to simulate shear behaviour of centre/non-coplanar rock joints. The DEM results show that the complete shear behaviour of jointed rock includes four stages: elastic shearing phase, crack propagation, the failure of rock bridges and the through-going discontinuity. The peak shear strength of centre joint increases as the joint connectivity rate decreases. For intermittent noncoplanar rock joints, as the inclination of the rock joints increases, its shear capacity decreases when the inclination angle is negative while increase when positive. Comparison with the experimental results proves the capability of this DEM model in capturing the mechanical properties of the jointed rocks.

  3. Assessment of rock wool as support material for on-site sanitation: hydrodynamic and mechanical characterization.

    PubMed

    Wanko, Adrien; Laurent, Julien; Bois, Paul; Mosé, Robert; Wagner-Kocher, Christiane; Bahlouli, Nadia; Tiffay, Serge; Braun, Bouke; Provo kluit, Pieter-Willem

    2016-01-01

    This study proposes mechanical and hydrodynamic characterization of rock wool used as support material in compact filter. A double-pronged approach, based on experimental simulation of various physical states of this material was done. First of all a scanning electron microscopy observation allows to highlight the fibrous network structure, the fibres sizing distribution and the atomic absorption spectrum. The material was essentially lacunar with 97 ± 2% of void space. Static compression tests on variably saturated rock wool samples provide the fact that the strain/stress behaviours depend on both the sample conditioning and the saturation level. Results showed that water exerts plastifying effect on mechanical behaviour of rock wool. The load-displacement curves and drainage evolution under different water saturation levels allowed exhibiting hydraulic retention capacities under stress. Finally, several tracer experiments on rock wool column considering continuous and batch feeding flow regime allowed: (i) to determine the flow model for each test case and the implications for water dynamic in rock wool medium, (ii) to assess the rock wool double porosity and discuss its advantages for wastewater treatment, (iii) to analyse the benefits effect for water treatment when the high level of rock wool hydric retention was associated with the plug-flow effect, and (iv) to discuss the practical contributions for compact filter conception and management. PMID:26165374

  4. Deformation mechanisms in a coal mine roadway in extremely swelling soft rock.

    PubMed

    Li, Qinghai; Shi, Weiping; Yang, Renshu

    2016-01-01

    The problem of roadway support in swelling soft rock was one of the challenging problems during mining. For most geological conditions, combinations of two or more supporting approaches could meet the requirements of most roadways; however, in extremely swelling soft rock, combined approaches even could not control large deformations. The purpose of this work was to probe the roadway deformation mechanisms in extremely swelling soft rock. Based on the main return air-way in a coal mine, deformation monitoring and geomechanical analysis were conducted, as well as plastic zone mechanical model was analysed. Results indicated that this soft rock was potentially very swelling. When the ground stress acted alone, the support strength needed in situ was not too large and combined supporting approaches could meet this requirement; however, when this potential released, the roadway would undergo permanent deformation. When the loose zone reached 3 m within surrounding rock, remote stress p ∞ and supporting stress P presented a linear relationship. Namely, the greater the swelling stress, the more difficult it would be in roadway supporting. So in this extremely swelling soft rock, a better way to control roadway deformation was to control the releasing of surrounding rock's swelling potential. PMID:27547684

  5. Centrifuge modeling of rocking-isolated inelastic RC bridge piers

    PubMed Central

    Loli, Marianna; Knappett, Jonathan A; Brown, Michael J; Anastasopoulos, Ioannis; Gazetas, George

    2014-01-01

    Experimental proof is provided of an unconventional seismic design concept, which is based on deliberately underdesigning shallow foundations to promote intense rocking oscillations and thereby to dramatically improve the seismic resilience of structures. Termed rocking isolation, this new seismic design philosophy is investigated through a series of dynamic centrifuge experiments on properly scaled models of a modern reinforced concrete (RC) bridge pier. The experimental method reproduces the nonlinear and inelastic response of both the soil-footing interface and the structure. To this end, a novel scale model RC (1:50 scale) that simulates reasonably well the elastic response and the failure of prototype RC elements is utilized, along with realistic representation of the soil behavior in a geotechnical centrifuge. A variety of seismic ground motions are considered as excitations. They result in consistent demonstrably beneficial performance of the rocking-isolated pier in comparison with the one designed conventionally. Seismic demand is reduced in terms of both inertial load and deck drift. Furthermore, foundation uplifting has a self-centering potential, whereas soil yielding is shown to provide a particularly effective energy dissipation mechanism, exhibiting significant resistance to cumulative damage. Thanks to such mechanisms, the rocking pier survived, with no signs of structural distress, a deleterious sequence of seismic motions that caused collapse of the conventionally designed pier. © 2014 The Authors Earthquake Engineering & Structural Dynamics Published by John Wiley & Sons Ltd. PMID:26300573

  6. Mechanical rock properties, fracture propagation and permeability development in deep geothermal reservoirs

    NASA Astrophysics Data System (ADS)

    Leonie Philipp, Sonja; Reyer, Dorothea

    2010-05-01

    layers. Thus, to minimise exploration risks and for effective stimulation, the geometry of the fracture system and the mechanical properties of the host rock must be known. Here we present first results of structural geological field studies of fracture systems in outcrop analogues studies of rocks that could be used to host man-made geothermal reservoirs in sedimentary rocks in the North German Basin. As examples, we show data from different lithologies, including Buntsandstein (Lower Triassic), a sandstone-shale succession and Muschelkalk (Middle Triassic), a limestone-marl succession. We analyse natural fracture systems and the effects of rock heterogeneities, particularly stiffness variations between layers (mechanical layering) on the propagation of natural fractures. Important fracture parameters include attitude, aperture and interconnectivity to fracture systems. The field studies are supplemented by laboratory measurements of the above mentioned rock mechanical properties. Our field studies indicate that many fractures become arrested at layer contacts, particularly at contacts between layers with contrasting mechanical properties. Measurements of thousands of fractures indicate that even very thin layers (mm to cm-scale thicknesses) of shale or marl may be responsible for the arrest of many fractures. Our results suggest that the propagation and aperture variation of fractures are important parameters in the permeability development of deep geothermal reservoirs. These studies provide a basis for models of fracture networks and fluid transport in future man-made reservoirs. We conclude that the likely permeability of a man-made geothermal reservoir can be inferred from field data from outcrop analogues, laboratory measurements, and numerical models.

  7. Characterizing and modelling 'ghost-rock' weathered limestones

    NASA Astrophysics Data System (ADS)

    Dubois, Caroline; Goderniaux, Pascal; Deceuster, John; Poulain, Angélique; Kaufmann, Olivier

    2016-04-01

    'Ghost-rock' karst aquifer has recently been highlighted. In this particular type of aquifer, the karst is not expressed as open conduits but consists in zones where the limestone is weathered. The in-situ weathering of limestone leaves a soft porous material called 'alterite'. The hydro-mechanical properties of this material differs significantly from those of the host rock: the weathering enhances the storage capacity and the conductivity of the rock. This type of weathered karst aquifer has never been studied from a hydrogeological point of view. In this study, we present the hydraulic characterization of such weathered zones. We also present a modelling approach derived from the common Equivalent Porous Medium (EPM) approach, but including the spatial distribution of hydrogeological properties through the weathered features, from the hard rock to the alterite, according to a weathering index. Unlike the Discrete Fracture Network (DFN) approaches, which enable to take into account a limited number of fractures, this new approach allows creating models including thousands of weathered features. As the properties of the alterite have to be considered at a centimeter scale, it is necessary to upscale these properties to carry out simulations over large areas. Therefore, an upscaling method was developed, taking into account the anisotropy of the weathered features. Synthetic models are built, upscaled and different hydrogeological simulations are run to validate the method. This methodology is finally tested on a real case study: the modelling of the dewatering drainage flow of an exploited quarry in a weathered karst aquifer in Belgium.

  8. Geometric Effect of Asperities on Shear Mechanism of Rock Joints

    NASA Astrophysics Data System (ADS)

    Fathi, Ali; Moradian, Zabihallah; Rivard, Patrice; Ballivy, Gérard; Boyd, Andrew J.

    2016-03-01

    Three-dimensional tracking of changes of asperities is one of the most important ways to illustrate shear mechanism of rock joints during testing. In this paper, the changes of the role of asperities during different stages of shearing are described by using a new methodology for the characterization of the asperities. The basis of the proposed method is the examination of the three-dimensional roughness of joint surfaces scanned before and after shear testing. By defining a concept named `tiny window', the geometric model of the joint surfaces is reconstructed. Tiny windows are expressed as a function of the x and y coordinates, the height (z coordinate), and the angle of a small area of the surface. Constant normal load (CNL) direct shear tests were conducted on replica joints and, by using the proposed method, the distribution and size of contact and damaged areas were identified. Image analysis of the surfaces was used to verify the results of the proposed method. The results indicated that the proposed method is suitable for determining the size and distribution of the contact and damaged areas at any shearing stage. The geometric properties of the tiny windows in the pre-peak, peak, post-peak softening, and residual shearing stages were investigated based on their angle and height. It was found that tiny windows that face the shear direction, especially the steepest ones, have a primary role in shearing. However, due to degradation of asperities at higher normal stresses and shear displacements, some of the tiny windows that do not initially face the shear direction also come in contact. It was also observed that tiny windows with different heights participate in the shearing process, not just the highest ones. Total contact area of the joint surfaces was considered as summation of just-in-contact areas and damaged areas. The results of the proposed method indicated that considering differences between just-in-contact areas and damaged areas provide useful

  9. Panel discussion on rock mechanics issues in repository design

    SciTech Connect

    Bieniawski, Z.T.; Kim, K.S.; Nataraja, M.

    1996-04-01

    The panel discussion was introduced by Mr. Z.T.(Richard) Bieniawski and then continued with five additional speakers. The topics covered in the discussion included rock mechanics pertaining to the design of underground facilities for the disposal of radioactive wastes and the safety of such facilities. The speakers included: Mr. Kun-Soo Kim who is a specialist in the area of rock mechanics testing during the Basalt Waste Isolation Project; Dr. Mysore Nataraja who is the senior project manager with the NRC; Dr. Michael Voegele who is the project manager for Science Applications International Corporation (SAIC) on the Yucca Mountain Project; Dr. Edward Cording who is a member of the Nuclear Waste Technical Review Board; and Dr. Hemendra Kalia who is employed by Los Alamos National Laboratory and coordinates various activities of testing programs at the Yucca Mountain Site.

  10. Subcritical crack growth and mechanical weathering: a new consideration of how moisture influences rock erosion rates.

    NASA Astrophysics Data System (ADS)

    Eppes, Martha-Cary; Keanini, Russell; Hancock, Gregory S.

    2016-04-01

    The contributions of moisture to the mechanical aspects of rock weathering and regolith production are poorly quantified. In particular, geomorphologists have largely overlooked the role of subcritical crack growth processes in physical weathering and the fact that moisture strongly influences the rates of those processes. This influence is irrespective of the function that moisture plays in stress loading mechanisms like freezing or hydration. Here we present a simple numerical model that explores the efficacy of subcritical crack growth in granite rock subaerially exposed under a range of moisture conditions. Because most weathering-related stress loading for rocks found at, or near, Earth's surface (hereafter surface rocks) is cyclic, we modeled crack growth using a novel combination of Paris' Law and Charles' Law. This combination allowed us to apply existing empirically-derived data for the stress corrosion index of Charles' Law to fatigue cracking. For stress, we focused on the relatively straightforward case of intergranular stresses that arise during solar-induced thermal cycling by conductive heat transfer, making the assumption that such stresses represent a universal minimum weathering stress experienced by all surface rocks. Because all other tensile weathering-related stresses would be additive in the context of crack growth, however, our model can be adapted to include other stress loading mechanisms. We validated our calculations using recently published thermal-stress-induced cracking rates. Our results demonstrate that 1) weathering-induced stresses as modeled herein, and as published by others, are sufficient to propagate fractures subcritically over long timescales with or without the presence of water 2) fracture propagation rates increase exponentially with respect to moisture, specifically relative humidity 3) fracture propagation rates driven by thermal cycling are strongly dependent on the magnitude of diurnal temperature ranges and the

  11. Modelling study on uranium migration in rocks under weathering condition

    SciTech Connect

    Ohnuki, Toshihiko; Isobe, Hiroshi; Sato, Tsutomu; Yanase, Nobuyuki; Murakami, Takashi

    1995-12-31

    A modelling study has been completed to understand the effect of rock alteration on uranium migration at the Koongarra ore deposit, Australia. The model considers the weathering process, the mechanism and rate of chlorite alteration, a major mineral of the host rock, and assumes the presence of reversible sorption sites of chlorite and the presence of reversible and irreversible sorption sites of the weathering products. One- and two-dimensional, calculated uranium concentrations were compared with those observed. Good agreement between the calculated and observed uranium concentration profiles was obtained only when an appropriate fraction of uranium is fixed to the irreversible sorption sites of Fe-minerals produced during weathering of chlorite. On the other hand, the conventional Kd model failed to estimate an adequate uranium concentration profile. The results suggest that the fixation of uranium to Fe-minerals has dominated the migration of uranium in the vicinity of the Koongarra ore deposit.

  12. Stress Evolution in Roadway Rock Bolts During Mining in a Fully Mechanized Longwall Face, and an Evaluation of Rock Bolt Support Design

    NASA Astrophysics Data System (ADS)

    Zhang, Kai; Zhang, Guimin; Hou, Rongbin; Wu, Yu; Zhou, Hongqi

    2015-01-01

    Rock bolts are widely used in coal mines throughout China. Approximately 8,000 km of roadways are excavated in coal mines every year in China, 80 % of which are supported by rock bolts. At present, the design of rock bolt support schemes is mainly based on analogies and experience from previous projects. In the present study, in order to evaluate the design of rock bolt support in roadways, several cross sections of a roadway were monitored for rock bolt stress during the roadway excavation and mining. The study results show that the stress in the rock bolts varied in the areas 20 m behind the excavating face and 30 m ahead of the mining face. For the rock bolts observed in this study, the max axial force was within the design limit of the bolts, thus the support design was shown to be acceptable. Then, numerical simulation was performed using FLAC3D to investigate the stress evolution in the rock bolts during the mining of the fully mechanized longwall face. The simulation results show an overall agreement with the in situ measurements. Finally, parametric study pertaining to length, anchorage length, and rock bolt spacing was carried out with the numerical model, and several suggestions for the support design were proposed.

  13. Textural remanence - A new model of lunar rock magnetism

    NASA Technical Reports Server (NTRS)

    Brecher, A.

    1976-01-01

    In reexamining the accumulated magnetic data on lunar rocks, several common patterns of magnetic behavior are recognized. Their joint occurrence strongly suggests a new model of lunar rock magnetism, which is based on partial preferred textural alignment of the spontaneous moments of magnetic grains without requiring the existence of ancient lunar magnetic fields. This magnetic fabric, mimetic to locally oriented petrofabric, gives rise to an apparent 'textural remanent magnetization'. In order to account for the observed intensity of 'stable remanence' in lunar rocks, only a minute fraction (0.001 to 0.00001) of the single-domain iron grains present need be preferentially aligned. Several mechanisms operating on the lunar surface, including shock and diurnal thermal cycling, appear adequate for producing the required type and degree of magnetic alignment in all lunar rock classes. The model is supported by a wide variety of direct and indirect evidence, and its predictions (e.g., regarding anisotropic susceptibility and remanence acquisition) can be experimentally tested.

  14. Failure Behavior and Constitutive Model of Weakly Consolidated Soft Rock

    PubMed Central

    Wang, Wei-ming; Zhao, Zeng-hui; Wang, Yong-ji; Gao, Xin

    2013-01-01

    Mining areas in western China are mainly located in soft rock strata with poor bearing capacity. In order to make the deformation failure mechanism and strength behavior of weakly consolidated soft mudstone and coal rock hosted in Ili No. 4 mine of Xinjiang area clear, some uniaxial and triaxial compression tests were carried out according to the samples of rocks gathered in the studied area, respectively. Meanwhile, a damage constitutive model which considered the initial damage was established by introducing a damage variable and a correction coefficient. A linearization process method was introduced according to the characteristics of the fitting curve and experimental data. The results showed that samples under different moisture contents and confining pressures presented completely different failure mechanism. The given model could accurately describe the elastic and plastic yield characteristics as well as the strain softening behavior of collected samples at postpeak stage. Moreover, the model could precisely reflect the relationship between the elastic modulus and confining pressure at prepeak stage. PMID:24489511

  15. Rock mechanics. Superplastic nanofibrous slip zones control seismogenic fault friction.

    PubMed

    Verberne, Berend A; Plümper, Oliver; de Winter, D A Matthijs; Spiers, Christopher J

    2014-12-12

    Understanding the internal mechanisms controlling fault friction is crucial for understanding seismogenic slip on active faults. Displacement in such fault zones is frequently localized on highly reflective (mirrorlike) slip surfaces, coated with thin films of nanogranular fault rock. We show that mirror-slip surfaces developed in experimentally simulated calcite faults consist of aligned nanogranular chains or fibers that are ductile at room conditions. These microstructures and associated frictional data suggest a fault-slip mechanism resembling classical Ashby-Verrall superplasticity, capable of producing unstable fault slip. Diffusive mass transfer in nanocrystalline calcite gouge is shown to be fast enough for this mechanism to control seismogenesis in limestone terrains. With nanogranular fault surfaces becoming increasingly recognized in crustal faults, the proposed mechanism may be generally relevant to crustal seismogenesis. PMID:25504714

  16. Modelling of reactive fluid transport in deformable porous rocks

    NASA Astrophysics Data System (ADS)

    Yarushina, V. M.; Podladchikov, Y. Y.

    2009-04-01

    One outstanding challenge in geology today is the formulation of an understanding of the interaction between rocks and fluids. Advances in such knowledge are important for a broad range of geologic settings including partial melting and subsequent migration and emplacement of a melt into upper levels of the crust, or fluid flow during regional metamorphism and metasomatism. Rock-fluid interaction involves heat and mass transfer, deformation, hydrodynamic flow, and chemical reactions, thereby necessitating its consideration as a complex process coupling several simultaneous mechanisms. Deformation, chemical reactions, and fluid flow are coupled processes. Each affects the others. Special effort is required for accurate modelling of the porosity field through time. Mechanical compaction of porous rocks is usually treated under isothermal or isoentropic simplifying assumptions. However, joint consideration of both mechanical compaction and reactive porosity alteration requires somewhat greater than usual care about thermodynamic consistency. Here we consider the modelling of multi-component, multi-phase systems, which is fundamental to the study of fluid-rock interaction. Based on the conservation laws for mass, momentum, and energy in the form adopted in the theory of mixtures, we derive a thermodynamically admissible closed system of equations describing the coupling of heat and mass transfer, chemical reactions, and fluid flow in a deformable solid matrix. Geological environments where reactive transport is important are located at different depths and accordingly have different rheologies. In the near surface, elastic or elastoplastic properties would dominate, whereas viscoplasticity would have a profound effect deeper in the lithosphere. Poorly understood rheologies of heterogeneous porous rocks are derived from well understood processes (i.e., elasticity, viscosity, plastic flow, fracturing, and their combinations) on the microscale by considering a

  17. A probabilistic approach to rock mechanical property characterization for nuclear waste repository design

    SciTech Connect

    Kim, Kunsoo; Gao, Hang

    1996-04-01

    A probabilistic approach is proposed for the characterization of host rock mechanical properties at the Yucca Mountain site. This approach helps define the probability distribution of rock properties by utilizing extreme value statistics and Monte Carlo simulation. We analyze mechanical property data of tuff obtained by the NNWSI Project to assess the utility of the methodology. The analysis indicates that laboratory measured strength and deformation data of Calico Hills and Bullfrog tuffs follow an extremal. probability distribution (the third type asymptotic distribution of the smallest values). Monte Carlo simulation is carried out to estimate rock mass deformation moduli using a one-dimensional tuff model proposed by Zimmermann and Finley. We suggest that the results of these analyses be incorporated into the repository design.

  18. Quantifying the impact of lithology upon the mechanical properties of rock

    NASA Astrophysics Data System (ADS)

    Weatherley, Dion

    2013-04-01

    The physical characteristics of rock, its lithology, undoubtedly influences its deformation under natural or engineering loads. Mineral texture, micro-damage, joints, bedding planes, inclusions, unconformities and faults are all postulated to alter the mechanical response of rock on different scales and under different stressing conditions. Whilst laboratory studies have elucidated some aspects of the relationship between lithology and mechanical properties, these small-scale results are difficult to extrapolate to lithospheric scales. To augment laboratory-derived knowledge, physics-based numerical modelling is a promising avenue [3]. Bonded particle models implemented using the Discrete Element Method (DEM [1]) are a practical numerical laboratory to investigate the interplay between lithology and the mechanical response of rock specimens [4]. Numerical rock specimens are represented as an assembly of indivisible spherical particles connected to nearest neighbours via brittle-elastic beams which impart forces and moments upon one-another as particles move relative to each other. By applying boundary forces and solving Newton's Laws for each particle, elastic deformation and brittle failure may be simulated [2]. Each beam interaction is defined by four model parameters: Young's modulus, Poisson's ratio, cohesive strength and internal friction angle. Beam interactions in different subvolumes of the specimen are assigned different parameters to model different rock types or mineral assemblages. Micro-cracks, joints, unconformities and faults are geometrically incorporated by fitting particles to either side of triangulated surfaces [5]. The utility of this modelling approach is verified by reproducing analytical results from fracture mechanics (Griffith crack propagation and wing-crack formation) and results of controlled laboratory investigations. To quantify the impact of particular lithologic structures on mechanical response, a range of control experiments are

  19. The interplay of predefined rock mechanics and permafrost forcing in a steep alpine rock crest (Steintaelli, Mattertal, Switzerland)

    NASA Astrophysics Data System (ADS)

    Halla, Christian; Krautblatter, Michael; Draebing, Daniel

    2014-05-01

    Freeze-thaw-processes in the active layer and degrading permafrost change ice and hydraulic pressures as well as rock- and ice-mechanical properties in rock masses which can cause instabilities. The characterization of the rock mass was determined by the geological strength index and a detailed discontinuity analysis along scanlines where the active layer reached depths of 5-15m (August 2012). Rock mass deformations and accordingly the divergence and convergence of deep reaching fractures were measured along 18 extensometer transects with various lengths from 2 to 27m. (1) The Geological Strength Index provided rock mechanical parameters which indicate stable conditions on the slope scale. The friction angle of the rock mass (44°) is higher than the mean slope inclination (37°). (2) The discontinuity analysis provided six joint sets, their geometries and mechanical properties. High roughness coefficients and wall strengths of the joints result in high total friction angles (>43°) and stable conditions on the block scale. However, the locations of several joints with wide apertures, ice fillings and joints influenced by snowmelt of the cornice at the crest are used as indicators for freeze-thaw related rock mass deformations. (3) The rock mass deformation rates during late summer (2012) are multiple times (>3) higher than deformation rates for several years (2008-2012). Furthermore, the direction of deformation changes between divergence and convergence over time. Both aspects indicate that seasonal and annual changes affect the rock mass deformations. We conclude, that rock mass deformation cannot be explained solely by rock mechanical parameters without freeze-thaw and permafrost influence. The higher deformation rates during the thaw season in late summer and the heterogeneity of deformation directions indicate that seasonal and annual changes of ice- and hydraulic pressures in discontinuities affect instabilities at the Steintaelli crestline. Here we show the

  20. Discrete fracture hydromechanical model for the disturbed rock zone in a clay rock

    NASA Astrophysics Data System (ADS)

    Asahina, D.; Houseworth, J. E.; Birkholzer, J. T.

    2013-12-01

    We have developed a coupled thermal-hydrological-mechanical (THM) fracture damage model, TOUGH-RBSN, to investigate the behavior of fracture generation and evolution in rock in the presence of perturbations to THM conditions. This model combines the capabilities of the TOUGH2 simulator to represent thermal-hydrological processes with a rigid-body-spring-network (RBSN) model, a type of discrete modeling, to treat geomechanical and fracture-damage processes. In particular, the development and evolution of fractures in the excavation damaged zone (EDZ) of a clay rock, with application to high-level nuclear waste disposal, is a focus for this model development. Previously, the TOUGH-RBSN approach has been used to model fracture damage under tensile conditions as a result of desiccation shrinkage. The next phase of model testing will be application to the HG-A test being conducted at the Mont Terri underground research laboratory (URL) near Saint-Ursanne, Switzerland. This test is being conducted in a 13-m long, 1-m diameter microtunnel in the Opalinus clay rock in which a test section at the far end of the microtunnel is isolated using a packer. The test is specifically targeted to observe how fluids injected into the test section penetrate into the rock, with particular emphasis on the EDZ. The HG-A microtunnel was excavated in 2005 and subsequent mapping of the tunnel surface shows preferential fracturing and tunnel breakouts along zones where bedding planes are tangential to the tunnel wall and where faults intercept the tunnel. It appears that the EDZ fracture damage can be attributed to both tensile and shear fracturing mechanisms. A series of injection tests with water and gas have been performed which also show preferential invasion of the fluid pressure along the observed damage zones, as well as fracture self-sealing over time. The TOUGH-RBSN approach has been successfully applied to modeling fracture driven by predominately tensile loading, whereas only

  1. Rocks.

    ERIC Educational Resources Information Center

    Lee, Alice

    This science unit is designed for limited- and non-English speaking students in a Chinese bilingual education program. The unit covers rock material, classification, characteristics of types of rocks, and rock cycles. It is written in Chinese and simple English. At the end of the unit there is a list of main terms in both English and Chinese, and…

  2. Effects of bioleaching on the mechanical and chemical properties of waste rocks

    NASA Astrophysics Data System (ADS)

    Yin, Sheng-Hua; Wu, Ai-Xiang; Wang, Shao-Yong; Ai, Chun-Ming

    2012-01-01

    Bioleaching processes cause dramatic changes in the mechanical and chemical properties of waste rocks, and play an important role in metal recovery and dump stability. This study focused on the characteristics of waste rocks subjected to bioleaching. A series of experiments were conducted to investigate the evolution of rock properties during the bioleaching process. Mechanical behaviors of the leached waste rocks, such as failure patterns, normal stress, shear strength, and cohesion were determined through mechanical tests. The results of SEM imaging show considerable differences in the surface morphology of leached rocks located at different parts of the dump. The mineralogical content of the leached rocks reflects the extent of dissolution and precipitation during bioleaching. The dump porosity and rock size change under the effect of dissolution, precipitation, and clay transportation. The particle size of the leached rocks decreased due to the loss of rock integrity and the conversion of dry precipitation into fine particles.

  3. Hitherto unknown shear rupture mechanism as a source of instability in intact hard rocks at highly confined compression

    NASA Astrophysics Data System (ADS)

    Tarasov, Boris G.

    2014-05-01

    Today, frictional shear resistance along pre-existing faults is considered to be the lower limit on rock shear strength for confined conditions corresponding to the seismogenic layer. This paper introduces a recently identified shear rupture mechanism providing a paradoxical feature of hard rocks - the possibility of shear rupture propagation through the highly confined intact rock mass at shear stress levels significantly less than frictional strength. In the new mechanism, the rock failure associated with consecutive creation of small slabs (known as ‘domino-blocks') from the intact rock in the rupture tip is driven by a fan-shaped domino structure representing the rupture head. The fan-head combines such unique features as: extremely low shear resistance, self-sustaining stress intensification, and self-unbalancing conditions. Due to this the failure process caused by the mechanism is very dynamic and violent. This makes it impossible to directly observe and study the mechanism and can explain why the mechanism has not been detected before. This paper provides physical motivation for the mechanism, based upon side effects accompanying the failure process. Physical and mathematical models of the mechanism presented in the paper explain unique and paradoxical features of the mechanism. The new shear rupture mechanism allows a novel point of view for understanding the nature of spontaneous failure processes in hard rocks including earthquakes.

  4. Rock shape, restitution coefficients and rockfall trajectory modelling

    NASA Astrophysics Data System (ADS)

    Glover, James; Christen, Marc; Bühler, Yves; Bartelt, Perry

    2014-05-01

    Restitution coefficients are used in rockfall trajectory modelling to describe the ratio between incident and rebound velocities during ground impact. They are central to the problem of rockfall hazard analysis as they link rock mass characteristics to terrain properties. Using laboratory experiments as a guide, we first show that restitution coefficients exhibit a wide range of scatter, although the material properties of the rock and ground are constant. This leads us to the conclusion that restitution coefficients are poor descriptors of rock-ground interaction. The primary problem is that "apparent" restitution coefficients are applied at the rock's centre-of-mass and do not account for rock shape. An accurate description of the rock-ground interaction requires the contact forces to be applied at the rock surface with consideration of the momentary rock position and spin. This leads to a variety of rock motions including bouncing, sliding, skipping and rolling. Depending on the impact configuration a wide range of motions is possible. This explains the large scatter of apparent restitution coefficients. We present a rockfall model based on newly developed hard-contact algorithms which includes the effects of rock shape and therefore is able to reproduce the results of different impact configurations. We simulate the laboratory experiments to show that it is possible to reproduce run-out and dispersion of different rock shapes using parameters obtained from independent tests. Although this is a step forward in rockfall trajectory modelling, the problem of parametersing real terrain remains.

  5. Deep-seated slowly moving rock slides in foliated metamorphic rock masses: New findings about kinematical and hydro-mechanical processes

    NASA Astrophysics Data System (ADS)

    Zangerl, Christian; Strauhal, Thomas; Holzmann, Michael

    2013-04-01

    Deep-seated slowly moving rock slides are characterised by deformation along one or several shear zones where most of the measured total slope displacement localizes. Many of these rock slides move downwards at mean annual rates of some centimetres or even less and do not show any evidence for non-reversible acceleration in the past or in the future. Whereas some of these rock slides are currently inactive (dormant) or have even reached a stabilised final state others show a temporally variable deformation behaviour characterized by low base activities superimposed by acceleration phases. The trigger for these phases can be manifold and include heavy rainfall, snow melt, water level fluctuations of reservoirs, changes in the slope's equilibrium state due to antecedent slow creeping processes, variations in the material behaviour within the shear/sliding zone, erosion along the foot of the slope, etc. In order to improve the understanding of the activity behaviour and trigger factors, to increase the quality of slope stability analyses and to assess the hazard potential detailed information about the rock slide geometry and kinematics are essential. Given that subsurface investigations such as boreholes and investigation adits on large-scale rock slides are costly, most published studies are related to investigations in the surroundings of infrastructures and human settlements. Within this study new field mapping, deformation monitoring, geophysical exploration and in-situ subsurface investigation data are presented which are obtained on case studies in paragneissic rock masses of the Austroalpine Ötztal-Stubai complex (Tyrol, Austria). The new investigations contribute to develop geometrical rock slide models, to study the internal deformation characteristics of the rock slide mass and to develop kinematical deformation models. In addition, results show that all case studies are characterised by slope deformation mechanisms due to shear/slide processes along dm to

  6. Mechanical effects associated with surface loading of dry rock due to glaciation

    SciTech Connect

    Wahi, K.K.; Hunter, R.L.

    1985-01-01

    Many scenarios of interest for a repository in the Pasco Basin begin with glaciation. Loading and unloading of joints and fractures due to the weight of ice sheets could affect the hydrologic properties of the host rock and surrounding units. Scoping calculations performed using two-dimensional numerical models with simplifying assumptions predict stress changes and uplift or subsidence caused by an advancing glacier. The magnitudes of surface uplift and subsidence predicted by the study agree well with previous independent predictions. Peak stress unloading near the repository horizon is a small fraction of the ambient stress. Any resultant aperture increase is likewise small. Based on the results of this study, mechanical loading caused by a glacier is expected to have a minimal effect on rock permeability, assuming that the excess compressive loads do not crush the rock. 13 refs., 3 figs., 1 tab.

  7. Modeling transport kinetics in clinoptilolite-phosphate rock systems

    NASA Technical Reports Server (NTRS)

    Allen, E. R.; Ming, D. W.; Hossner, L. R.; Henninger, D. L.

    1995-01-01

    Nutrient release in clinoptilolite-phosphate rock (Cp-PR) systems occurs through dissolution and cation-exchange reactions. Investigating the kinetics of these reactions expands our understanding of nutrient release processes. Research was conducted to model transport kinetics of nutrient release in Cp-PR systems. The objectives were to identify empirical models that best describe NH4, K, and P release and define diffusion-controlling processes. Materials included a Texas clinoptilolite (Cp) and North Carolina phosphate rock (PR). A continuous-flow thin-disk technique was used. Models evaluated included zero order, first order, second order, parabolic diffusion, simplified Elovich, Elovich, and power function. The power-function, Elovich, and parabolic-diffusion models adequately described NH4, K, and P release. The power-function model was preferred because of its simplicity. Models indicated nutrient release was diffusion controlled. Primary transport processes controlling nutrient release for the time span observed were probably the result of a combination of several interacting transport mechanisms.

  8. Mechanical twinning as stress indicator in fault rocks

    NASA Astrophysics Data System (ADS)

    Wenk, H.

    2011-12-01

    At low stresses and elevated temperatures rocks deform by dislocation movements and diffusion. At very high stresses they undergo brittle failure. For many minerals there is an intermediate regime where mechanical twinning occurs. This has been studied extensively in calcite (Turner, Griggs and Heard, GSA Mem. 1954) and also documented for quartz (Tullis, Science, 1972). In this study we use twinning microstructures to characterize rocks that were subjected to seismic stresses, specifically pseudotachylites and samples from the San Andreas Fault Observatory at Depth. For calcite in SAFOD samples, dislocation densities derived from TEM images as well as twin densities measured by optical microscopy, indicate stresses between 50 and 200 MPa. Similar residual stress magnitudes were obtained from preserved lattice distortion determined by synchrotron X-ray microdiffraction. Also quartz shows characteristic twin microstructures. Orientation maps with SEM-EBSD reveal that quartz associated with pseudotachylite veins is profusely twinned, similar to structures observed in quartz subjected to meteorite impacts. It suggests that local dynamic seismic stresses during earthquakes are responsible for Dauphiné twinning. Thus microstructures in minerals add information to constrain macroscopic conditions during faulting.

  9. Rock thermal conductivity as key parameter for geothermal numerical models

    NASA Astrophysics Data System (ADS)

    Di Sipio, Eloisa; Chiesa, Sergio; Destro, Elisa; Galgaro, Antonio; Giaretta, Aurelio; Gola, Gianluca; Manzella, Adele

    2013-04-01

    The geothermal energy applications are undergoing a rapid development. However, there are still several challenges in the successful exploitation of geothermal energy resources. In particular, a special effort is required to characterize the thermal properties of the ground along with the implementation of efficient thermal energy transfer technologies. This paper focuses on understanding the quantitative contribution that geosciences can receive from the characterization of rock thermal conductivity. The thermal conductivity of materials is one of the main input parameters in geothermal modeling since it directly controls the steady state temperature field. An evaluation of this thermal property is required in several fields, such as Thermo-Hydro-Mechanical multiphysics analysis of frozen soils, designing ground source heat pumps plant, modeling the deep geothermal reservoirs structure, assessing the geothermal potential of subsoil. Aim of this study is to provide original rock thermal conductivity values useful for the evaluation of both low and high enthalpy resources at regional or local scale. To overcome the existing lack of thermal conductivity data of sedimentary, igneous and metamorphic rocks, a series of laboratory measurements has been performed on several samples, collected in outcrop, representative of the main lithologies of the regions included in the VIGOR Project (southern Italy). Thermal properties tests were carried out both in dry and wet conditions, using a C-Therm TCi device, operating following the Modified Transient Plane Source method.Measurements were made at standard laboratory conditions on samples both water saturated and dehydrated with a fan-forced drying oven at 70 ° C for 24 hr, for preserving the mineral assemblage and preventing the change of effective porosity. Subsequently, the samples have been stored in an air-conditioned room while bulk density, solid volume and porosity were detected. The measured thermal conductivity

  10. Fundamental Research on Percussion Drilling: Improved rock mechanics analysis, advanced simulation technology, and full-scale laboratory investigations

    SciTech Connect

    Michael S. Bruno

    2005-12-31

    This report summarizes the research efforts on the DOE supported research project Percussion Drilling (DE-FC26-03NT41999), which is to significantly advance the fundamental understandings of the physical mechanisms involved in combined percussion and rotary drilling, and thereby facilitate more efficient and lower cost drilling and exploration of hard-rock reservoirs. The project has been divided into multiple tasks: literature reviews, analytical and numerical modeling, full scale laboratory testing and model validation, and final report delivery. Literature reviews document the history, pros and cons, and rock failure physics of percussion drilling in oil and gas industries. Based on the current understandings, a conceptual drilling model is proposed for modeling efforts. Both analytical and numerical approaches are deployed to investigate drilling processes such as drillbit penetration with compression, rotation and percussion, rock response with stress propagation, damage accumulation and failure, and debris transportation inside the annulus after disintegrated from rock. For rock mechanics modeling, a dynamic numerical tool has been developed to describe rock damage and failure, including rock crushing by compressive bit load, rock fracturing by both shearing and tensile forces, and rock weakening by repetitive compression-tension loading. Besides multiple failure criteria, the tool also includes a damping algorithm to dissipate oscillation energy and a fatigue/damage algorithm to update rock properties during each impact. From the model, Rate of Penetration (ROP) and rock failure history can be estimated. For cuttings transport in annulus, a 3D numerical particle flowing model has been developed with aid of analytical approaches. The tool can simulate cuttings movement at particle scale under laminar or turbulent fluid flow conditions and evaluate the efficiency of cutting removal. To calibrate the modeling efforts, a series of full-scale fluid hammer

  11. Analysis on the Rock-Cutter Interaction Mechanism During the TBM Tunneling Process

    NASA Astrophysics Data System (ADS)

    Yang, Haiqing; Wang, He; Zhou, Xiaoping

    2016-03-01

    The accurate prediction of rock cutting forces of disc cutters is crucial for tunnel boring machine (TBM) design and construction. Disc cutter wear, which affects TBM penetration performance, has frequently been found at TBM sites. By considering the operating path and wear of the disc cutter, a new model is proposed for evaluating the cutting force and wear of the disc cutter in the tunneling process. The circular path adopted herein, which is the actual running path of the TBM disc cutter, shows that the lateral force of the disc cutter is asymmetric. The lateral forces on the sides of the disc cutter are clearly different. However, traditional solutions are obtained by assuming a linear path, where the later forces are viewed as equal. To simulate the interaction between the rock and disc cutter, a simple brittle damage model for rock mass is introduced here. Based on the explicit dynamic finite element method, the cutting force acting on the rock generated by a single disc cutter is simulated. It is shown that the lateral cutting force of the disc cutter strongly affects the wear extent of disc cutter. The wear mechanism is thus underestimated by the classical model, which was obtained by linear cutting tests. The simulation results are discussed and compared with other models, and these simulation results agree well with the results of present ones.

  12. Modeling temperature and stress in rocks exposed to the sun

    NASA Astrophysics Data System (ADS)

    Hallet, B.; Mackenzie, P.; Shi, J.; Eppes, M. C.

    2012-12-01

    The potential contribution of solar-driven thermal cycling to the progressive breakdown of surface rocks on the Earth and other planets is recognized but under studied. To shed light on this contribution we have launched a collaborative study integrating modern instrumental and numerical approaches to define surface temperatures, stresses, strains, and microfracture activity in exposed boulders, and to shed light on the thermo-mechanical response of boulders to diurnal solar exposure. The instrumental portion of our study is conducted by M. Eppes and coworkers who have monitored the surface and environmental conditions of two ~30 cm dia. granite boulders (one in North Carolina, one in New Mexico) in the field for one and tow years, respectively. Each boulder is instrumented with 8 thermocouples, 8 strain gauges, a surface moisture sensor and 6 acoustic emission (AE) sensors to monitor microfracture activity continuously and to locate it within 2.5 cm. Herein, we focus on the numerical modeling. Using a commercially available finite element program, MSC.Marc®2008r1, we have developed an adaptable, realistic thermo-mechanical model to investigate quantitatively the temporal and spatial distributions of both temperature and stress throughout a boulder. The model accounts for the effects of latitude and season (length of day and the sun's path relative to the object), atmospheric damping (reduction of solar radiation when traveling through the Earth's atmosphere), radiative interaction between the boulder and its surrounding soil, secondary heat exchange of the rock with air, and transient heat conduction in both rock and soil. Using representative thermal and elastic rock properties, as well as realistic representations of the size, shape and orientation of a boulder instrumented in the field in North Carolina, the model is validated by comparison with direct measurements of temperature and strain on the surface of one boulder exposed to the sun. Using the validated

  13. Numerical modelling of fluid-rock interactions: Lessons learnt from carbonate rocks diagenesis studies

    NASA Astrophysics Data System (ADS)

    Nader, Fadi; Bachaud, Pierre; Michel, Anthony

    2015-04-01

    Quantitative assessment of fluid-rock interactions and their impact on carbonate host-rocks has recently become a very attractive research topic within academic and industrial realms. Today, a common operational workflow that aims at predicting the relevant diagenetic processes on the host rocks (i.e. fluid-rock interactions) consists of three main stages: i) constructing a conceptual diagenesis model including inferred preferential fluids pathways; ii) quantifying the resulted diagenetic phases (e.g. depositing cements, dissolved and recrystallized minerals); and iii) numerical modelling of diagenetic processes. Most of the concepts of diagenetic processes operate at the larger, basin-scale, however, the description of the diagenetic phases (products of such processes) and their association with the overall petrophysical evolution of sedimentary rocks remain at reservoir (and even outcrop/ well core) scale. Conceptual models of diagenetic processes are thereafter constructed based on studying surface-exposed rocks and well cores (e.g. petrography, geochemistry, fluid inclusions). We are able to quantify the diagenetic products with various evolving techniques and on varying scales (e.g. point-counting, 2D and 3D image analysis, XRD, micro-CT and pore network models). Geochemical modelling makes use of thermodynamic and kinetic rules as well as data-bases to simulate chemical reactions and fluid-rock interactions. This can be through a 0D model, whereby a certain process is tested (e.g. the likelihood of a certain chemical reaction to operate under specific conditions). Results relate to the fluids and mineral phases involved in the chemical reactions. They could be used as arguments to support or refute proposed outcomes of fluid-rock interactions. Coupling geochemical modelling with transport (reactive transport model; 1D, 2D and 3D) is another possibility, attractive as it provides forward simulations of diagenetic processes and resulting phases. This

  14. An Experimental Study and Constitutive Modeling of Saturated Porous Rocks

    NASA Astrophysics Data System (ADS)

    Xie, S. Y.; Shao, J. F.

    2015-01-01

    This paper is devoted to the experimental characterization and constitutive modeling of saturated porous rocks. A typical porous chalk is investigated. Drained hydrostatic and triaxial compression tests are first performed to characterize the basic mechanical behavior of chalk. Drained triaxial tests with constant interstitial pressure are then carried out to study the effects of interstitial pressure on the plastic deformation and failure criterion. Finally, undrained triaxial compression tests are performed to investigate poromechanical coupling in saturated conditions. Based on the experimental data and some relevant micromechanical considerations, a micromechanics-based plastic model is proposed and extended to poroplastic coupling using the effective stress concept. The proposed model is verified through comparisons between the numerical results and experimental data for both drained and undrained tests.

  15. Simulation of sedimentary rock deformation: Lab-scale model calibration and parameterization

    NASA Astrophysics Data System (ADS)

    Boutt, David F.; McPherson, Brian J. O. L.

    2002-02-01

    Understanding the mechanical behavior of rock is critical for researchers and decision-makers in fields from petroleum recovery to hazardous waste disposal. Traditional continuum-based numerical models are hampered by inadequate constitutive relationships governing fracture initiation and growth. To overcome limits associated with continuum models we employed a discrete model based on the fundamental laws of contact physics to calibrate triaxial tests. Results from simulations of triaxial compression tests on a suite of sedimentary rocks indicate that the basic physics of rock behavior are clearly captured. Evidence for this conclusion lie in the fact that one set of model parameters describes rock behavior at many confining pressures. The use of both inelastic and elastic parameters for comparison yields insight concerning the uniqueness of these models. These tests will facilitate development and calibration of larger scale discrete element models, which may be applied to a wide range of geological problems.

  16. Numerical Modeling of Jointed Rock Under Compressive Loading Using X-ray Computerized Tomography

    NASA Astrophysics Data System (ADS)

    Yu, Qinglei; Yang, Shengqi; Ranjith, P. G.; Zhu, Wancheng; Yang, Tianhong

    2016-03-01

    As jointed rocks consist of joints embedded within intact rock blocks, the presence and geometrical fabric of joints have a great influence on the mechanical behavior of rock. With consideration of the actual spatial shape of joints, a numerical model is proposed to investigate the fracture evolution mechanism of jointed rocks. In the proposed model, computerized tomography (CT) scanning is first used to capture the microstructure of a jointed sandstone specimen, which is artificially fabricated by loading the intact sample until the residual strength, and then digital image processing (DIP) techniques are applied to characterize the geometrical fabric of joints from the CT images. A simple vectorization method is used to convert the microstructure based on a cross-sectional image into a layer of 3-D vectorized microstructure and the overall 3-D model of the jointed sandstone including the real spatial shape of the joints is established by stacking the layers in a specific sequence. The 3-D model is then integrated into a well-established code [three-dimensional Rock Failure Process Analysis, (RFPA3D)]. Using the proposed model, a uniaxial compression test of the jointed sandstone is simulated. The results show that the presence of joints can produce tensile stress zones surrounding them, which result in the fracture of jointed rocks under a relatively small external load. In addition, the spatial shape of the joints has a great influence on the fracture process of jointed rocks.

  17. Petrological modeling of basaltic rocks from Venus: A case for the presence of silicic rocks

    NASA Astrophysics Data System (ADS)

    Shellnutt, J. Gregory

    2013-06-01

    presence of highly evolved igneous rocks on Venus is a controversial issue. The formations of highland terranes and pancake domes are the two principal tectonic and volcanic features which argue in favor of the presence of silicic igneous rocks; however, the lack of water on Venus casts doubt on whether or not granites and rhyolites can form. Data returned to Earth from the Venera 13 and 14 landers show that the surface of Venus is composed of basaltic rocks similar in composition to those found on Earth. Here it is shown that anhydrous and hydrous fractional crystallization modeling using the Venera 13 and 14 data as starting materials can produce compositions similar to terrestrial phonolites and rhyolites. It is suggested that at shallow crustal levels (i.e., ≤ 0.1 GPa), mafic magmas can differentiate into silicic magmas resembling phonolites or rhyolites which may or may not erupt. Furthermore, the hydrous equilibrium partial melting models can produce rocks similar to terrestrial andesites and rhyolites, whereas anhydrous models suggest that there may be a uniquely Venusian type of silicic rock. The silicic rocks, if present, could act as "continental nucleation" sites and/or their presence may facilitate preferential sites of shearing and deformation of the Venusian crust.

  18. Petrological modeling of basaltic rocks from Venus: a case for the presence of silicic rocks

    NASA Astrophysics Data System (ADS)

    Shellnutt, J. G.

    2013-12-01

    The presence of highly evolved igneous rocks on Venus is debated. The formation of highland terranes and pancake domes are the two principle tectonic and volcanic features which argue in favor of the presence of silicic igneous rocks; however, the lack of water on Venus casts doubt on whether or not granites and rhyolites can form. Data returned to Earth from the Venera 13 and 14 landers show that the surface of Venus is comprised of basaltic rocks similar in composition to those found on Earth. Here is it shown that anhydrous and hydrous fractional crystallization modeling using the Venera 13 and 14 data as starting materials can produce compositions similar to terrestrial phonolites and rhyolites. It is suggested that at shallow crustal levels (i.e. ≤ 0.1 GPa) mafic magmas can differentiate into silicic magmas resembling phonolites or rhyolites which may or may not erupt. Furthermore, the hydrous equilibrium partial melting models can produce rocks similar to terrestrial andesites and rhyolites whereas anhydrous models suggest there may be a uniquely Venusian type of silicic rock. The silicic rocks, if present, could act as ';continental nucleation' sites and/or their presence may facilitate preferential sites of shearing and deformation of the Venusian crust.

  19. Mapping the mechanical properties of rocks using automated microindentation tests

    NASA Astrophysics Data System (ADS)

    Masson, Yder; Pride, Steven R.

    2015-10-01

    A microindentation scanner is constructed that measures the spatial fluctuation in the elastic properties of natural rocks. This novel instrument performs automated indentation tests on the surface of a rock slab and outputs 2-D maps of the indentation modulus at submillimeter resolution. Maps obtained for clean, well-consolidated, sandstone are presented and demonstrate the capabilities of the instrument. We observe that the elastic structure of sandstones correlates well with their visual appearance. Further, we show that the probability distribution of the indentation modulus fluctuations across the slab surfaces can be modeled using a lognormal probability density function. To illustrate possible use of the data obtained with the microindentation scanner, we use roughly 10 cm × 10 cm scans with millimeter resolution over four sandstone planar slabs to numerically compute the overall drained elastic moduli for each sandstone sample. We show that such numerically computed moduli are well modeled using the multicomponent form of the Hashin-Shtrikman lower bound that employs the observed lognormal probability distribution for the mesoscopic-scale moduli (the geometric mean works almost the same). We also compute the seismic attenuation versus frequency associated with wave-induced fluid flow between the heterogeneities in the scanned sandstones and observe relatively small values for the inverse quality factor (Q-1<10-2) in the seismic frequency band 102 Hz

  20. Modeling of the nonlinear resonant response in sedimentary rocks

    SciTech Connect

    Ten Cate, James A; Shankland, Thomas J; Vakhnenko, Vyacheslav O; Vakhnenko, Oleksiy

    2009-04-03

    We suggest a model for describing a wide class of nonlinear and hysteretic effects in sedimentary rocks at longitudinal bar resonance. In particular, we explain: hysteretic behaviour of a resonance curve on both its upward and downward slopes; linear softening of resonant frequency with increase of driving level; gradual (almost logarithmic) recovery of resonant frequency after large dynamical strains; and temporal relaxation of response amplitude at fixed frequency. Starting with a suggested model, we predict the dynamical realization of end-point memory in resonating bar experiments with a cyclic frequency protocol. These theoretical findings were confirmed experimentally at Los Alamos National Laboratory. Sedimentary rocks, particularly sandstones, are distinguished by their grain structure in which each grain is much harder than the intergrain cementation material. The peculiarities of grain and pore structures give rise to a variety of remarkable nonlinear mechanical properties demonstrated by rocks, both at quasistatic and alternating dynamic loading. Thus, the hysteresis earlier established for the stress-strain relation in samples subjected to quasistatic loading-unloading cycles has also been discovered for the relation between acceleration amplitude and driving frequency in bar-shaped samples subjected to an alternating external drive that is frequency-swept through resonance. At strong drive levels there is an unusual, almost linear decrease of resonant frequency with strain amplitude, and there are long-term relaxation phenomena such as nearly logarithmic recovery (increase) of resonant frequency after the large conditioning drive has been removed. In this report we present a short sketch of a model for explaining numerous experimental observations seen in forced longitudinal oscillations of sandstone bars. According to our theory a broad set of experimental data can be understood as various aspects of the same internally consistent pattern. Furthermore

  1. A modelling study of the effect of rock alteration on the redistribution of uranium

    SciTech Connect

    Murakami, Takashi; Kimura, Hideo

    1993-12-31

    A modelling study was carried out to understand the effect of rock alteration on uranium radionuclide concentrations in rocks, in the vicinity of the Koongarra ore deposit, Australia. The one-dimensional, advection-dispersion-sorption model considers two important factors resulting from the process, mechanism and rate of chlorite weathering, one type of rock alteration that has occurred over a million-year period; (a) the changes in this distribution coefficients of uranium and thorium over time, and (b) that in rock porosity. The distribution coefficient and rock porosity at a given time, are expressed as average values of those of the coexisting minerals. By assuming a Darcy velocity of 0.9 m/y and an initial uranium concentration in the groundwater of 0.5 ppm, similar to present-day values, we were able to derive, from the finite-element modelling, uranium concentrations in the rock which are in good agreement with observed values. The calculated values of {sup 230}Th,{sup 234}U activity ratios show a similar trend to those observed, although more experimental data are necessary to confirm the similarity. However, in the absence of the rock alteration, the calculated results did not agree with those observed. The present study suggests that rock alteration should be included in models used to predict uranium migration over long-geologic timescales.

  2. APPLICATIONS OF BOREHOLE-ACOUSTIC METHODS IN ROCK MECHANICS.

    USGS Publications Warehouse

    Paillet, Frederick L.

    1985-01-01

    Acoustic-logging methods using a considerable range of wavelengths and frequencies have proven very useful in the in situ characterization of deeply buried crystalline rocks. Seismic velocities are useful in investigating the moduli of unfractured rock, and in producing a continuous record of rock quality for comparison with discontinuous intervals of core. The considerable range of frequencies makes the investigation of scale effects possible in both fractured and unfractured rock. Several specific methods for the characterization of in situ permeability have been developed and verified in the field.

  3. Estimating changes in rock permeability due to thermal-mechanical effects

    SciTech Connect

    Wang, H.F.; Blair, S.C.; Berge, P.A.

    1997-10-01

    This paper presents results of a modeling study of changes in fracture permeability due to thermal-mechanical effects associated with the potential geological repository at Yucca Mountain. a methodology for estimating changes in permeability is developed and applied to the Drift Scale Test (DST) now being conducted in the Exploratory Studies Facility (ESF) at Yucca Mountain. Temperature, stress, and displacement of rock in the heated zone are presented along with predicted zones where slip on fractures may occur. The zones of predicted fracture slip are used as a basis for predicting where permeability may be changed. this new procedure goes beyond previous models that relate stress to strain or displacement, and provides information about rock response that is needed for design of future tests at Yucca Mountain. Our results also contribute to the understanding of coupled processes in the near-field environment of a repository.

  4. An experimental study on fracture mechanical behavior of rock-like materials containing two unparallel fissures under uniaxial compression

    NASA Astrophysics Data System (ADS)

    Huang, Yan-Hua; Yang, Sheng-Qi; Tian, Wen-Ling; Zeng, Wei; Yu, Li-Yuan

    2016-06-01

    Strength and deformability characteristics of rock with pre-existing fissures are governed by cracking behavior. To further research the effects of pre-existing fissures on the mechanical properties and crack coalescence process, a series of uniaxial compression tests were carried out for rock-like material with two unparallel fissures. In the present study, cement, quartz sand, and water were used to fabricate a kind of brittle rock-like material cylindrical model specimen. The mechanical properties of rock-like material specimen used in this research were all in good agreement with the brittle rock materials. Two unparallel fissures (a horizontal fissure and an inclined fissure) were created by inserting steel during molding the model specimen. Then all the pre-fissured rock-like specimens were tested under uniaxial compression by a rock mechanics servo-controlled testing system. The peak strength and Young's modulus of pre-fissured specimen all first decreased and then increased when the fissure angle increased from 0° to 75°. In order to investigate the crack initiation, propagation and coalescence process, photographic monitoring was adopted to capture images during the entire deformation process. Moreover, acoustic emission (AE) monitoring technique was also used to obtain the AE evolution characteristic of pre-fissured specimen. The relationship between axial stress, AE events, and the crack coalescence process was set up: when a new crack was initiated or a crack coalescence occurred, the corresponding axial stress dropped in the axial stress-time curve and a big AE event could be observed simultaneously. Finally, the mechanism of crack propagation under microscopic observation was discussed. These experimental results are expected to increase the understanding of the strength failure behavior and the cracking mechanism of rock containing unparallel fissures.

  5. Deformation and stabilisation mechanisms of slow rock slides in crystalline bedrock

    NASA Astrophysics Data System (ADS)

    Zangerl, C.; Prager, C.

    2009-04-01

    . On a regional scale several valleys located in amphibolites, ortho- and paragneisses of the Ötztal-Stubai crystalline basement (i.e. Kaunertal, Pitztal, Ötztal, Lüsenstal, all located in North Tyrol, Austria) were investigated. Therefore geological and morphological basis data were compiled and re-evaluated, remote sensing methods (i.e. airborne laser scanning terrain models and orthofotos) applied and field mapping campaigns performed. On a local scale several rock slides were investigated and analysed in high detail with regard to their lithological and structural inventory, geometry of sliding masses and -zones, failure mechanisms, kinematics and temporal deformation characteristics. Field data clearly show that competent rock masses, e.g. orthogneisses and amphibolites, are affected by rapid failure events and therefore are characterised by "brittle" rock mass behaviour. In contrast, the majority of the slowly moving and "self-stabilising" rock slides are located totally or partly in mica-rich incompetent crystalline rock masses, e.g. paragneisses and micaschists, and are characterised by moderately dipping sliding zones. Apart from a causal lithological influence, numerous field observations demonstrate a major influence of pre-existing geological structures on the formation and deformation behaviour of these rock slides. The nature of rock slides implies that the temporal deformation behaviour is primarily dominated by two key-features of the sliding zone i.e. the mechanical properties (shear strain strengthening or weakening) and the effective in-situ stresses. The in-situ stresses along a sliding zone are influenced by the geometry of both the sliding mass and sliding zone, the internal deformation of the sliding mass and the pore pressures. All these properties can vary during progressive shear displacements. Especially large shear displacements in the range of tens to hundreds of metres along a distinct sliding zone can cause significant in-situ stress

  6. Mechanical Behaviour of Reservoir Rock Under Brine Saturation

    NASA Astrophysics Data System (ADS)

    Shukla, Richa; Ranjith, P. G.; Choi, S. K.; Haque, A.; Yellishetty, Mohan; Hong, Li

    2013-01-01

    Acoustic emissions (AE) and stress-strain curve analysis are well accepted ways of analysing crack propagation and monitoring the various failure stages (such as crack closure, crack initiation level during rock failure under compression) of rocks and rock-like materials. This paper presents details and results of experimental investigations conducted for characterizing the brittle failure processes induced in a rock due to monocyclic uniaxial compression on loading of two types of sandstone core samples saturated in NaCl brines of varying concentration (0, 2, 5, 10 and 15 % NaCl by weight). The two types of sandstone samples were saturated under vacuum for more than 45 days with the respective pore fluid to allow them to interact with the rocks. It was observed that the uniaxial compressive strength and stress-strain behaviour of the rock specimens changed with increasing NaCl concentration in the saturating fluid. The acoustic emission patterns also varied considerably for increasing ionic strength of the saturating brines. These observations can be attributed to the deposition of NaCl crystals in the rock's pore spaces as well some minor geo-chemical interactions between the rock minerals and the brine. The AE pattern variations could also be partly related to the higher conductivity of the ionic strength of the high-NaCl concentration brine as it is able to transfer more acoustic energy from the cracks to the AE sensors.

  7. Analysis of a microcrack model and constitutive equations for time-dependent dilatancy of rocks

    NASA Astrophysics Data System (ADS)

    Chen, Zuan

    2003-11-01

    Based on experimental observations and theoretical analyses, the author introduces an ideal microcrack model in which an array of cracks with the same shape and initial size is distributed evenly in rocks. The mechanism of creep dilatancy for rocks is analysed theoretically. Initiation, propagation and linkage of pre-existing microcracks during creep are well described. Also, the relationship between the velocity of microcrack growth and the duration of the creep process is derived numerically. The relationship agrees well with the character of typical experimental creep curves, and includes three stages of creep. Then the damage constitutive equations and damage evolution equations, which describe the dilatant behaviour of rocks, are presented. Because the dilatant estimated value is taken as the damage variable, the relationship between the microscopic model and the macroscopic constitutive equations is established. In this way the mechanical behaviour of rocks can be predicted.

  8. Event triggered data acquisition in the Rock Mechanics Laboratory

    SciTech Connect

    Hardy, R.D.

    1993-03-01

    Increasing complexity of experiments coupled with limitations of the previously used computers required improvements in both hardware and software in the Rock Mechanics Laboratories. Increasing numbers of input channels and the need for better graphics could no longer be supplied by DATAVG, an existing software package for data acquisition and display written by D. J. Holcomb in 1983. After researching the market and trying several alternatives, no commercial program was found which met our needs. The previous version of DATAVG had the basic features needed but was tied to obsolete hardware. Memory limitations on the previously used PDP-11 made it impractical to upgrade the software further. With the advances in IBM compatible computers it is now desirable to use them as data recording platforms. With this information in mind, it was decided to write a new version of DATAVG which would take advantage of newer hardware. The new version had to support multiple graphic display windows and increased channel counts. It also had to be easier to use.

  9. A THC Simulator for Modeling Fluid-Rock Interactions

    NASA Astrophysics Data System (ADS)

    Hamidi, Sahar; Galvan, Boris; Heinze, Thomas; Miller, Stephen

    2014-05-01

    Fluid-rock interactions play an essential role in many earth processes, from a likely influence on earthquake nucleation and aftershocks, to enhanced geothermal system, carbon capture and storage (CCS), and underground nuclear waste repositories. In THC models, two-way interactions between different processes (thermal, hydraulic and chemical) are present. Fluid flow influences the permeability of the rock especially if chemical reactions are taken into account. On one hand solute concentration influences fluid properties while, on the other hand, heat can affect further chemical reactions. Estimating heat production from a naturally fractured geothermal systems remains a complex problem. Previous works are typically based on a local thermal equilibrium assumption and rarely consider the salinity. The dissolved salt in fluid affects the hydro- and thermodynamical behavior of the system by changing the hydraulic properties of the circulating fluid. Coupled thermal-hydraulic-chemical models (THC) are important for investigating these processes, but what is needed is a coupling to mechanics to result in THMC models. Although similar models currently exist (e.g. PFLOTRAN), our objective here is to develop algorithms for implementation using the Graphics Processing Unit (GPU) computer architecture to be run on GPU clusters. To that aim, we present a two-dimensional numerical simulation of a fully coupled non-isothermal non-reactive solute flow. The thermal part of the simulation models heat transfer processes for either local thermal equilibrium or nonequilibrium cases, and coupled to a non-reactive mass transfer described by a non-linear diffusion/dispersion model. The flow process of the model includes a non-linear Darcian flow for either saturated or unsaturated scenarios. For the unsaturated case, we use the Richards' approximation for a mixture of liquid and gas phases. Relative permeability and capillary pressure are determined by the van Genuchten relations

  10. Sensitivity Analysis of Mechanical Parameters of Different Rock Layers to the Stability of Coal Roadway in Soft Rock Strata

    PubMed Central

    Zhao, Zeng-hui; Wang, Wei-ming; Gao, Xin; Yan, Ji-xing

    2013-01-01

    According to the geological characteristics of Xinjiang Ili mine in western area of China, a physical model of interstratified strata composed of soft rock and hard coal seam was established. Selecting the tunnel position, deformation modulus, and strength parameters of each layer as influencing factors, the sensitivity coefficient of roadway deformation to each parameter was firstly analyzed based on a Mohr-Columb strain softening model and nonlinear elastic-plastic finite element analysis. Then the effect laws of influencing factors which showed high sensitivity were further discussed. Finally, a regression model for the relationship between roadway displacements and multifactors was obtained by equivalent linear regression under multiple factors. The results show that the roadway deformation is highly sensitive to the depth of coal seam under the floor which should be considered in the layout of coal roadway; deformation modulus and strength of coal seam and floor have a great influence on the global stability of tunnel; on the contrary, roadway deformation is not sensitive to the mechanical parameters of soft roof; roadway deformation under random combinations of multi-factors can be deduced by the regression model. These conclusions provide theoretical significance to the arrangement and stability maintenance of coal roadway. PMID:24459447

  11. Microstructures and flow mechanisms in regional metamorphic rocks of Japan

    NASA Astrophysics Data System (ADS)

    Toriumi, Mitsuhiro; Teruya, Jun; Masui, Megumi; Kuwahara, Hidesato

    1986-09-01

    A number of microstructural features indicate a difference in the dominant deformation mechanism between the higher temperature Ryoke and the lower temperature Sambagawa and Shimanto metamorphic belts of Japan. The microstructures of metacherts containing deformed radiolaria are divided into two types: in both the Sambagawa and Shimanto belts the quartz grains are tabular while in the Ryoke belt they are equiaxed. TEM studies of these metacherts revealed that the tabular grains contain abundant subboundaries consisting of large numbers of network dislocations and bowe-out dislocations, while the equiaxed grains contain no subboundaries and have low densities of dislocations which are not bowed-out. There is a corresponding difference in the textures (lattice preferred orientation of quartz): the Ryoke metacherts display randomly distributed c-axes of quartz, while the Sambagawa and Shimanto metacherts show conspicuous crossed girdle patterns with some asymmetry. There is a third difference between these regions: in the metacherts of the Ryoke metamorphic belt, the strain magnitudes determined from deformed radiolaria increase with increasing volume fraction of mica in the same metamorphic P and T conditions, while in the Sambagawa and the Shimanto metamorphic cherts the strain magnitudes decrease with increasing the mica fraction. These microstructures, textures, and rheological behaviours of quartz-mica rocks suggest a change of deformation mechanism between the lower temperature Sambagawa and Shimanto, and the higher temperature Ryoke metamorphic belts. Since random fabrics of c-axes of quartz are inconsistent with lattice rotation due to dislocation glide, the Ryoke metacherts may have deformed by pressure-solution.

  12. The geology and mechanics of formation of the Fort Rock Dome, Yavapai County, Arizona

    USGS Publications Warehouse

    Fuis, Gary S.

    1996-01-01

    The Fort Rock Dome, a craterlike structure in northern Arizona, is the erosional product of a circular domal uplift associated with a Precambrian shear zone exposed within the crater and with Tertiary volcanism. A section of Precambrian to Quaternary rocks is described, and two Tertiary units, the Crater Pasture Formation and the Fort Rock Creek Rhyodacite, are named. A mathematical model of the doming process is developed that is consistent with the history of the Fort Rock Dome.

  13. Core-log integration for rock mechanics using borehole breakouts and rock strength experiments: Recent results from plate subduction margins

    NASA Astrophysics Data System (ADS)

    Saito, S.; Lin, W.

    2014-12-01

    Core-log integration has been applied for rock mechanics studies in scientific ocean drilling since 2007 in plate subduction margins such as Nankai Trough, Costa Rica margin, and Japan Trench. State of stress in subduction wedge is essential for controlling dynamics of plate boundary fault. One of the common methods to estimate stress state is analysis of borehole breakouts (drilling induced borehole wall compressive failures) recorded in borehole image logs to determine the maximum horizontal principal stress orientation. Borehole breakouts can also yield possible range of stress magnitude based on a rock compressive strength criterion. In this study, we constrained the stress magnitudes based on two different rock failure criteria, the Mohr-Coulomb (MC) criteria and the modified Wiebols-Cook (mWC) criteria. As the MC criterion is the same as that under unconfined compression state, only one rock parameter, unconfined compressive strength (UCS) is needed to constrain stress magnitudes. The mWC criterion needs the UCS, Poisson's ratio and internal frictional coefficient determined by triaxial compression experiments to take the intermediate principal stress effects on rock strength into consideration. We conducted various strength experiments on samples taken during IODP Expeditions 334/344 (Costa Rica Seismogenesis Project) to evaluate reliable method to estimate stress magnitudes. Our results show that the effects of the intermediate principal stress on the rock compressive failure occurred on a borehole wall is not negligible.

  14. DEM modeling of fracture propagation in veined rock

    NASA Astrophysics Data System (ADS)

    Virgo, S.; Abe, S.; Urai, J. L.

    2012-04-01

    One fundamental aspect of crack seal veins is that an existing vein can act as a heterogeneity in the rock which controls the localization of successive fracturing at unchanged mean stress orientations. Observations from crack-seal vein systems suggest that existing veins fundamentally influence the fracture behavior of a rock even in cases where the orientation of the stress field is highly incompatible with the orientation of the vein. We used a series of 3D Discrete Element Simulations to systematically investigate the influence of existing veins with varying orientation and mechanical properties on an approaching fracture. The models consist of a tabular heterogeneity within a bonded particle volume fractured under uniaxial tension. The parameters varied in the study are the orientation of the heterogeneity relative to the direction of uniaxial extension and therefore relative to the orientation of the favorable fracture plane as well as the fracture strength ratio between the matrix material, the vein material and the interface between vein and matrix material. The elastic parameters (e.g. Young's modulus) are kept homogeneous throughout the model. Thereby it is ensured that the results are not altered by stress field perturbation induced by stiffness contrasts. The model materials used were carefully tested and calibrated to ensure comparability with natural examples in terms of their fracture-mechanical properties. The simulations were repeated for several random particle packings to eliminate the effect of heterogeneities in the packing on the results. The results show a strong influence of the tabular heterogeneity on the fracture propagation for all orientations and at cohesion ratios within the range of natural systems. Besides curving and deflection of the fracture path associated with changes in fracture mode, bifurcation of fractures as well as arrest of propagation and nucleation of new fractures can be observed.

  15. Geophysical analysis of rock glacier internal structure and implications for deformation mechanics

    NASA Astrophysics Data System (ADS)

    Florentine, C. E.; Skidmore, M. L.; Speece, M. A.; Link, C. A.; Locke, W. W.; Carr, C. G.; Shaw, C. A.

    2011-12-01

    invoked to explain either a direct connection between individual transverse ridges to sub-surface structures or a specific structural regime. Our passive roof duplex faulting interpretation of GPR data at the LPRG is consistent with findings from previous studies on the internal composition and structure of rock glaciers and thus provides a testable model for improved understanding of rock glacier deformation mechanics.

  16. A coupled model of fluid flow in jointed rock

    SciTech Connect

    Swenson, Daniel; Martineau, Rick; James, Mark; Brown, Don

    1991-01-01

    We present a fully coupled model of fluid flow in jointed rock, where the fluid flow depends on the joint openings and the joint openings depend on the fluid pressure. The joints and rock blocks are modeled discretely using the finite element method. Solutions for the fluid and rock are obtained and iteration is performed until both solutions converge. Example applications include an examination of the effects of back-pressure on flow in a geothermal reservoir and transient fluid injection into a reservoir.

  17. MODELING UNDERGROUND STRUCTURE VULNERABILITY IN JOINTED ROCK

    SciTech Connect

    R. SWIFT; D. STEEDMAN

    2001-02-01

    The vulnerability of underground structures and openings in deep jointed rock to ground shock attack is of chief concern to military planning and security. Damage and/or loss of stability to a structure in jointed rock, often manifested as brittle failure and accompanied with block movement, can depend significantly on jointed properties, such as spacing, orientation, strength, and block character. We apply a hybrid Discrete Element Method combined with the Smooth Particle Hydrodynamics approach to simulate the MIGHTY NORTH event, a definitive high-explosive test performed on an aluminum lined cylindrical opening in jointed Salem limestone. Representing limestone with discrete elements having elastic-equivalence and explicit brittle tensile behavior and the liner as an elastic-plastic continuum provides good agreement with the experiment and damage obtained with finite-element simulations. Extending the approach to parameter variations shows damage is substantially altered by differences in joint geometry and liner properties.

  18. Determination of basic physical and mechanical properties of basaltic rocks from P-wave velocity

    NASA Astrophysics Data System (ADS)

    Karakuş, Askeri; Akatay, Mahmut

    2013-12-01

    Physical and mechanical properties of basaltic rocks used as main building material in historical buildings in Diyarbakir show great diversity depending on the place of origin. Especially, earthquake studies as well as restoration jobs and civil engineers and architects who work on building dynamics need to know basic material properties of basaltic rocks that are the main building material. In this study, the basalt samples obtained from 18 different locations of the Diyarbakir area were tested in order to estimate the main material properties of basalts used in historical buildings without collecting samples from them. Subsequently, statistical relationships between the nondestructive P-wave velocity and other properties of basalts were investigated. Consequently, highly correlated models (R2 = 0.717-0.890) were obtained between P-wave velocity and density, porosity, uniaxial compressive strength, Brazilian tensile strength, modulus of elasticity and Poisson's ratio.

  19. An aerodynamic model for one and two degree of freedom wing rock of slender delta wings

    NASA Technical Reports Server (NTRS)

    Hong, John

    1993-01-01

    The unsteady aerodynamic effects due to the separated flow around slender delta wings in motion were analyzed. By combining the unsteady flow field solution with the rigid body Euler equations of motion, self-induced wing rock motion is simulated. The aerodynamic model successfully captures the qualitative characteristics of wing rock observed in experiments. For the one degree of freedom in roll case, the model is used to look into the mechanisms of wing rock and to investigate the effects of various parameters, like angle of attack, yaw angle, displacement of the separation point, and wing inertia. To investigate the roll and yaw coupling for the delta wing, an additional degree of freedom is added. However, no limit cycle was observed in the two degree of freedom case. Nonetheless, the model can be used to apply various control laws to actively control wing rock using, for example, the displacement of the leading edge vortex separation point by inboard span wise blowing.

  20. Rock Degradation by Alkali Metals: A Possible Lunar Erosion Mechanism.

    PubMed

    Naughton, J J; Barnes, I L; Hammond, D A

    1965-08-01

    When rocks melt under ultrahigh-vacuum conditions, their alkali components volatilize as metals. These metal vapors act to comminute polycrystalline rocks to their component minerals. The resultant powder is porous and loosely packed and its characteristics may be compatible with the lunar surface as revealed by the Ranger photographs. If meteorite impact or lunar volcanism has produced vaporization or areas of molten lava, alkali erosion may have given dust of this character in adjacent solid areas. PMID:17747570

  1. Statistical mechanics of fragmentation processes of ice and rock bodies

    NASA Astrophysics Data System (ADS)

    Bashkirov, A. G.; Vityazev, A. V.

    1996-09-01

    It is a well-known experimental fact that impact fragmentation, specifically of ice and rock bodies, causes a two-step ("knee"-shaped) power distribution of fragment masses with exponent values within the limits -4 and -1.5 (here and henceforth the differential distribution is borne in mind). A new theoretical approach is proposed to determine the exponent values, a minimal fracture mass, and properties of the knee. As a basis for construction of non-equilibrium statistical mechanics of condensed matter fragmentation the maximum-entropy variational principle is used. In contrast to the usual approach founded on the Boltzmann entropy the more general Tsallis entropy allowing stationary solutions not only in the exponential Boltzmann-Gibbs form but in the form of the power (fractal) law distribution as well is invoked. Relying on the analysis of a lot of published experiments a parameter β is introduced to describe an inhomogeneous distribution of the impact energy over the target. It varies from 0 (for an utterly inhomogeneous distribution of the impact energy) to 1 (for a homogeneous distribution). The lower limit of fragment masses is defined as a characteristic fragment mass for which the energy of fragment formation is minimal. This mass value depends crucially on the value of β. It is shown that for β≪1 only small fragments can be formed, and the maximal permitted fragment (of mass m1) is the upper boundary of the first stage of the fracture process and the point where the knee takes place. The second stage may be realized after a homogeneous redistribution of the remainder of the impact energy over the remainder of the target (when β→1). Here, the formation of great fragments is permitted only and the smallest of them (of mass m2) determines a lower boundary of the second stage. Different forms of the knee can be observed depending on relations between m1 and m2.

  2. Linear Elastic and Cohesive Fracture Analysis to Model Hydraulic Fracture in Brittle and Ductile Rocks

    NASA Astrophysics Data System (ADS)

    Yao, Yao

    2012-05-01

    Hydraulic fracturing technology is being widely used within the oil and gas industry for both waste injection and unconventional gas production wells. It is essential to predict the behavior of hydraulic fractures accurately based on understanding the fundamental mechanism(s). The prevailing approach for hydraulic fracture modeling continues to rely on computational methods based on Linear Elastic Fracture Mechanics (LEFM). Generally, these methods give reasonable predictions for hard rock hydraulic fracture processes, but still have inherent limitations, especially when fluid injection is performed in soft rock/sand or other non-conventional formations. These methods typically give very conservative predictions on fracture geometry and inaccurate estimation of required fracture pressure. One of the reasons the LEFM-based methods fail to give accurate predictions for these materials is that the fracture process zone ahead of the crack tip and softening effect should not be neglected in ductile rock fracture analysis. A 3D pore pressure cohesive zone model has been developed and applied to predict hydraulic fracturing under fluid injection. The cohesive zone method is a numerical tool developed to model crack initiation and growth in quasi-brittle materials considering the material softening effect. The pore pressure cohesive zone model has been applied to investigate the hydraulic fracture with different rock properties. The hydraulic fracture predictions of a three-layer water injection case have been compared using the pore pressure cohesive zone model with revised parameters, LEFM-based pseudo 3D model, a Perkins-Kern-Nordgren (PKN) model, and an analytical solution. Based on the size of the fracture process zone and its effect on crack extension in ductile rock, the fundamental mechanical difference of LEFM and cohesive fracture mechanics-based methods is discussed. An effective fracture toughness method has been proposed to consider the fracture process zone

  3. Modeling of ductile deformation in anisotropic rocks with slip surfaces

    NASA Astrophysics Data System (ADS)

    Dabrowski, Marcin

    2013-04-01

    Flanking structures and sheath folds can develop in layered rocks due to flow perturbation around slip surfaces in shear zones (Exner and Dabrowski, 2010; Reber et al., submitted). Mechanical anisotropy of the host rock has been shown to play a major role in determining the slip rate and the flow pattern around it (Kocher and Mancktelow, 2006; Fletcher, 2011). In addition, anisotropic fluids such as ductile foliated rocks have a 'memory' of deformation due to evolving microstructure. For example, the rotation of a rigid circular inclusion embedded in a layered host in layer-parallel shear results in the structural reorganization around it, which leads to the modification of the flow pattern in the host and in consequence to a massive reduction of the inclusion rotation rate (Dabrowski and Schmid, 2011). Willis (1964) derived an analytical elastic solution for an elliptical inclusion in a homogeneous anisotropic matrix subject to a uniform load in the far field. The solution can be reduced to the case of an incompressible viscous medium. The case of an arbitrarily oriented inviscid slit under shear parallel to the principal axis of anisotropy can be obtained by reducing it even further. Although derived for the initial state of homogeneous planar anisotropy, the solution provides useful insights into the large deformation behavior of the system. In this study, I will use different models and numerical modeling techniques to assess the impact of mechanical anisotropy and structural development on the perturbing flow around an inviscid slit (slip surface) embedded in a host comprising discrete isotropic layers in layer-parallel simple shear. In the limit of thin layers (the number of layers intercepting the slit tends to infinity), the host is modeled as an anisotropic fluid. The anisotropic viscosity is determined by the bulk anisotropic viscosity of the layered system. The layering is initially planar or equivalently the anisotropy is initially homogeneous. Both non

  4. Is rock slope instability in high-mountain systems driven by topo-climatic, paraglacial or rock mechanical factors? - A question of scale!

    NASA Astrophysics Data System (ADS)

    Messenzehl, Karoline; Dikau, Richard

    2016-04-01

    Due to the emergent and (often non-linear) complex nature of mountain systems the key small-scale system properties responsible for rock slope instability contrast to those being dominant at larger spatial scales. This geomorphic system behaviour has major epistemological consequences for the study of rockfalls and associated form-process-relationships. As each scale requires its own scientific explanation, we cannot simply upscale bedrock-scale findings and, in turn, we cannot downscale the valley-scale knowledge to smaller phenomena. Here, we present a multi-scale study from the Turtmann Valley (Swiss Alps), that addresses rock slope properties at three different geomorphic levels: (i) regional valley scale, (ii) the hillslope scale and (iii) the bedrock scale. Using this hierarchical approach, we aim to understand the key properties of high-mountain systems responsible for rockfall initiation with respect to the resulting form-process-relationship at each scale. (i) At the valley scale (110 km2) rock slope instability was evaluated using a GIS-based modelling approach. Topo-climatic parameters, i.e. the permafrost distribution and the time since deglaciation after LGM were found to be the key variables causative for the regional-scale bedrock erosion and the storage of 62.3 - 65.3 x 106 m3 rockfall sediments in the hanging valleys (Messenzehl et al. 2015). (ii) At the hillslope scale (0.03 km2) geotechnical scanline surveys of 16 rock slopes and one-year rock temperature data of 25 ibuttons reveal that the local rockfall activity and the resulting deposition of individual talus slope landforms is mainly controlled by the specific rock mass strength with respect to the slope aspect, than being a paraglacial reaction. Permafrost might be only of secondary importance for the present-day rock mechanical state as geophysical surveys disprove the existence of frozen bedrock below 2600 m asl. (Messenzehl & Draebing 2015). (iii) At the bedrock scale (0.01 mm - 10 m) the

  5. Experimental study on the mechanical properties of simulated columnar jointed rock masses

    NASA Astrophysics Data System (ADS)

    Xiao, Wei-min; Deng, Rong-gui; Zhong, Zhi-bin; Fu, Xiao-min; Wang, Cong-yan

    2015-02-01

    Columnar jointed rock mass is a kind of structural rock mass commonly encountered in igneous rocks. Due to the effects of columnar joint networks, anisotropy is the typical mechanical property of columnar jointed rock mass, i.e. deformation and strength varying with loading direction. Correct understanding of the mechanical anisotropy of columnar jointed rock mass is a key problem that should be solved for demonstration and design of large scale rock mass projects such as dams and underground cavern excavations constructed in it. Plaster simulated columnar jointed rock mass specimens at dip angles varying from 0° to 90° with respect to the axial stress were tested under uniaxial compression conditions to investigate the mechanical anisotropy and failure modes. Based on analyses of experimental results, it was found that the strength and deformation of columnar jointed rock masses had pronounced ‘U-shaped’ anisotropy. In the anisotropic curves, the maximum and minimum values occurred at β = 90° and β = 45°, respectively. It was also shown that the lateral strain ratio was relatively high, especially when the dip angle was close to (45° - φj/2), where φj was the joint friction angle. An empirical expression was adopted to predict the ‘U-shaped’ anisotropy of deformation and strength and the predicted anisotropic curves agreed reasonably well with experimental data. Furthermore, four types of failure modes were summarized based on experimental results and corresponding mechanisms were also discussed.

  6. Damage-based long-term modelling of a large alpine rock slope

    NASA Astrophysics Data System (ADS)

    Riva, Federico; Agliardi, Federico; Amitrano, David; Crosta, Giovanni B.

    2016-04-01

    The morphology and stability of large alpine rock slopes result from the long-term interplay of different factors, following a complex history spanning several glacial cycles over thousands of years in changing morpho-climatic settings. Large rock slopes often experience slow long-term, creep-like movements interpreted as the macroscopic evidence of progressive failure in subcritically stressed rock masses. Slope damage and rock mass weakening associated to deglaciation are considered major triggers of these processes in alpine environments. Depending on rock mass properties, slope topography and removed ice thickness, valley flanks can progressively evolve over time into rockslides showing seasonal displacement trends, interpreted as evidence of hydro-mechanically coupled responses to hydrologic perturbations. The processes linking the long-term evolution of deglaciated rock slopes and their changing sensitivity to hydrologic triggers until rockslide failure, with significant implications in risk management and Early Warning, are not fully understood. We suggest that modelling long-term rock mass damage under changing conditions may provide such a link. We simulated the evolution of the Spriana rock slope (Italian Central Alps). This is affected by a 50 Mm3 rockslide, significantly active since the late 19th century and characterized by massive geological and geotechnical investigations and monitoring during the last decades. Using an improved version of the 2D Finite-Element, damage-based brittle creep model proposed by Amitrano and Helmstetter (2006) and Lacroix and Amitrano (2013), we combined damage and time-to-failure laws to reproduce diffused damage, strain localization and the long-term creep deformation of the slope. The model was implemented for application to real slopes, by accounting for: 1) fractured rock mass properties upscaling based on site characterization data; 2) fluid pressures in a progressive failure context, relating fluid occurrence to

  7. Digital Rock Physics: Mechanical Properties of Carbonate Core Plug at Different Resolutions

    NASA Astrophysics Data System (ADS)

    Jouini, M. S.; Faisal, T. F.; Islam, A.; Chevalier, S.; Jouiad, M.; Sassi, M.

    2014-12-01

    Digital Rock Physics (DRP) is a novel technology that could be used to generate accurate, fast and cost effective special core analysis (SCAL) properties to support reservoir characterization and simulation tools. For this work, Micro-CT images at different resolutions have been used to run simulations to determine elastic properties like bulk, shear, Young's Modulus and Poisson's ratio of a dry carbonate core plug from Abu Dhabi reservoirs. Pre processing and segmentation of raw images is performed in FEI 3D visualization and analysis tool Avizo. Carbonates are characterized by a very complex pore-space structure and so a high degree of heterogeneity. Abaqus that is based on Finite Element Method is used to run 2D and 3D elastic simulations. Results will be compared by simulating the same core-plug in an alternative segmentation and FEM modeling environment used previously by Jouini & Vega et al. 2012 [1]. Acoustic wave propagation experiments at different confining pressures are performed in the laboratory Triaxial machine to determine the dynamic Young's modulus and Poisson's ratio for the same core plug. Expeirmental results are compared with numerical results. [1] Jouini, M.S. and Vega, S. 2012. Simulation of carbonate rocks elastic properties using 3D X-Ray computed tomography images based on Discrete Element Method and Finite Element Method. 46th US Rock Mechanics / Geomechanics Symposium, Chicago, Il, USA, 24-27 June 2012.

  8. Study on the Constitutive Model for Jointed Rock Mass

    PubMed Central

    Xu, Qiang; Chen, Jianyun; Li, Jing; Zhao, Chunfeng; Yuan, Chenyang

    2015-01-01

    A new elasto-plastic constitutive model for jointed rock mass, which can consider the persistence ratio in different visual angle and anisotropic increase of plastic strain, is proposed. The proposed the yield strength criterion, which is anisotropic, is not only related to friction angle and cohesion of jointed rock masses at the visual angle but also related to the intersection angle between the visual angle and the directions of the principal stresses. Some numerical examples are given to analyze and verify the proposed constitutive model. The results show the proposed constitutive model has high precision to calculate displacement, stress and plastic strain and can be applied in engineering analysis. PMID:25885695

  9. Rock mass mechanical property estimations for the Yucca Mountain Site Characterization Project; Yucca Mountain Site Characterization Project

    SciTech Connect

    Lin, M.; Hardy, M.P.; Bauer, S.J.

    1993-06-01

    Rock mass mechanical properties are important in the design of drifts and ramps. These properties are used in evaluations of the impacts of thermomechanical loading of potential host rock within the Yucca Mountain Site Characterization Project. Representative intact rock and joint mechanical properties were selected for welded and nonwelded tuffs from the currently available data sources. Rock mass qualities were then estimated using both the Norwegian Geotechnical Institute (Q) and Geomechanics Rating (RMR) systems. Rock mass mechanical properties were developed based on estimates of rock mass quality, the current knowledge of intact properties, and fracture/joint characteristics. Empirical relationships developed to correlate the rock mass quality indices and the rock mass mechanical properties were then used to estimate the range of rock mass mechanical properties.

  10. A two-stage model of fracture of rocks

    USGS Publications Warehouse

    Kuksenko, V.; Tomilin, N.; Damaskinskaya, E.; Lockner, D.

    1996-01-01

    In this paper we propose a two-stage model of rock fracture. In the first stage, cracks or local regions of failure are uncorrelated occur randomly throughout the rock in response to loading of pre-existing flaws. As damage accumulates in the rock, there is a gradual increase in the probability that large clusters of closely spaced cracks or local failure sites will develop. Based on statistical arguments, a critical density of damage will occur where clusters of flaws become large enough to lead to larger-scale failure of the rock (stage two). While crack interaction and cooperative failure is expected to occur within clusters of closely spaced cracks, the initial development of clusters is predicted based on the random variation in pre-existing Saw populations. Thus the onset of the unstable second stage in the model can be computed from the generation of random, uncorrelated damage. The proposed model incorporates notions of the kinetic (and therefore time-dependent) nature of the strength of solids as well as the discrete hierarchic structure of rocks and the flaw populations that lead to damage accumulation. The advantage offered by this model is that its salient features are valid for fracture processes occurring over a wide range of scales including earthquake processes. A notion of the rank of fracture (fracture size) is introduced, and criteria are presented for both fracture nucleation and the transition of the failure process from one scale to another.

  11. Discrete element thermomechanical modelling of rock cutting with valuation of tool wear

    NASA Astrophysics Data System (ADS)

    Rojek, Jerzy

    2014-05-01

    The paper presents a thermomechanical discrete element model of rock cutting process. The thermomechanical formulation of the discrete element method considers mechanical and thermal phenomena and their reciprocal influence. The thermal model developed for transient heat conduction problems takes into account conductive heat transfer at the contact between particles and convection on the free surface. The thermal and mechanical problems are coupled by consideration of: (1) heat generated due to friction which is calculated in the mechanical problem and passed to the thermal solution, (2) influence of thermal expansion on mechanical interaction between particles. Estimation of temperature dependent wear has been included into the contact model. The coupled problem is solved using the staggered scheme.The thermomechanical algorithm has been implemented in a discrete element program and applied to simulation of rock cutting with single pick of a dredge cutter head. Numerical results confirm good performance of the developed algorithm.

  12. Mechanical Characteristics of rocks cored from Hanging Wall of Chelungpu Fault

    NASA Astrophysics Data System (ADS)

    Lu, C. Y.; Chen, C. W.; Hu, J. C.; Tsai, L. S.; Lin, M. L.; Jeng, F. S.

    2005-12-01

    Representatives mechanical parameters of the strata involved in fault movements are essentially needed when conducting numerical simulation of fault movements. To study the mechanical characteristic of the rocks, to conduct the subsequent numberical analyses, and to interpret the Chelungpu fault movements, this research systematically studied the mechanical properties of the rocks cored from the TCDP program. Accordingly, the samples sampled from the cores are subjected to a series of mechanical experiments, including stress-path controlled pure shear tests to study their physical properties as well as their strengths and deformability with elastic and plastic strained being distinguished. Meanwhile, temperature factor are also conducted so as to study its influence on the rocks. It was found that the rocks, with a depth ranging from 450 to 1300m, mainly comprise of silt stone, fossil-riched sandstone and sandstone, with a porosities of 2%, 6% and 15% and uniaxial compressive strengths of 67~73 MPa, 61~65 MPa and 8~11 MPa, respectively. Accordingly, we can find that the strength and deformation of rocks relate to their porosities. According to the pure shear path test results, it was also found that the behavior of the last two types rocks, including the non-linear elastic deformation, plastic strain locking and the elastic strain, was coupled with shear stress. As to the temperature , its influence on the rocks was found to be not significant. On the other hand, the experimental results indicate that the Drucker-Prager failure criteria can describe the failure envelope line of these three type rocks. The strength of these rocks is independent with stress-path. Finally, the constitutive parameters of these sandstones were obtained, which enables realistic prediction the deformational behavior of the rocks in the near future.

  13. A damage-softening statistical constitutive model considering rock residual strength

    NASA Astrophysics Data System (ADS)

    Wang, Zhi-liang; Li, Yong-chi; Wang, J. G.

    2007-01-01

    Under stress, the microcracks in rock evolve (initiation, growth and coalescence) from damage to fracture with a continuous process. In order to describe this continuous process, a damage-softening statistical constitutive model for rock was proposed based on the Weibull distribution of mesoscopic element strength. This model usually adopts the Drucker-Prager criterion as its distribution parameter of mesoscopic element strength, which may produce larger damage zone in numerical simulations. This paper mainly studies the effects of strength criteria and residual strength on the performance of this damage-softening statistical constitutive model of rock. Main works include following three aspects: Firstly, the mechanical behaviors of rock are comparatively studied when the Drucker-Prager and the Mohr-Coulomb criteria are employed, respectively, as the distribution parameter. Then, a coefficient is introduced to make this constitutive model be capable of describing the residual strength of rock. Finally, a user-defined subroutine is concisely developed for this model and checked through typical strain paths. The current work lays a good foundation for further application of this model in geotechnics and geosciences.

  14. Formation and exhumation mechanisms of high-grade rocks: Sagduction and Subduction processes during the Archean

    NASA Astrophysics Data System (ADS)

    François, C.; Philippot, P.; Rey, P.

    2012-04-01

    The interpretation of high-grade rocks in the Archean is controversial. Mid- to high-pressure assemblages are commonly interpreted in terms of plate tectonic processes including subduction. In the Archean however, mid- to high-pressure assemblages could have been also produced during the sagduction of greenstone covers into their crustal basement. Often put in opposition, sagduction and subduction are not incompatible processes. In order to better documents the P-T-t signatures of both processes we are conducting a comparative study - structural, metamorphic and numerical - of supposedly subduction-related metamorphic rocks described in ~3.5-3.2 Ga old Barberton greenstones (Kaapvaal Craton, South Africa) (Moyen et al., 2006), and supposedly sagduction-related high-grade rocks in the 3.5-3.2 Ga old East Pilbara Craton (Western Australia) (Delor et al., 1991). Interestingly, these two terranes display dome-and-keel structure in which narrow belts of greenstone (ultramafic and mafic metabasalts) and overlying sedimentary rocks occur in association with broad TTG (trondhjemite-tonalite-granodiorite) granitoids. We present here preliminary results from fieldwork, metamorphic investigations and numerical experiments. Petrological analyses have been conducted on metabasalts and metasediments in enclaves in migmatitic and granitic rocks, both inside and outside granitic domes. We sampled high-grade mafic rocks in enclaves within the NNE trending steeply deeping migmatitic Inyoni shear zone located between the 3.45 Ga Stolzburg pluton and the 3.2 Ga Badplass gneisses in the southern Barberton terrane (Moyen et al., 2006). Preliminary P-T estimations have been performed with multi-equilibrium approach using Thermocalc and with thermodynamic modeling using PerpleX on garnet-amphibole-clinopyroxene-epidote-plagioclase assemblage reveals pressures of 12-14 kbar at temperatures of 600-650°C for the metamorphic peak. Maximum temperature is reached at the beginning of exhumation

  15. Frictional sliding in layered rock model: Preliminary experiments. Yucca Mountain Site Characterization Project

    SciTech Connect

    Perry, K.E. Jr.; Buescher, B.J.; Anderson, D.; Epstein, J.S.

    1995-09-01

    An important aspect of determining the suitability of Yucca Mountain as a possible nuclear waste repository requires understanding the mechanical behavior of jointed rock-masses. To this end we have studied the frictional sliding between simulated rock joints in the laboratory using the technique of phase shifting moire interferometry. The models were made from stacks of Lexan plates and contained a central hole to induce slip between the plates when the models were loaded in compression. These preliminary results confirm the feasibility of the approach and show a clear evolution of slip as function of load.

  16. Rock fracture by ice segregation: linking laboratory modelling and rock slope erosion

    NASA Astrophysics Data System (ADS)

    Murton, J.

    2009-04-01

    It has been unclear until recently if ice can fracture intact bedrock subject to natural freezing regimes, or whether it simply enlarges existing fractures or does both. This question is important, because if ice segregation in bedrock permafrost is widespread, then there may be considerable potential for significantly increased rock slope instability as ice-cemented discontinuities warm and active layers thicken. Laboratory modelling has now begun to elucidate the process of ice segregation in bedrock. Laboratory experiments indicate that moist, porous rock behaves remarkably like moist, frost-susceptible soil, with both substrates experiencing ice enrichment and fracture / fissuring of near-surface permafrost. It appears that significant concentrations of segregated ice are most likely in the transition zone between the active layer and the permafrost, as a result of downward migration of water in summer and upward advance of freezing at the beginning of the winter. Laboratory modelling indicates that given adequate water supply, ice segregation produces a zone of ice-bonded fractured bedrock immediately below the permafrost table. In general, the importance of ice segregation relative to in situ volume expansion increases with decreasing thermal gradients and increasing duration of freezing. Recent theoretical developments suggest that the maximum possible disjointing pressure is governed by the temperature depression below the bulk-melting point, even in the absence of large temperature gradients, and therefore slow ice segregation in bedrock may be possible at greater depths where the frozen permeability of rock limits the actual amount of heave produced. Thus, over long timescales, ice segregation may be highly significant in frozen steep bedrock slopes where the presence of ice-rich fractured bedrock may be critically important in releasing rock falls and rock slides during climate-induced warming and permafrost degradation. With recent climate warming

  17. Rock-Mechanical Constraints on SPH Applications to Asteroid Impact Evolution

    NASA Astrophysics Data System (ADS)

    Bruesch, L.; Asphaug, E.

    2002-12-01

    The smooth particle hydrodynamics (SPH) code as adapted for dynamic brittle fragmentation (Benz and Asphaug 1994, 1995) has become a leading technique for modeling meteoroid collisions into asteroids with realistic geologies and shapes (Asphaug et al., Icarus 1996; Nature 1998). Together with earlier techniques relying on the same Weibull-based Grady-Kipp fracture model (e.g. Melosh et al. 1992), it has been used to establish that asteroids larger than a few hundred meters diameter are rubble piles (Benz and Asphaug, Icarus 1999), and is applied for learning how binary asteroids form during tidal events and collisions (Michel et al., Science 2001) and how craters and regolith form on irregular, rotating bodies. But all of these applications, especially when the outcome involves a consideration of mechanical strength, rely upon the assumption that flaws are distributed according to a Weibull distribution throughout a rock mass, and that those flaws are activated dynamically and relieve local stress in a circumscribing volume. Our SPH fragmentation code has been calibrated against a variety of laboratory impact experiments, but never for a suite of experiments spanning size scale and energy scale. It would be prudent, before relying upon model predictions for km-scale rock masses, to benchmark our code against laboratory experiments, in particular the controlled set of laboratory experiments of Housen and Holsapple (Icarus 1999) designed to examine the dependence of a body's strength on its size. On our new computers we are running suites of numerical simulations to reproduce these experiments in which our only varied parameters are the laboratory impact conditions themselves: projectile and target diameter, for identical rock types (granite and basalt).

  18. Coupling photogrammetric data with DFN-DEM model for rock slope hazard assessment

    NASA Astrophysics Data System (ADS)

    Donze, Frederic; Scholtes, Luc; Bonilla-Sierra, Viviana; Elmouttie, Marc

    2013-04-01

    Structural and mechanical analyses of rock mass are key components for rock slope stability assessment. The complementary use of photogrammetric techniques [Poropat, 2001] and coupled DFN-DEM models [Harthong et al., 2012] provides a methodology that can be applied to complex 3D configurations. DFN-DEM formulation [Scholtès & Donzé, 2012a,b] has been chosen for modeling since it can explicitly take into account the fracture sets. Analyses conducted in 3D can produce very complex and unintuitive failure mechanisms. Therefore, a modeling strategy must be established in order to identify the key features which control the stability. For this purpose, a realistic case is presented to show the overall methodology from the photogrammetry acquisition to the mechanical modeling. By combining Sirovision and YADE Open DEM [Kozicki & Donzé, 2008, 2009], it can be shown that even for large camera to rock slope ranges (tested about one kilometer), the accuracy of the data are sufficient to assess the role of the structures on the stability of a jointed rock slope. In this case, on site stereo pairs of 2D images were taken to create 3D surface models. Then, digital identification of structural features on the unstable block zone was processed with Sirojoint software [Sirovision, 2010]. After acquiring the numerical topography, the 3D digitalized and meshed surface was imported into the YADE Open DEM platform to define the studied rock mass as a closed (manifold) volume to define the bounding volume for numerical modeling. The discontinuities were then imported as meshed planar elliptic surfaces into the model. The model was then submitted to gravity loading. During this step, high values of cohesion were assigned to the discontinuities in order to avoid failure or block displacements triggered by inertial effects. To assess the respective role of the pre-existing discontinuities in the block stability, different configurations have been tested as well as different degree of

  19. Geotechnical modeling of high-level nuclear waste disposal by rock melting

    SciTech Connect

    Heuze, F.E.

    1981-12-01

    A new strategy has been developed for the geotechnical modeling of nuclear waste disposal by rock melting (DRM). Three seeparate tasks were performed to reach this objective: a review of the four scenarios which have been proposed for DRM, to date; an evaluation of computer-based numerical models which could be used to analyze the mechanical, thermal, and hydraulic processes involved in DRM; and a critical review of rock mass properties which are relevant to the design and safety of waste disposal by rock melting. It is concluded that several geotechnical aspects of DRM can be studied realistically with current state-of-the-art model capabilities and knowledge of material properties. The next step in the feasibility study of DRM should be a best-estimate calculation of the four cavity-melt and canister-burial concepts. These new analyses will indicate the most critical areas for subsequent research.

  20. Modeling rock specimens through 3D printing: Tentative experiments and prospects

    NASA Astrophysics Data System (ADS)

    Jiang, Quan; Feng, Xiating; Song, Lvbo; Gong, Yahua; Zheng, Hong; Cui, Jie

    2016-02-01

    Current developments in 3D printing (3DP) technology provide the opportunity to produce rock-like specimens and geotechnical models through additive manufacturing, that is, from a file viewed with a computer to a real object. This study investigated the serviceability of 3DP products as substitutes for rock specimens and rock-type materials in experimental analysis of deformation and failure in the laboratory. These experiments were performed on two types of materials as follows: (1) compressive experiments on printed sand-powder specimens in different shapes and structures, including intact cylinders, cylinders with small holes, and cuboids with pre-existing cracks, and (2) compressive and shearing experiments on printed polylactic acid cylinders and molded shearing blocks. These tentative tests for 3DP technology have exposed its advantages in producing complicated specimens with special external forms and internal structures, the mechanical similarity of its product to rock-type material in terms of deformation and failure, and its precision in mapping shapes from the original body to the trial sample (such as a natural rock joint). These experiments and analyses also successfully demonstrate the potential and prospects of 3DP technology to assist in the deformation and failure analysis of rock-type materials, as well as in the simulation of similar material modeling experiments.

  1. Brittle Rock Modeling Approach and its Validation Using Excavation-Induced Micro-Seismicity

    NASA Astrophysics Data System (ADS)

    Ma, Chun-Chi; Li, Tian-Bin; Xing, Hui-Lin; Zhang, Hang; Wang, Min-Jie; Liu, Tian-Yi; Chen, Guo-Qing; Chen, Zi-Quan

    2016-08-01

    With improvements to the bonded-particle model, a custom indicator of crack intensity is introduced to grade rock fractures accurately. Brittle fracturing of rock mass is studied using the bonded-particle model; here, "brittle" refers to the process where more energy is released towards making particles collide and disperse, and hence results in the quick emergence of "chain cracks". Certain principles concerning how to construct brittle rock are then proposed. Furthermore, a modeling approach for brittle rocks based on the adaptive continuum/discontinuum (AC/DC) method is proposed to aid the construction of large-scale models of tunnel excavations. To connect with actual tunneling conditions, fundamental mechanical properties, the mechanism for brittle fracturing, the joint distribution, and the initial stress field are considered in the modeling approach. Results from micro-seismic monitoring of a tunnel excavation confirmed the suitability of this modeling approach to simulate crack behavior, and results show that simulated cracking exhibit similar trends (evolution, location, and intensity) with micro-seismic cracking.

  2. VNIR spectral modeling of Mars analogue rocks: first results

    NASA Astrophysics Data System (ADS)

    Pompilio, L.; Roush, T.; Pedrazzi, G.; Sgavetti, M.

    Knowledge regarding the surface composition of Mars and other bodies of the inner solar system is fundamental to understanding of their origin, evolution, and internal structures. Technological improvements of remote sensors and associated implications for planetary studies have encouraged increased laboratory and field spectroscopy research to model the spectral behavior of terrestrial analogues for planetary surfaces. This approach has proven useful during Martian surface and orbital missions, and petrologic studies of Martian SNC meteorites. Thermal emission data were used to suggest two lithologies occurring on Mars surface: basalt with abundant plagioclase and clinopyroxene and andesite, dominated by plagioclase and volcanic glass [1,2]. Weathered basalt has been suggested as an alternative to the andesite interpretation [3,4]. Orbital VNIR spectral imaging data also suggest the crust is dominantly basaltic, chiefly feldspar and pyroxene [5,6]. A few outcrops of ancient crust have higher concentrations of olivine and low-Ca pyroxene, and have been interpreted as cumulates [6]. Based upon these orbital observations future lander/rover missions can be expected to encounter particulate soils, rocks, and rock outcrops. Approaches to qualitative and quantitative analysis of remotely-acquired spectra have been successfully used to infer the presence and abundance of minerals and to discover compositionally associated spectral trends [7-9]. Both empirical [10] and mathematical [e.g. 11-13] methods have been applied, typically with full compositional knowledge, to chiefly particulate samples and as a result cannot be considered as objective techniques for predicting the compositional information, especially for understanding the spectral behavior of rocks. Extending the compositional modeling efforts to include more rocks and developing objective criteria in the modeling are the next required steps. This is the focus of the present investigation. We present results of

  3. Effective stress model for partially and fully saturated rocks

    SciTech Connect

    Dey, T.N.

    1989-01-01

    An effective stress model which calculates the pressure-volume (P-V) and deviatoric stress response of partially and fully saturated rocks is described here. The model includes pore pressure effects on pore crushing and shear strength as well as effects of shear enhanced void collapse and shear caused dilatancy. The model can directly use tabular data for the P-V behavior of the rock solids and the water, and for the drained pore crushing behavior and shear strength, which simplifies model fitting. Phase transitions in the solids and vaporization of the water are also allowed. Use of the model is illustrated by an example of wave propagation in limestone. 6 refs., 4 figs.

  4. The Evolution of Fracture Systems in Rocks with Veins: Insights from 3D Discrete Element Models

    NASA Astrophysics Data System (ADS)

    Virgo, S.; Urai, J. L.; Abe, S.

    2014-12-01

    Observations from natural vein systems suggest that preexisting veins can strongly influence orientation, continuity and connectivity of fractures in a rock even in cases where the orientation of the veins is incompatible with the orientation of the stress field. We present a numerical method to model cycles of fracturing and sealing in a rotating stress field to simulate such systems, for different strength ratios of host rock and vein. We study a layered model under vertical stress and uniaxial horizontal extension. This represents common conditions in sedimentary basins with layers of varying composition. The model with fractures that form during the first deformation phase is sealed and deformed again in a different direction to model the effect of a changing horizontal stress field. We find different types of fracture interaction with veins, depending on the strength contrast between veins and host rock and amount of rotation. The crack-seal and crack-jump mechanisms ensue naturally from the models as a result of the strength of the vein material relative to the host rock. Weak veins localize fracturing and reactivate, even in high misorientation to the extension direction. Connecting fractures between reactivated veins form at a higher angle to the veins than expected. In these systems, the connectivity of the fracture network is dramatically increased. Veins stronger than the host rock have less influence on the new fractures. Most fractures crosscut the veins by the step-over mechanism. Deflection occurs for favorable vein orientations but the deflection length is very short. The results are in good agreement with natural crack seal vein networks found in carbonate rocks of the Oman Mountains. We find that preexisting veins can change the fracture behavior of a rock in a way that new fractures do not necessarily align with the principle extension direction and form a highly connected network with reactivated veins that dramatically enhances lateral

  5. Numerical investigation of the hydro-mechanical contribution to seismic attenuation in damaged rocks

    NASA Astrophysics Data System (ADS)

    Pollmann, Nele; Jänicke, Ralf; Renner, Jörg; Steeb, Holger

    2016-04-01

    The investigation of hydro-mechanical processes, in particular the modeling of seismic waves in fractured porous media, is essential for the physical interpretation of data obtained from seismic exploration. Here, we specifically investigate attenuation processes in fluid-saturated porous rock containing fracture networks to identify effective hydro-mechanical properties by numerical simulation. The main purpose of this work is the characterization of the overall hydro-mechanical properties by computational homogenization. We determine an effective Skempton coefficient by investigating the fluid pressure and the solid displacement of the skeleton saturated by compressible fluids. Fracture networks are stochastically generated to mimic geological in-situ situations. The fractures are approximated as ellipses with aspect ratios up to 1/100, i.e. they constitute thin and long hydraulic conduits with high permeabilities. Simulations are designed on the material scale with and without conservation of fluid mass in the control volume. Using computational homogenization approaches, we define an effective Skempton coefficient. A range of fracture networks with different characteristic properties is studied for different varieties of fractures. On the material scale we find strongly heterogeneous pressure propagation in the fracture network and the surrounding rock, respectively. The pressure diffusion is much faster in the fracture network than in the matrix, rendering the macroscopic hydro-mechanical behavior strongly time dependent. The effective Skempton coefficient converges to an ensemble-specific instantaneous value and to 1 for long-time studies. The ultimate objective of our study is to evaluate whether constraints on the structure of fracture networks can be deduced from observations of attenuation and its frequency dependence.

  6. Hydrological modelling in sandstone rocks watershed

    NASA Astrophysics Data System (ADS)

    Ponížilová, Iva; Unucka, Jan

    2015-04-01

    The contribution is focused on the modelling of surface and subsurface runoff in the Ploučnice basin. The used rainfall-runoff model is HEC-HMS comprising of the method of SCS CN curves and a recession method. The geological subsurface consisting of sandstone is characterised by reduced surface runoff and, on the contrary, it contributes to subsurface runoff. The aim of this paper is comparison of the rate of influence of sandstone on reducing surface runoff. The recession method for subsurface runoff was used to determine the subsurface runoff. The HEC-HMS model allows semi- and fully distributed approaches to schematisation of the watershed and rainfall situations. To determine the volume of runoff the method of SCS CN curves is used, which results depend on hydrological conditions of the soils. The rainfall-runoff model assuming selection of so-called methods of event of the SCS-CN type is used to determine the hydrograph and peak flow rate based on simulation of surface runoff in precipitation exceeding the infiltration capacity of the soil. The recession method is used to solve the baseflow (subsurface) runoff. The method is based on the separation of hydrograph to direct runoff and subsurface or baseflow runoff. The study area for the simulation of runoff using the method of SCS CN curves to determine the hydrological transformation is the Ploučnice basin. The Ploučnice is a hydrologically significant river in the northern part of the Czech Republic, it is a right tributary of the Elbe river with a total basin area of 1.194 km2. The average value of CN curves for the Ploučnice basin is 72. The geological structure of the Ploučnice basin is predominantly formed by Mesozoic sandstone. Despite significant initial loss of rainfall the basin response to the causal rainfall was demonstrated by a rapid rise of the surface runoff from the watershed and reached culmination flow. Basically, only surface runoff occures in the catchment during the initial phase of

  7. Micromechanical Tests and Geochemical Modeling to Evaluate Evolution of Rock Alteration by CO2-Water Mixtures

    NASA Astrophysics Data System (ADS)

    Aman, M.; Sun, Y.; Ilgen, A.; Espinoza, N.

    2015-12-01

    Injection of large volumes of CO2 into geologic formations can help reduce the atmospheric CO2 concentration and lower the impact of burning fossil fuels. However, the injection of CO2 into the subsurface shifts the chemical equilibrium between the mineral assemblage and the pore fluid. This shift will situationally facilitate dissolution and reprecipitation of mineral phases, in particular intergranular cements, and can potentially affect the long term mechanical stability of the host formation. The study of these coupled chemical-mechanical reservoir rock responses can help identify and control unexpected emergent behavior associated with geological CO2 storage.Experiments show that micro-mechanical methods are useful in capturing a variety of mechanical parameters, including Young's modulus, hardness and fracture toughness. In particular, micro-mechanical measurements are well-suited for examining thin altered layers on the surfaces of rock specimens, as well as capturing variability on the scale of lithofacies. We performed indentation and scratching tests on sandstone and siltstone rocks altered in natural CO2-brine environments, as well as on analogous samples altered under high pressure, temperature, and dissolved CO2 conditions in a controlled laboratory experiment. We performed geochemical modeling to support the experimental observations, in particular to gain the insight into mineral dissolution/precipitation as a result of the rock-water-CO2reactions. The comparison of scratch measurements performed on specimens both unaltered and altered by CO2 over geologic time scales results in statistically different values for fracture toughness and scratch hardness, indicating that long term exposure to CO2 caused mechanical degradation of the reservoir rock. Geochemical modeling indicates that major geochemical change caused by CO2 invasion of Entrada sandstone is dissolution of hematite cement, and its replacement with siderite and dolomite during the

  8. Influence of Water Content on the Mechanical Properties of an Argillaceous Swelling Rock

    NASA Astrophysics Data System (ADS)

    Vergara, Maximiliano R.; Triantafyllidis, Theodoros

    2016-07-01

    This work presents the results of an experimental investigation aimed to study the effect of water on the mechanical properties of a partially saturated argillaceous swelling rock. The mineralogical composition of the rock, in particular the portion of swelling clays, was determined with X-ray diffraction. The water retention curve was estimated using a dew-point potential meter and the swelling behavior of the studied rock was examined by unconfined and oedometric swelling tests. The influence of water on the rock mechanical properties was assessed by means of triaxial tests. The experimental results indicate a strong decrease of strength and stiffness with increasing saturation or decreasing suction. This occurs only within a certain range of saturation. Degradation of the rock properties can be expected for small increments in the water content within this range. At low suction and close to the air-entry value, the stiffness remained constant. As the rock desaturates, the strength and stiffness increase approaching constant values. For suction greater than about 76 MPa, low increase of strength and stiffness was observed. The specimens in the swelling tests reached a saturation degree of 70 % which corresponds to a decrease of strength and stiffness of approximately 80 %. Rock swelling occurring simultaneously with reduction of strength and stiffness, increases deformations and it is an important issue for the stability of excavations.

  9. Multiphysics processes in partially saturated fractured rock: Experiments and models from Yucca Mountain

    NASA Astrophysics Data System (ADS)

    Rutqvist, Jonny; Tsang, Chin-Fu

    2012-09-01

    The site investigations at Yucca Mountain, Nevada, have provided us with an outstanding data set, one that has significantly advanced our knowledge of multiphysics processes in partially saturated fractured geological media. Such advancement was made possible, foremost, by substantial investments in multiyear field experiments that enabled the study of thermally driven multiphysics and testing of numerical models at a large spatial scale. The development of coupled-process models within the project have resulted in a number of new, advanced multiphysics numerical models that are today applied over a wide range of geoscientific research and geoengineering applications. Using such models, the potential impact of thermal-hydrological-mechanical (THM) multiphysics processes over the long-term (e.g., 10,000 years) could be predicted and bounded with some degree of confidence. The fact that the rock mass at Yucca Mountain is intensively fractured enabled continuum models to be used, although discontinuum models were also applied and are better suited for analyzing some issues, especially those related to predictions of rockfall within open excavations. The work showed that in situ tests (rather than small-scale laboratory experiments alone) are essential for determining appropriate input parameters for multiphysics models of fractured rocks, especially related to parameters defining how permeability might evolve under changing stress and temperature. A significant laboratory test program at Yucca Mountain also made important contributions to the field of rock mechanics, showing a unique relation between porosity and mechanical properties, a time dependency of strength that is significant for long-term excavation stability, a decreasing rock strength with sample size using very large core experiments, and a strong temperature dependency of the thermal expansion coefficient for temperatures up to 200°C. The analysis of in situ heater experiments showed that fracture

  10. Primary migration by diffusion through kerogen: I. Model experiments with organic-coated rocks

    SciTech Connect

    Thomas, M.M.; Clouse, J.A. )

    1990-10-01

    Laboratory experiments were performed to assess diffusion through kerogen as a mechanism of hydrocarbon transport through fine-grained rocks. Such transport can be important in primary migration within source rocks and in leakage through seals. To test the concept of diffusion through organic matter networks, model experiments were performed in which hydrocarbon diffusion was measured through Austin chalk cores that had been coated with a monolayer of fatty acids. Hydrocarbon fluxes through the coated cores were compared to hydrocarbon fluxes through uncoated Austin chalk cores. Results showed that the organic coating enhanced transport through the core by a hundredfold over diffusion through its water-filled pore space alone.

  11. Committee neural network model for rock permeability prediction

    NASA Astrophysics Data System (ADS)

    Bagheripour, Parisa

    2014-05-01

    Quantitative formulation between conventional well log data and rock permeability, undoubtedly the most critical parameter of hydrocarbon reservoir, could be a potent tool for solving problems associated with almost all tasks involved in petroleum engineering. The present study proposes a novel approach in charge of the quest for high-accuracy method of permeability prediction. At the first stage, overlapping of conventional well log data (inputs) was eliminated by means of principal component analysis (PCA). Subsequently, rock permeability was predicted from extracted PCs using multi-layer perceptron (MLP), radial basis function (RBF), and generalized regression neural network (GRNN). Eventually, a committee neural network (CNN) was constructed by virtue of genetic algorithm (GA) to enhance the precision of ultimate permeability prediction. The values of rock permeability, derived from the MPL, RBF, and GRNN models, were used as inputs of CNN. The proposed CNN combines results of different ANNs to reap beneficial advantages of all models and consequently producing more accurate estimations. The GA, embedded in the structure of the CNN assigns a weight factor to each ANN which shows relative involvement of each ANN in overall prediction of rock permeability from PCs of conventional well logs. The proposed methodology was applied in Kangan and Dalan Formations, which are the major carbonate reservoir rocks of South Pars Gas Field-Iran. A group of 350 data points was used to establish the CNN model, and a group of 245 data points was employed to assess the reliability of constructed CNN model. Results showed that the CNN method performed better than individual intelligent systems performing alone.

  12. Fractured rock stress-permeability relationships from in situ data and effects of temperature and chemical-mechanical couplings

    SciTech Connect

    Rutqvist, J.

    2014-09-19

    The purpose of this paper is to (i) review field data on stress-induced permeability changes in fractured rock; (ii) describe estimation of fractured rock stress-permeability relationships through model calibration against such field data; and (iii) discuss observations of temperature and chemically mediated fracture closure and its effect on fractured rock permeability. The field data that are reviewed include in situ block experiments, excavation-induced changes in permeability around tunnels, borehole injection experiments, depth (and stress) dependent permeability, and permeability changes associated with a large-scale rock-mass heating experiment. Data show how the stress-permeability relationship of fractured rock very much depends on local in situ conditions, such as fracture shear offset and fracture infilling by mineral precipitation. Field and laboratory experiments involving temperature have shown significant temperature-driven fracture closure even under constant stress. Such temperature-driven fracture closure has been described as thermal overclosure and relates to better fitting of opposing fracture surfaces at high temperatures, or is attributed to chemically mediated fracture closure related to pressure solution (and compaction) of stressed fracture surface asperities. Back-calculated stress-permeability relationships from field data may implicitly account for such effects, but the relative contribution of purely thermal-mechanical and chemically mediated changes is difficult to isolate. Therefore, it is concluded that further laboratory and in situ experiments are needed to increase the knowledge of the true mechanisms behind thermally driven fracture closure, and to further assess the importance of chemical-mechanical coupling for the long-term evolution of fractured rock permeability.

  13. Fractured rock stress-permeability relationships from in situ data and effects of temperature and chemical-mechanical couplings

    DOE PAGESBeta

    Rutqvist, J.

    2014-09-19

    The purpose of this paper is to (i) review field data on stress-induced permeability changes in fractured rock; (ii) describe estimation of fractured rock stress-permeability relationships through model calibration against such field data; and (iii) discuss observations of temperature and chemically mediated fracture closure and its effect on fractured rock permeability. The field data that are reviewed include in situ block experiments, excavation-induced changes in permeability around tunnels, borehole injection experiments, depth (and stress) dependent permeability, and permeability changes associated with a large-scale rock-mass heating experiment. Data show how the stress-permeability relationship of fractured rock very much depends on localmore » in situ conditions, such as fracture shear offset and fracture infilling by mineral precipitation. Field and laboratory experiments involving temperature have shown significant temperature-driven fracture closure even under constant stress. Such temperature-driven fracture closure has been described as thermal overclosure and relates to better fitting of opposing fracture surfaces at high temperatures, or is attributed to chemically mediated fracture closure related to pressure solution (and compaction) of stressed fracture surface asperities. Back-calculated stress-permeability relationships from field data may implicitly account for such effects, but the relative contribution of purely thermal-mechanical and chemically mediated changes is difficult to isolate. Therefore, it is concluded that further laboratory and in situ experiments are needed to increase the knowledge of the true mechanisms behind thermally driven fracture closure, and to further assess the importance of chemical-mechanical coupling for the long-term evolution of fractured rock permeability.« less

  14. Quantifying Biofilm in Porous Media Using Rock Physics Models

    NASA Astrophysics Data System (ADS)

    Alhadhrami, F. M.; Jaiswal, P.; Atekwana, E. A.

    2012-12-01

    Biofilm formation and growth in porous rocks can change their material properties such as porosity, permeability which in turn will impact fluid flow. Finding a non-intrusive method to quantify biofilms and their byproducts in rocks is a key to understanding and modeling bioclogging in porous media. Previous geophysical investigations have documented that seismic techniques are sensitive to biofilm growth. These studies pointed to the fact that microbial growth and biofilm formation induces heterogeneity in the seismic properties. Currently there are no rock physics models to explain these observations and to provide quantitative interpretation of the seismic data. Our objectives are to develop a new class of rock physics model that incorporate microbial processes and their effect on seismic properties. Using the assumption that biofilms can grow within pore-spaces or as a layer coating the mineral grains, P-wave velocity (Vp) and S-wave (Vs) velocity models were constructed using travel-time and waveform tomography technique. We used generic rock physics schematics to represent our rock system numerically. We simulated the arrival times as well as waveforms by treating biofilms either as fluid (filling pore spaces) or as part of matrix (coating sand grains). The preliminary results showed that there is a 1% change in Vp and 3% change in Vs when biofilms are represented discrete structures in pore spaces. On the other hand, a 30% change in Vp and 100% change in Vs was observed when biofilm was represented as part of matrix coating sand grains. Therefore, Vp and Vs changes are more rapid when biofilm grows as grain-coating phase. The significant change in Vs associated with biofilms suggests that shear velocity can be used as a diagnostic tool for imaging zones of bioclogging in the subsurface. The results obtained from this study have significant implications for the study of the rheological properties of biofilms in geological media. Other applications include

  15. Water and CO2 chemistry influences on the mechanical integrity of rocks

    NASA Astrophysics Data System (ADS)

    Darling, T. W.; Le Bas, P.; Carey, J. W.; Johnson, P. A.

    2010-12-01

    The internal bonds in porous rocks expose large surface areas to penetrating fluids which can carry chemical reactants leading to dissolution, accretion, and reactions to new compounds or phases. This may alter the bond mechanical properties to change the bulk dynamical behavior, the overall elasticity and ultimate strength of the rock. Large volume mechanical changes in rock may lead to unexpected seismic activity. We are examining the chemical-mechanical response of porous rocks under a range of controlled environments to learn how vast volumes of porous rock will behave with changing H2O/CO2 concentrations as envisaged under various carbon sequestration schemes. Our program is to measure the linear and non-linear elastic properties of lab-sized samples over the range of UHV (almost no free water or CO2) to supercritical CO2 saturated water at high temperature and pressures. Nonlinear elastic properties are linked to the material integrity while the linear properties are related to the complex modulus and density. We present data on core samples of Berea sandstone, carbonate/arenite mixtures and limestones from dynamical nonlinearity measurements. We find that a carbonate rich sandstone exposed to supercritical CO2 in the presence of water shows significant decrease in both its linear wavespeed and nonlinear elastic response, and there is a marked decrease in mass density. Removal of both water and H2O to partial pressures below 10-8 Torr changes, but does not remove the essential nonlinear behavior of the rock.

  16. A huge deep-seated ancient rock landslide: recognition, mechanism and stability

    NASA Astrophysics Data System (ADS)

    Tang, M. G.; Xu, Q.; Li, Y. S.; Huang, R. Q.; Zheng, G.

    2015-11-01

    The identification of deep-seated landslides is a difficult problem and its failure mechanism is a research hotspot. This paper mainly discusses a very attractive huge deep-seated ancient landslide, it is a very good case to go further research. About 15 years ago a large-scale abnormal geomorphy and geological phenomenon, containing a discontinuous stratum in output and color, was found in the new city of Fengjie, Three Gorges Project Reservoir, China. Two hypotheses for the interpretation of the abnormal phenomenon are a fault graben or a large-scale landslide. From then on continue collecting and analyzing relevant information, field investigation and test, now the results show that the fault graben, consisting of normal faults, could not have been formed under the north-south compressive structure stress of the local region. Meanwhile, a lot of unique geological features, interesting sliding trails and marks of the ancient landslide are discovered and identified in field and experiments. The deformation process and failure mechanism of the ancient landslide are clearly reappeared by a large centrifuge model experiment. Its failure mechanism can be analyzed as "creep-crack-cut". The experiment strongly confirms that it is a huge deep-seated ancient rock landslide. And the failure precursor and key factors of rock slope are discussed. At last, the stability analysis shows that the landslide as a whole is stable and the secondary landslides at the front are basically stable. The results provide a technical support for decision making of the land use planning and construction of the new city, Fengjie.

  17. A rock physics model for analysis of anisotropic parameters in a shale reservoir in Southwest China

    NASA Astrophysics Data System (ADS)

    Qian, Keran; Zhang, Feng; Chen, Shuangquan; Li, Xiangyang; Zhang, Hui

    2016-02-01

    A rock physics model is a very effective tool to describe the anisotropy and mechanical properties of rock from a seismology perspective. Compared to a conventional reservoir, modelling a shale reservoir requires us to face two main challenges in modelling: the existence of organic matter and strong anisotropy. We construct an anisotropic rock physics workflow for a typical shale reservoir in Southwest China, in which the organic matter is treated separately from other minerals by using a combination of anisotropic self-consistent approximation and the differential effective medium method. The standard deviation of the distribution function is used to model the degree of lamination of clay and kerogen. A double scan workflow is introduced to invert the probability of pore aspect ratio and lamination simultaneously, which can give us a better understanding of the shale formation. The anisotropic properties of target formation have been analysed based on the proposed model. Inverted Thomsen parameters, especially the sign of delta, are analysed in terms of the physical properties of rock physics modelling.

  18. Fault damage zones in mechanically layered rocks: The effects of planar anisotropy

    NASA Astrophysics Data System (ADS)

    Misra, Santanu; Ellis, Susan; Mandal, Nibir

    2015-08-01

    This study shows how inherited strength anisotropy influences damage localization at both the tip and wall regions of a fault or fracture. We performed analogue and numerical compression experiments on transversely isotropic models with single and multiple cuts of finite length, simulating the propagation of preexisting faults and cracks in layered rock. The stress-strain curves from the analogue experiments show a change in bulk yield behavior with fault inclination and anisotropy orientation with respect to the stress direction. Earlier isotropic models demonstrated a brittle (wing fracturing) to ductile (shear-zone formation) transition as the fault angle (α) to the principal compression direction increased. The experiments with anisotropic models show patterns of damage localization change dramatically with the orientation of transversely isotropic planes (θ, measured with respect to principal extension direction). Under layer-normal (θ = 0°) and layer-parallel compression (θ = 90°), preexisting faults undergo significant reactivation when 0 < α < 90°, and fault slip eventually leads to mechanical instabilities within the anisotropic layering, causing damage zones in the tip regions. For layer-normal (θ = 0°) compression, the damage processes involve intense extensional shear localization, whereas for layer-parallel compression, contractional shear localization and tensile opening result in characteristic internal shear-band structures. In contrast, for 0 < θ < 90°, the faults undergo little or no reactivation, irrespective of α. In this case, bulk compression leads to an interlayer slip-mediated global deformation. Obliquely anisotropic models thus produce weak or no fault damage zones. We also show that the fault-parallel principal damage localized at the tips can be coupled with transversely oriented, antithetic secondary damage in the wall regions. However, secondary damage develops predominantly when θ = 90°. Field examples of fault damage

  19. Critical review of the state-of-the-art of fracture mechanics with emphasis on layered rocks

    SciTech Connect

    Kuruppu, M.D.; Cheng, K.P.; Edl, J.N. Jr.

    1983-07-01

    Results are presented of a literature survey of over 70 pertinent publications and critical reviews of fracture mechanics emphasizing the fracture behavior of layered rocks. Historical perspective, fracture mechanisms, linear and nonlinear fracture mechanics, energy theories, ductile and brittle fractures, process regions, specific work of fracture, J-integrals, failure theories, static and dynamic fractures, rock fracture mechanics, fracture toughness of layered rocks (e.g., oil shale), experimental and numerical methods are reviewed and discussed. Innovative and promising methods tailored for the fracture mechanics of layered rocks are recommended.

  20. Earthquake-induced collapse mechanism of two types of dangerous rock masses

    NASA Astrophysics Data System (ADS)

    Wang, Wei; Yuan, Wei; Wang, Qizhi; Xue, Kang

    2016-06-01

    As the economy of China develops, an increasing number of key traffic projects have been undertaken in the west of China, where there are high, steep rock slopes. The collapse of dangerous rock masses, especially following a strong earthquake, is one of common geological disasters known in rock slope engineering. Therefore, it is important to study the collapse mechanism of dangerous rock masses induced by an earthquake and the analysis approach of its stability. This study provides a simple and convenient method to determine the collapse mechanisms of two types of dangerous rock masses (i.e. cantilever and upright) associated with the definition and calculation of the safety factor, which is based on the flexure theory of a constant-section beam by combining with the maximum tensile-stress criterion to depict the process of crack propagation caused by seismic waves. The calculation results show that there are critical crack depths in each form of the dangerous rock masses. Once the accumulated depth of the crack growth during an earthquake exceeds the critical depth, the collapse will occur. It is also demonstrated that the crack extension amount of each step is not a constant value, and is closely associated with the current accumulated crack depth. The greater the cumulative crack depth, the more easily the crack propagates. Finally, the validity and applicability of the proposed method are verified through two actual engineering examples.

  1. Integrating rock mechanics issues with repository design through design process principles and methodology

    SciTech Connect

    Bieniawski, Z.T.

    1996-04-01

    A good designer needs not only knowledge for designing (technical know-how that is used to generate alternative design solutions) but also must have knowledge about designing (appropriate principles and systematic methodology to follow). Concepts such as {open_quotes}design for manufacture{close_quotes} or {open_quotes}concurrent engineering{close_quotes} are widely used in the industry. In the field of rock engineering, only limited attention has been paid to the design process because design of structures in rock masses presents unique challenges to the designers as a result of the uncertainties inherent in characterization of geologic media. However, a stage has now been reached where we are be able to sufficiently characterize rock masses for engineering purposes and identify the rock mechanics issues involved but are still lacking engineering design principles and methodology to maximize our design performance. This paper discusses the principles and methodology of the engineering design process directed to integrating site characterization activities with design, construction and performance of an underground repository. Using the latest information from the Yucca Mountain Project on geology, rock mechanics and starter tunnel design, the current lack of integration is pointed out and it is shown how rock mechanics issues can be effectively interwoven with repository design through a systematic design process methodology leading to improved repository performance. In essence, the design process is seen as the use of design principles within an integrating design methodology, leading to innovative problem solving. In particular, a new concept of {open_quotes}Design for Constructibility and Performance{close_quotes} is introduced. This is discussed with respect to ten rock mechanics issues identified for repository design and performance.

  2. The variation of the mechanical properties of rock on spatial scales from the laboratory to outcrop

    NASA Astrophysics Data System (ADS)

    Gage, J.; Wang, H. F.; Fratta, D.; Maclaughlin, M.; Turner, A. L.; GEOX^TM

    2011-12-01

    We have installed a dense array of Fiber Bragg Grating (FBG) strain and temperature sensors on the 4100'-level (1250 m) at the site of the former Homestake gold mine in Lead, SD. The sensor installation site is composed of the Precambrian Poorman formation that contains deformed and metamorphosed Precambrian sediments that is anisotropic including a well-developed foliation, quartz veins, and several joint sets. We have installed nine Micron Optics Inc. OS3600 tube gages. Four of these gages are mounted on the surface of the rock mass and attached to rock bolts that extend 2 m into the rock mass. The other five OS3600 sensors are embedded in drill holes into the rock mass. Additionally, we have developed a new method for measuring in situ strain and temperature in intact rock masses. Fiber optically instrumented rock strain and temperature strips (FROSTS) are 2 m-long strips of 304 stainless steel specially designed to measure temperature and both shortening and elongation in an intact rock mass. FROSTS have FBG strain and temperature sensors mounted on them at 30 cm interval and are grouted into a drill hole in a rock mass. In May 2011, we performed an active loading experiment that consisted of using two hydraulic rams to apply over 200 kN of force to the rock mass. Elastic strain was measured with the fiber optic sensor array. A one-dimensional Boussinesq solution calculates a Young's Modulus of 6.25 GPa for the rock mass. The laboratory-determined values for Young's Modulus in the Poorman formation vary between 49.6 and 94.5 GPa. The difference between the laboratory and field values can be attributed to the closing of fractures and microcracks in the rock mass making the rock mass more compliant than the smaller specimens used for the laboratory experiments. The results of the active loading experiment have implications for the up-scaling of rock mechanical properties between the laboratory and field scales.

  3. Modelling surface roughness and rocks in LRO Diviner observations

    NASA Astrophysics Data System (ADS)

    Williams, J.-P.; Hayne, P. O.; Paige, D. A.

    2012-09-01

    The Diviner Lunar Radiometer Experiment on NASA's Lunar Reconnaissance Orbiter (LRO) observes radiance in 7 infrared spectral channels from which brightness temperatures of the lunar surface are derived. In general, Diviner's surface footprint contains small scale variations in temperature. This anisothermality results in different observed brightness temperatures in Diviner's individual channels. A three-dimensional heat diffusion model is used to explore anisothermality in Diviner observations resulting from surface roughness and rocks at multiple length-scales and illumination conditions.

  4. Nonlinear creep damage constitutive model for soft rocks

    NASA Astrophysics Data System (ADS)

    Liu, H. Z.; Xie, H. Q.; He, J. D.; Xiao, M. L.; Zhuo, L.

    2016-06-01

    In some existing nonlinear creep damage models, it may be less rigorous to directly introduce a damage variable into the creep equation when the damage variable of the viscous component is a function of time or strain. In this paper, we adopt the Kachanov creep damage rate and introduce a damage variable into a rheological differential constitutive equation to derive an analytical integral solution for the creep damage equation of the Bingham model. We also propose a new nonlinear viscous component which reflects nonlinear properties related to the axial stress of soft rock in the steady-state creep stage. Furthermore, we build an improved Nishihara model by using this new component in series with the correctional Nishihara damage model that describes the accelerating creep, and deduce the rheological constitutive relation of the improved model. Based on superposition principle, we obtain the damage creep equation for conditions of both uniaxial and triaxial compression stress, and study the method for determining the model parameters. Finally, this paper presents the laboratory test results performed on mica-quartz schist in parallel with, or vertical to the schistosity direction, and applies the improved Nishihara model to the parameter identification of mica-quartz schist. Using a comparative analysis with test data, results show that the improved model has a superior ability to reflect the creep properties of soft rock in the decelerating creep stage, the steady-state creep stage, and particularly within the accelerating creep stage, in comparison with the traditional Nishihara model.

  5. Laboratory tools to quantify biogenic dissolution of rocks and minerals: a model rock biofilm growing in percolation columns

    NASA Astrophysics Data System (ADS)

    Seiffert, Franz; Bandow, Nicole; Kalbe, Ute; Milke, Ralf; Gorbushina, Anna

    2016-04-01

    Sub-aerial biofilms (SAB) are ubiquitous, self-sufficient microbial ecosystems found on mineral surfaces at all altitudes and latitudes. SABs, which are the principal causes of weathering on exposed terrestrial surfaces, are characterised by patchy growth dominated by associations of algae, cyanobacteria, fungi and heterotrophic bacteria. A recently developed in vitro system to study colonisation of rocks exposed to air included two key SAB participants - the rock-inhabiting ascomycete Knufia petricola (CBS 123872) and the phototrophic cyanobacterium Nostoc punctiforme ATCC29133. Both partners are genetically tractable and we used them here to study weathering of granite, K-feldspar and plagioclase. Small fragments of the various rocks or minerals (1 to 6 mm) were packed into flow-through columns and incubated with 0.1% glucose and 10 µM thiamine-hydrochloride (90 µL.min-1) to compare weathering with and without biofilms. Dissolution of the minerals was followed by: analysing (i) the degradation products in the effluent from the columns via Inductively Coupled Plasma Spectroscopy and (ii) by studying polished sections of the incubated mineral fragment/grains using scanning electron microscopy, transmission electron microscopy and energy dispersive X-ray analyses. K. petricola/N. punctiforme stimulated release of Ca, Na, Mg and Mn. Analyses of the polished sections confirmed depletion of Ca, Na and K near the surface of the fragments. The abrupt decrease in Ca concentration observed in peripheral areas of plagioclase fragments favoured a dissolution-reprecipitation mechanism. Percolation columns in combination with a model biofilm can thus be used to study weathering in closed systems. Columns can easily be filled with different minerals and biofilms, the effluent as well as grains can be collected after long-term exposure under axenic conditions and easily analysed.

  6. Rock penetration : finite element sensitivity and probabilistic modeling analyses.

    SciTech Connect

    Fossum, Arlo Frederick

    2004-08-01

    This report summarizes numerical analyses conducted to assess the relative importance on penetration depth calculations of rock constitutive model physics features representing the presence of microscale flaws such as porosity and networks of microcracks and rock mass structural features. Three-dimensional, nonlinear, transient dynamic finite element penetration simulations are made with a realistic geomaterial constitutive model to determine which features have the most influence on penetration depth calculations. A baseline penetration calculation is made with a representative set of material parameters evaluated from measurements made from laboratory experiments conducted on a familiar sedimentary rock. Then, a sequence of perturbations of various material parameters allows an assessment to be made of the main penetration effects. A cumulative probability distribution function is calculated with the use of an advanced reliability method that makes use of this sensitivity database, probability density functions, and coefficients of variation of the key controlling parameters for penetration depth predictions. Thus the variability of the calculated penetration depth is known as a function of the variability of the input parameters. This simulation modeling capability should impact significantly the tools that are needed to design enhanced penetrator systems, support weapons effects studies, and directly address proposed HDBT defeat scenarios.

  7. Finite element modelling of frictional instability between deformable rocks

    NASA Astrophysics Data System (ADS)

    Xing, H. L.; Makinouchi, A.

    2003-10-01

    Earthquakes are recognized as resulting from a stick-slip frictional instability along faults. Based on the node-to-point contact element strategy (an arbitrarily shaped contact element strategy applied with the static-explicit algorithm for modelling non-linear frictional contact problems proposed by authors), a finite element code for modelling the 3-D non-linear friction contact between deformable bodies has been developed and extended here to analyse the non-linear stick-slip frictional instability between deformable rocks with a rate- and state-dependent friction law. A typical fault bend model is taken as an application example to be analysed here. The variations of the normal contact force, the frictional force, the transition of stick-slip instable state and the related relative slip velocity along the fault between the deformable rocks and the stress evolution in the total bodies during the different stages are investigated, respectively. The calculated results demonstrate the usefulness of this code for simulating the non-linear frictional instability between deformable rocks. Copyright

  8. Influence of Water Content on Mechanical Properties of Rock in Both Saturation and Drying Processes

    NASA Astrophysics Data System (ADS)

    Zhou, Zilong; Cai, Xin; Cao, Wenzhuo; Li, Xibing; Xiong, Cheng

    2016-08-01

    Water content has a pronounced influence on the properties of rock materials, which is responsible for many rock engineering hazards, such as landslides and karst collapse. Meanwhile, water injection is also used for the prevention of some engineering disasters like rock-bursts. To comprehensively investigate the effect of water content on mechanical properties of rocks, laboratory tests were carried out on sandstone specimens with different water contents in both saturation and drying processes. The Nuclear Magnetic Resonance technique was applied to study the water distribution in specimens with variation of water contents. The servo-controlled rock mechanics testing machine and Split Hopkinson Pressure Bar technique were used to conduct both compressive and tensile tests on sandstone specimens with different water contents. From the laboratory tests, reductions of the compressive and tensile strength of sandstone under static and dynamic states in different saturation processes were observed. In the drying process, all of the saturated specimens could basically regain their mechanical properties and recover its strength as in the dry state. However, for partially saturated specimens in the saturation and drying processes, the tensile strength of specimens with the same water content was different, which could be related to different water distributions in specimens.

  9. Rock deformation models and fluid leak-off in hydraulic fracturing

    NASA Astrophysics Data System (ADS)

    Yarushina, Viktoriya M.; Bercovici, David; Oristaglio, Michael L.

    2013-09-01

    Fluid loss into reservoir rocks during hydraulic fracturing is modelled via a poro-elastoplastic pressure diffusion equation in which the total compressibility is a sum of fluid, rock and pore space compressibilities. Inclusion of pore compressibility and porosity-dependent permeability in the model leads to a strong pressure dependence of leak-off (i.e. drainage rate). Dilation of the matrix due to fluid invasion causes higher rates of fluid leak-off. The present model is appropriate for naturally fractured and tight gas reservoirs as well as for soft and poorly consolidated formations whose mechanical behaviour departs from simple elastic laws. Enhancement of the leak-off coefficient by dilation, predicted by the new model, may help explain the low percentage recovery of fracturing fluid (usually between 5 and 50 per cent) in shale gas stimulation by hydraulic fracturing.

  10. Rock Physics Models of Biofilm Growth in Porous Media

    NASA Astrophysics Data System (ADS)

    Jaiswal, P.; alhadhrami, F. M.; Atekwana, E. A.

    2013-12-01

    Recent studies suggest the potential to use acoustic techniques to image biofilm growth in porous media. Nonetheless the interpretation of the seismic response to biofilm growth and development remains speculative because of the lack of quantitative petrophysical models that can relate changes in biofilm saturation to changes in seismic attributes. Here, we report our efforts in developing quantitative rock physics models to biofilm saturation with increasing and decreasing P-wave velocity (VP) and amplitudes recorded in the Davis et al. [2010] physical scale experiment. We adapted rock physics models developed for modeling gas hydrates in unconsolidated sediments. Two distinct growth models, which appear to be a function of pore throat size, are needed to explain the experimental data. First, introduction of biofilm as an additional mineral grain in the sediment matrix (load-bearing mode) is needed to explain the increasing time-lapse VP. Second, introduction of biofilm as part of the pore fluid (pore-filling mode) is required to explain the decreasing time-lapse VP. To explain the time-lapse VP, up to 15% of the pore volume was required to be saturated with biofilm. The recorded seismic amplitudes, which can be expressed as a function of porosity, permeability and grain size, showed a monotonic time-lapse decay except on Day 3 at a few selected locations, where it increased. Since porosity changes are constrained by VP, amplitude increase could be modeled by increasing hydraulic conductivity. Time lapse VP at locations with increasing amplitudes suggest that these locations have a load-bearing growth style. We conclude that permeability can increase by up to 10% at low (~2%) biofilm saturation in load-bearing growth style due to the development of channels within the biofilm structure. Developing a rock physics model for the biofilm growth in general may help create a field guide for interpreting porosity and permeability changes in bioremediation, MEOR and

  11. Hot Dry Rock Geothermal Reservoir Model Development at Los Alamos

    SciTech Connect

    Robinson, Bruce A.; Birdsell, Stephen A.

    1989-03-21

    Discrete fracture and continuum models are being developed to simulate Hot Dry Rock (HDR) geothermal reservoirs. The discrete fracture model is a two-dimensional steady state simulator of fluid flow and tracer transport in a fracture network which is generated from assumed statistical properties of the fractures. The model's strength lies in its ability to compute the steady state pressure drop and tracer response in a realistic network of interconnected fractures. The continuum approach models fracture behavior by treating permeability and porosity as functions of temperature and effective stress. With this model it is practical to model transient behavior as well as the coupled processes of fluid flow, heat transfer, and stress effects in a three-dimensional system. The model capabilities being developed will also have applications in conventional geothermal systems undergoing reinjection and in fractured geothermal reservoirs in general.

  12. Hot Dry Rock geothermal reservoir model development at Los Alamos

    SciTech Connect

    Robinson, B.A.; Birdsell, S.A.

    1989-01-01

    Discrete fracture and continuum models are being developed to simulate Hot Dry Rock (HDR) geothermal reservoirs. The discrete fracture model is a two-dimensional steady state simulator of fluid flow and tracer transport in a fracture network which is generated from assumed statistical properties of the fractures. The model's strength lies in its ability to compute the steady state pressure drop and tracer response in a realistic network of interconnected fractures. The continuum approach models fracture behavior by treating permeability and porosity as functions of temperature and effective stress. With this model it is practical to model transient behavior as well as the coupled processes of fluid flow, heat transfer, and stress effects in a three-dimensional system. The model capabilities being developed will also have applications in conventional geothermal systems undergoing reinjection and in fractured geothermal reservoirs in general. 15 refs., 7 figs.

  13. A damage mechanics approach for quantifying stress changes due to brittle failure of porous rocks

    NASA Astrophysics Data System (ADS)

    Jacquey, Antoine B.; Cacace, Mauro; Blöcher, Guido; Milsch, Harald; Scheck-Wenderoth, Magdalena

    2016-04-01

    Natural fault zones or man-made injection or production of fluid impact the regional stress distribution in Earth's crust and can be responsible for localized stress discontinuities. Understanding the processes controlling fracturing of the porous rocks and mechanical behaviour of fault zones is therefore of interest for several applications including geothermal energy production. In this contribution, we will present a thermodynamically consistent visco-poroelastic damage model which can deal with the multi-scale and multi-physics nature of the physical processes controlling the deformation of porous rocks during and after brittle failure. Deformation of a porous medium is crucially influenced by the changes in the effective stress. Considering a strain-formulated yield cap and the compaction-dilation transition, three different regimes can be identified: quasi-elastic deformation, cataclastic compaction with microcracking (damage accumulation) and macroscopic brittle failure with dilation. The governing equations for deformation, damage accumulation/healing and fluid flow have been implemented in a fully-coupled finite-element-method based framework (MOOSE). The MOOSE framework provides a powerful and flexible platform to solve multiphysics problems implicitly and in a tightly coupled manner on unstructured meshes which is of interest for such non-linear context. To illustrate the model, simulation of a compaction experiment of a sandstone leading to shear failure will be presented which allows to quantify the stress drop accompanying the failure. Finally, we will demonstrate that this approach can also be used at the field scale to simulate hydraulic fracturing and assess the resulting changes in the stress field.

  14. Modelling Cochlear Mechanics

    PubMed Central

    Elliott, Stephen J.; Teal, Paul D.

    2014-01-01

    The cochlea plays a crucial role in mammal hearing. The basic function of the cochlea is to map sounds of different frequencies onto corresponding characteristic positions on the basilar membrane (BM). Sounds enter the fluid-filled cochlea and cause deflection of the BM due to pressure differences between the cochlear fluid chambers. These deflections travel along the cochlea, increasing in amplitude, until a frequency-dependent characteristic position and then decay away rapidly. The hair cells can detect these deflections and encode them as neural signals. Modelling the mechanics of the cochlea is of help in interpreting experimental observations and also can provide predictions of the results of experiments that cannot currently be performed due to technical limitations. This paper focuses on reviewing the numerical modelling of the mechanical and electrical processes in the cochlea, which include fluid coupling, micromechanics, the cochlear amplifier, nonlinearity, and electrical coupling. PMID:25136555

  15. Modelling of nuclear explosions in hard rock sites

    SciTech Connect

    Brunish, W.M.; App, F.N.

    1993-01-01

    This study represents part of a larger effort to systematically model the effects of differing source region properties on ground motion from underground nuclear explosions at the Nevada Test Site. In previous work by the authors the primary emphasis was on alluvium and both saturated and unsaturated tuff. We have attempted to model events on Pahute Mesa, where either the working point medium, or some of the layers above the working point, or both, are hard rock. The complex layering at these sites, however, has prevented us from drawing unambiguous conclusions about modelling hard rock. In order to learn more about the response of hard rock to underground nuclear explosions, we have attempted to model the PILEDRIVER event. PILEDRIVER was fired on June 2, 1966 in the granite stock of Area 15 at the Nevada Test Site. The working point was at a depth of 462.7 m and the yield was determined to be 61 kt. Numerous surface, sub-surface and free-field measurements were made and analyzed by SRI. An attempt was made to determine the contribution of spall to the teleseismic signal, but proved unsuccessful because most of the data from below-shot-level gauges was lost. Nonetheless, there is quite a bit of good quality data from a variety of locations. We have been able to obtain relatively good agreement with the experimental PILEDRIVER waveforms. In order to do so, we had to model the granodiorite as being considerably weaker than good quality'' granite, and it had to undergo considerable weakening due to shock damage as well. In addition, the near-surface layers had to be modeled as being weak and compressible and as have a much lower sound speed than the material at depth. The is consistent with a fractured and jointed material at depth, and a weathered material near the surface.

  16. Modelling of nuclear explosions in hard rock sites

    SciTech Connect

    Brunish, W.M.; App, F.N.

    1993-05-01

    This study represents part of a larger effort to systematically model the effects of differing source region properties on ground motion from underground nuclear explosions at the Nevada Test Site. In previous work by the authors the primary emphasis was on alluvium and both saturated and unsaturated tuff. We have attempted to model events on Pahute Mesa, where either the working point medium, or some of the layers above the working point, or both, are hard rock. The complex layering at these sites, however, has prevented us from drawing unambiguous conclusions about modelling hard rock. In order to learn more about the response of hard rock to underground nuclear explosions, we have attempted to model the PILEDRIVER event. PILEDRIVER was fired on June 2, 1966 in the granite stock of Area 15 at the Nevada Test Site. The working point was at a depth of 462.7 m and the yield was determined to be 61 kt. Numerous surface, sub-surface and free-field measurements were made and analyzed by SRI. An attempt was made to determine the contribution of spall to the teleseismic signal, but proved unsuccessful because most of the data from below-shot-level gauges was lost. Nonetheless, there is quite a bit of good quality data from a variety of locations. We have been able to obtain relatively good agreement with the experimental PILEDRIVER waveforms. In order to do so, we had to model the granodiorite as being considerably weaker than ``good quality`` granite, and it had to undergo considerable weakening due to shock damage as well. In addition, the near-surface layers had to be modeled as being weak and compressible and as have a much lower sound speed than the material at depth. The is consistent with a fractured and jointed material at depth, and a weathered material near the surface.

  17. ACOUSTICAL IMAGING AND MECHANICAL PROPERTIES OF SOFT ROCK AND MARINE SEDIMENTS

    SciTech Connect

    Thurman E. Scott, Jr.; Younane Abousleiman

    2004-04-01

    The research during this project has concentrated on developing a correlation between rock deformation mechanisms and their acoustic velocity signature. This has included investigating: (1) the acoustic signature of drained and undrained unconsolidated sands, (2) the acoustic emission signature of deforming high porosity rocks (in comparison to their low porosity high strength counterparts), (3) the effects of deformation on anisotropic elastic and poroelastic moduli, and (4) the acoustic tomographic imaging of damage development in rocks. Each of these four areas involve triaxial experimental testing of weak porous rocks or unconsolidated sand and involves measuring acoustic properties. The research is directed at determining the seismic velocity signature of damaged rocks so that 3-D or 4-D seismic imaging can be utilized to image rock damage. These four areas of study are described in the report: (1) Triaxial compression experiments have been conducted on unconsolidated Oil Creek sand at high confining pressures. (2) Initial experiments on measuring the acoustic emission activity from deforming high porosity Danian chalk were accomplished and these indicate that the AE activity was of a very low amplitude. (3) A series of triaxial compression experiments were conducted to investigate the effects of induced stress on the anisotropy developed in dynamic elastic and poroelastic parameters in rocks. (4) Tomographic acoustic imaging was utilized to image the internal damage in a deforming porous limestone sample. Results indicate that the deformation damage in rocks induced during laboratory experimentation can be imaged tomographically in the laboratory. By extension the results also indicate that 4-D seismic imaging of a reservoir may become a powerful tool for imaging reservoir deformation (including imaging compaction and subsidence) and for imaging zones where drilling operation may encounter hazardous shallow water flows.

  18. Modelling fault surface roughness and fault rocks thickness evolution with slip: calibration based on field and laboratory data

    NASA Astrophysics Data System (ADS)

    Bistacchi, A.; Tisato, N.; Spagnuolo, E.; Nielsen, S. B.; Di Toro, G.

    2012-12-01

    deformation processes (e.g. frictional melting vs. cataclasis) and experimental conditions (unconfined vs. confined). Since the model is based on geometrical and volume-conservation considerations (and not on a particular deformation mechanism), we conclude that the surface roughness and fault-rock thickness after some slip is mostly determined by the initial roughness (measured over several orders of magnitude in wavelength), rather than the particular deformation process (cataclasis, melting, etc.) activated during faulting. Conveniently, since the model can be applied (under certain conditions) to surfaces which depart from self-affine roughness, the model parameters can be calibrated with laboratory experiments. If this conclusion will be confirmed by a larger dataset, the forward model proposed here will provide realistic fault roughness and fault rock thickness predictions to be used in the mechanics of earthquakes and faulting, oil and water exploration, and underground engineering projects.

  19. REDBACK: an Open-Source Highly Scalable Simulation Tool for Rock Mechanics with Dissipative Feedbacks

    NASA Astrophysics Data System (ADS)

    Poulet, T.; Veveakis, M.; Paesold, M.; Regenauer-Lieb, K.

    2014-12-01

    Multiphysics modelling has become an indispensable tool for geoscientists to simulate the complex behaviours observed in their various fields of study where multiple processes are involved, including thermal, hydraulic, mechanical and chemical (THMC) laws. This modelling activity involves simulations that are computationally expensive and its soaring uptake is tightly linked to the increasing availability of supercomputing power and easy access to powerful nonlinear solvers such as PETSc (http://www.mcs.anl.gov/petsc/). The Multiphysics Object-Oriented Simulation Environment (MOOSE) is a finite-element, multiphysics framework (http://mooseframework.org) that can harness such computational power and allow scientists to develop easily some tightly-coupled fully implicit multiphysics simulations that run automatically in parallel on large clusters. This open-source framework provides a powerful tool to collaborate on numerical modelling activities and we are contributing to its development with REDBACK (https://github.com/pou036/redback), a module for Rock mEchanics with Dissipative feedBACKs. REDBACK builds on the tensor mechanics finite strain implementation available in MOOSE to provide a THMC simulator where the energetic formulation highlights the importance of all dissipative terms in the coupled system of equations. We show first applications of fully coupled dehydration reactions triggering episodic fluid transfer through shear zones (Alevizos et al, 2014). The dimensionless approach used allows focusing on the critical underlying variables which are driving the resulting behaviours observed and this tool is specifically designed to study material instabilities underpinning geological features like faulting, folding, boudinage, shearing, fracturing, etc. REDBACK provides a collaborative and educational tool which captures the physical and mathematical understanding of such material instabilities and provides an easy way to apply this knowledge to realistic

  20. Discrete element modeling of rock deformation, fracture network development and permeability evolution under hydraulic stimulation

    SciTech Connect

    Shouchun Deng; Robert Podgorney; Hai Huang

    2011-02-01

    Key challenges associated with the EGS reservoir development include the ability to reliably predict hydraulic fracturing and the deformation of natural fractures as well as estimating permeability evolution of the fracture network with time. We have developed a physics-based rock deformation and fracture propagation simulator by coupling a discrete element model (DEM) for fracturing with a network flow model. In DEM model, solid rock is represented by a network of discrete elements (often referred as particles) connected by various types of mechanical bonds such as springs, elastic beams or bonds that have more complex properties (such as stress-dependent elastic constants). Fracturing is represented explicitly as broken bonds (microcracks), which form and coalesce into macroscopic fractures when external and internal load is applied. The natural fractures are represented by a series of connected line segments. Mechanical bonds that intersect with such line segments are removed from the DEM model. A network flow model using conjugate lattice to the DEM network is developed and coupled with the DEM. The fluid pressure gradient exerts forces on individual elements of the DEM network, which therefore deforms the mechanical bonds and breaks them if the deformation reaches a prescribed threshold value. Such deformation/fracturing in turn changes the permeability of the flow network, which again changes the evolution of fluid pressure, intimately coupling the two processes. The intimate coupling between fracturing/deformation of fracture networks and fluid flow makes the meso-scale DEM- network flow simulations necessary in order to accurately evaluate the permeability evolution, as these methods have substantial advantages over conventional continuum mechanical models of elastic rock deformation. The challenges that must be overcome to simulate EGS reservoir stimulation, preliminary results, progress to date and near future research directions and opportunities will be

  1. Stress-dependent voltage offsets from polymer insulators used in rock mechanics and material testing

    NASA Astrophysics Data System (ADS)

    Carlson, G. G.; Dahlgren, R.; Vanderbilt, V. C.; Johnston, M. J.; Dunson, C.; Gray, A.; Freund, F.

    2013-12-01

    Dielectric insulators are used in a variety of laboratory settings when performing experiments in rock mechanics, petrology, and electromagnetic studies of rocks in the fields of geophysics, material science, and civil engineering. These components may be used to electrically isolate geological samples from the experimental equipment, to perform a mechanical compliance function between brittle samples and the loading equipment, to match ultrasonic transducers, or perform other functions. In many experimental configurations the insulators bear the full brunt of force applied to the sample but do not need to withstand high voltages, therefore the insulators are often thin sheets of mechanically tough polymers. From an instrument perspective, transduction from various types of mechanical perturbation has been qualitatively compared for a number of polymers [1, 2] and these error sources are readily apparent during high-impedance measurements if not mitigated. However even when following best practices, a force-dependent voltage signal still remains and its behavior is explored in this presentation. In this experiment two thin sheets (0.25 mm) of high-density polyethylene (HDPE) were set up in a stack, held alternately between three aluminum bars; this stack was placed on the platen of a 60T capacity hydraulic testing machine. The surface area, A, over which the force is applied to the PE sheets in this sandwich is roughly 40 square cm, each sheet forming a parallel-plate capacitor having roughly 320 pF [3], assuming the relative dielectric permittivity of PE is ~2.3. The outer two aluminum bars were connected to the LO input of the electrometer and the central aluminum bar was connected to the HI input of a Keithley model 617 electrometer. Once the stack is mechanically well-seated with no air gaps, the voltage offset is observed to be a linear function of the baseline voltage for a given change in applied force. For a periodically applied force of 66.7 kN the voltage

  2. VNIR spectral modeling of Mars analogue rocks: first results

    NASA Astrophysics Data System (ADS)

    Pompilio, L.; Roush, T.; Pedrazzi, G.; Sgavetti, M.

    Knowledge regarding the surface composition of Mars and other bodies of the inner solar system is fundamental to understanding of their origin, evolution, and internal structures. Technological improvements of remote sensors and associated implications for planetary studies have encouraged increased laboratory and field spectroscopy research to model the spectral behavior of terrestrial analogues for planetary surfaces. This approach has proven useful during Martian surface and orbital missions, and petrologic studies of Martian SNC meteorites. Thermal emission data were used to suggest two lithologies occurring on Mars surface: basalt with abundant plagioclase and clinopyroxene and andesite, dominated by plagioclase and volcanic glass [1,2]. Weathered basalt has been suggested as an alternative to the andesite interpretation [3,4]. Orbital VNIR spectral imaging data also suggest the crust is dominantly basaltic, chiefly feldspar and pyroxene [5,6]. A few outcrops of ancient crust have higher concentrations of olivine and low-Ca pyroxene, and have been interpreted as cumulates [6]. Based upon these orbital observations future lander/rover missions can be expected to encounter particulate soils, rocks, and rock outcrops. Approaches to qualitative and quantitative analysis of remotely-acquired spectra have been successfully used to infer the presence and abundance of minerals and to discover compositionally associated spectral trends [7-9]. Both empirical [10] and mathematical [e.g. 11-13] methods have been applied, typically with full compositional knowledge, to chiefly particulate samples and as a result cannot be considered as objective techniques for predicting the compositional information, especially for understanding the spectral behavior of rocks. Extending the compositional modeling efforts to include more rocks and developing objective criteria in the modeling are the next required steps. This is the focus of the present investigation. We present results of

  3. Modelling of rock temperatures for deep alpine tunnel projects

    NASA Astrophysics Data System (ADS)

    Goy, L.; Fabre, D.; Menard, G.

    1996-01-01

    The construction of deep railway tunnels requires the prediction of natural temperatures at depth. Geothermal data for the Alps are presented and principles of previously employed methods to predict temperatures, using Andreae's analytical approach, are discussed. We then use a finite element numerical model based on pure conduction to calculate temperatures at depth. This method allows rock heterogeneity and anisotropy to be taken into account. This model is applied to the Maurienne-Ambin tunnel project, a 55 km long tunnel between St-Jean-de-Maurienne (France) and Susa (Italy), which will be the longest tunnel for the planned TGV (high speed train) Lyon-Torino link. Data from several deep boreholes (10 total, with 3>1000 m) are used to provide essential parameters for the model, i.e.: - geological structure; - geothermal gradients; - rock conductivities from cores; - geothermal deep heat flow. Modelling is done in two dimensions, but the effect of surface topography (3 D) is considered. Results are given in the form of a geothermal cross-section along the tunnel axis that provides maximum temperatures and lengths of zones of high temperature encountered (for instance, zones where θ is ≥40°C). In general, differences between calculated and measured temperatures are less than 1°C at great depth. At shallow depth, differences are sometimes higher and probably best explained by water circulation connected to the surface. The modelling of temperatures, in relation to the geological structure, rock properties, and geothermal data for this area, appears to be a very useful tool for comparing alternative routes for deep tunnel projects and, during construction, to predict potential local geological or hydrological anomalies.

  4. Albite [yields] jadeite + quartz transformation in rock: Mechanism and kinetics

    SciTech Connect

    Bohlen, S.R.; Kirby, S.H. ); Hacker, B.R.

    1992-01-01

    Recent work on the calcite [yields] aragonite transformation using fully dense marble revealed significant differences from earlier experiments on powders and single-crystals. The reaction rate is retarded by a factor of > 1,000 and reaction mechanisms and resultant textures are considerably more complex. Stimulated by this, the authors conducted a study of the albite [yields] jadeite + quartz/coesite transformation in a fully dense albitite. Again the results are in marked contrast with previous powder-based studies of this archetypal metamorphic reaction. Solid cores of albitite were held at temperatures of 500-1,200 C and at pressure oversteps of 500 MPa into the jadeite + quartz stability field for 1--8 days in piston-cylinder apparatus. Samples that were dried in vacuum transformed appreciably only at temperatures in excess of 1,000 C. At all grain boundaries there is subequal transformation to micron-scale intergrowths of jadeite + quartz. Samples that were vacuum-impregnated with 1 wt% water contain jadeite + quartz to temperatures as low as 600 C. In contrast to the dried samples, transformation is much less homogeneous. The jadeite + quartz intergrowths do not form rows of subparallel crystals on grain boundaries, but rather are flower-shaped clusters that radiate outward from single nucleation sites at 3-grain edges and 4-grain corners. Compared to powders, pressure oversteps a factor of 10 greater are necessary to induce equivalent reaction in albitite. The sluggishness of this reaction has important implications for the evolution of the lower continental crust and subducting oceanic crust in terms of their (1) seismic velocity profiles, (2) petrological evolution, and (3) buoyancy forces, stresses and vertical crustal movements connected with densification and dilatational reactions.

  5. Thermo-mechanical pressurization of experimental faults in cohesive rocks during seismic slip

    NASA Astrophysics Data System (ADS)

    Violay, M.; Di Toro, G.; Nielsen, S.; Spagnuolo, E.; Burg, J. P.

    2015-11-01

    Earthquakes occur because fault friction weakens with increasing slip and slip rates. Since the slipping zones of faults are often fluid-saturated, thermo-mechanical pressurization of pore fluids has been invoked as a mechanism responsible for frictional dynamic weakening, but experimental evidence is lacking. We performed friction experiments (normal stress 25 MPa, maximal slip-rate ∼3 ms-1) on cohesive basalt and marble under (1) room-humidity and (2) immersed in liquid water (drained and undrained) conditions. In both rock types and independently of the presence of fluids, up to 80% of frictional weakening was measured in the first 5 cm of slip. Modest pressurization-related weakening appears only at later stages of slip. Thermo-mechanical pressurization weakening of cohesive rocks can be negligible during earthquakes due to the triggering of more efficient fault lubrication mechanisms (flash heating, frictional melting, etc.).

  6. A unified asperity-deformation model for cracked rocks

    NASA Astrophysics Data System (ADS)

    Gao, K.; Gibson, R. L.; Ge, J.

    2010-12-01

    Seismic velocities in rocks increase with pressure, a pattern often explained by changes in cracks within the rock volume. Specifically, the width of microcracks will decrease, increasing the contact area of crack surfaces, which in turn leads to increased stiffness of the crack. This causes P- and S-wave velocities to increase, and having accurate models of this behavior is important in many applications where it is important to related stress changes and observed seismic velocities, such as investigations of both fault zones and changing conditions in geothermal fields or hydrocarbon reservoirs. Several different models have been proposed to explain and quantify the relationship between confining pressure and the physical properties of cracked rocks, including a number of solutions describing cracks as ellipsoidal voids (penny-shaped cracks) of with varying aspect ratios that close at different pressures. Differential effective medium theories based on this model will typically use a fairly large number of parameters to fit measured data. An alternative approach describes cracks or fractures as rough surfaces that do not fully close with increasing pressure, and estimates the increase in contact area as asperities come into contact. An asperity-deformation model of this type was developed by, e.g., Gangi and Carlson (1985), and it parameterizes asperities on crack surfaces in terms of a set of cylindrical rods with heights following power-law distribution. The model results in a simple expression for the increase in velocity caused by increased number of asperities (rods) in contact with increasing pressure, an expression that accurately models increases in velocity using only three parameters. However, the model was formulated only for normal modulus of the cracks, the modulus relevant for P-wave propagation across the crack. Since S-waves are controlled by the tangential or shear modulus of the crack, this model cannot be used to jointly invert both P- and S

  7. Stress-Dependent Voltage Offsets From Polymer Insulators Used in Rock Mechanics and Material Testing

    NASA Technical Reports Server (NTRS)

    Carlson, G. G.; Dahlgren, Robert; Gray, Amber; Vanderbilt, V. C.; Freund, F.; Johnston, M. J.; Dunson, C.

    2013-01-01

    Dielectric insulators are used in a variety of laboratory settings when performing experiments in rock mechanics, petrology, and electromagnetic studies of rocks in the fields of geophysics,material science, and civil engineering. These components may be used to electrically isolate geological samples from the experimental equipment, to perform a mechanical compliance function between brittle samples and the loading equipment, to match ultrasonic transducers, or perform other functions. In manyexperimental configurations the insulators bear the full brunt of force applied to the sample but do not need to withstand high voltages, therefore the insulators are often thin sheets of mechanically tough polymers. From an instrument perspective, transduction from various types of mechanical perturbation has beenqualitatively compared for a number of polymers [1, 2] and these error sources are readily apparent duringhigh-impedance measurements if not mitigated. However even when following best practices, a force dependent voltage signal still remains and its behavior is explored in this presentation. In this experimenttwo thin sheets (0.25 mm) of high-density polyethylene (HDPE) were set up in a stack, held alternatelybetween three aluminum bars; this stack was placed on the platen of a 60T capacity hydraulic testingmachine. The surface area, A, over which the force is applied to the PE sheets in this sandwich is roughly 40 square cm, each sheet forming a parallel-plate capacitor having roughly 320 pF [3], assuming therelative dielectric permittivity of PE is approximately 2.3. The outer two aluminum bars were connected to the LO input ofthe electrometer and the central aluminum bar was connected to the HI input of a Keithley model 617 electrometer. Once the stack is mechanically well-seated with no air gaps, the voltage offset is observed tobe a linear function of the baseline voltage for a given change in applied force. For a periodically appliedforce of 66.7 kN the

  8. Thermo-osmosis coupled-flow characterization in clay-rocks: experiments and modeling

    NASA Astrophysics Data System (ADS)

    Tremosa, J.; Goncalves, J.; Matray, J.; Violette, S.

    2009-12-01

    Water flow in clay-rocks is not only driven by a hydraulic gradient but also by chemical, thermal or electrical gradients. It implies a re-evaluation of the Darcy law by considering all gradients occurring in the clay-rock and their associated coupling coefficients (e.g. the osmotic efficiency to link a chemical gradient to a water flow). The occurrence of such processes in clay-rocks is due to the low hydraulic conductivity of this media and because of electrical charges at the clay minerals surface. Here, we focused on the thermo-osmosis process, a water flow under a temperature gradient, which is poorly characterized in spite of its implications in nuclear waste storage in clay-rocks. A set of thermo-osmotic experiments was performed in an equipped borehole installed in a Toarcian compacted clay at the IRSN’s Underground Research Laboratory in the south of France. The water flow induced by a temperature gradient (from the hotter towards the colder zone) was reproduced by the help of a numerical model, including coupled-flow processes, mass conservation laws and hydro-thermo-mechanical changes (see Figure). A range of thermo-osmotic permeability (kT), between 6.10-12 and 2.10-10 m2.K-1.s-2, was obtained during the experiments depending on the temperature gradient and uncertainties on the model parameters. Values obtained for the Tournemire’s argillite are in the high range of thermo-osmotic permeabilities for argillaceous materials and suggest an effect of pore size on the thermo-osmotic permeability of a clay-rock (kT being higher with little pore size). Another dependence of thermo-osmotic permeability with temperature is observed, with kT decreasing when the temperature increases. These experiments and modeling indicate thermo-osmosis will have an influence on water flow in presence of a temperature gradient and this process is to consider in water flow studies in clay-rocks. Reference: Tremosa et al. Estimating thermo-osmotic coefficients in clay-rocks

  9. Leaching of boron, arsenic and selenium from sedimentary rocks: II. pH dependence, speciation and mechanisms of release.

    PubMed

    Tabelin, Carlito Baltazar; Hashimoto, Ayaka; Igarashi, Toshifumi; Yoneda, Tetsuro

    2014-03-01

    Sedimentary rocks excavated in Japan from road- and railway-tunnel projects contain relatively low concentrations of hazardous trace elements like boron (B), arsenic (As) and selenium (Se). However, these seemingly harmless waste rocks often produced leachates with concentrations of hazardous trace elements that exceeded the environmental standards. In this study, the leaching behaviors and release mechanisms of B, As and Se were evaluated using batch leaching experiments, sequential extraction and geochemical modeling calculations. The results showed that B was mostly partitioned with the residual/crystalline phase that is relatively stable under normal environmental conditions. In contrast, the majority of As and Se were associated with the exchangeable and organics/sulfides phases that are unstable under oxidizing conditions. Dissolution of water-soluble phases controlled the leaching of B, As and Se from these rocks in the short term, but pyrite oxidation, calcite dissolution and adsorption/desorption reactions became more important in the long term. The mobilities of these trace elements were also strongly influenced by the pH of the rock-water system. Although the leaching of Se only increased in the acidic region, those of B and As were enhanced under both acidic and alkaline conditions. Under strongly acidic conditions, the primarily release mechanism of B, As and Se was the dissolution of mineral phases that incorporated and/or adsorbed these elements. Lower concentrations of these trace elements in the circumneutral pH range could be attributed to their strong adsorption onto minerals like Al-/Fe-oxyhydroxides and clays, which are inherently present and/or precipitated in the rock-water system. The leaching of As and B increased under strongly alkaline conditions because of enhanced desorption and pyrite oxidation while that of Se remained minimal due to its adsorption onto Fe-oxyhydroxides and co-precipitation with calcite. PMID:24370699

  10. Three-Dimensional Numerical Investigations of the Failure Mechanism of a Rock Disc with a Central or Eccentric Hole

    NASA Astrophysics Data System (ADS)

    Wang, S. Y.; Sloan, S. W.; Tang, C. A.

    2014-11-01

    The diametrical compression of a circular disc (Brazilian test) or cylinder with a small eccentric hole is a simple but important test to determine the tensile strength of rocks. This paper studies the failure mechanism of circular disc with an eccentric hole by a 3D numerical model (RFPA3D). A feature of the code RFPA3D is that it can numerically simulate the evolution of cracks in three-dimensional space, as well as the heterogeneity of the rock mass. First, numerically simulated Brazilian tests are compared with experimental results. Special attention is given to the effect of the thickness to radius ratio on the failure modes and the peak stress of specimens. The effects of the compressive strength to tensile strength ratio ( C/T), the loading arc angle (2 α), and the homogeneity index ( m) are also studied in the numerical simulations. Secondly, the failure process of a rock disc with a central hole is studied. The effects of the ratio of the internal hole radius ( r) to the radius of the rock disc ( R) on the failure mode and the peak stress are investigated. Thirdly, the influence of the vertical and horizontal eccentricity of an internal hole on the initiation and propagation of cracks inside a specimen are simulated. The effect of the radius of the eccentric hole and the homogeneity index ( m) are also investigated.

  11. Subcritical fracture propagation in rocks: An examination using the methods of fracture mechanics and non-destructive testing. Ph.D. Thesis

    NASA Technical Reports Server (NTRS)

    Swanson, P. L.

    1984-01-01

    An experimental investigation of tensile rock fracture is presented with an emphasis on characterizing time dependent crack growth using the methods of fracture mechanics. Subcritical fracture experiments were performed in moist air on glass and five different rock types at crack velocities using the double torsion technique. The experimental results suggest that subcritical fracture resistance in polycrystals is dominated by microstructural effects. Evidence for gross violations of the assumptions of linear elastic fracture mechanics and double torsion theory was found in the tests on rocks. In an effort to obtain a better understanding of the physical breakdown processes associated with rock fracture, a series of nondestructive evaluation tests were performed during subcritical fracture experiments on glass and granite. Comparison of the observed process zone shape with that expected on the basis of a critical normal principal tensile stress criterion shows that the zone is much more elongated in the crack propagation direction than predicted by the continuum based microcracking model alone.

  12. Hydromechanical coupling in fractured rock masses: mechanisms and processes of selected case studies

    NASA Astrophysics Data System (ADS)

    Zangerl, Christian

    2015-04-01

    in the range of millimetres to a very few centimetres and can be linked to annual groundwater fluctuations. Due to pore pressure variations HM coupling can influence seepage forces and effective stresses in the rock mass. Effective stress changes can adversely affect the stability and deformation behaviour of deep-seated rock slides by influencing the shear strength or the time dependent (viscous) material behaviour of the basal shear zone. The shear strength of active shear zones is often reasonably well described by Coulomb's law. In Coulomb's law the operative normal stresses to the shear surface/zone are effective stresses and hence pore pressures which should be taken into account reduces the shear strength. According to the time dependent material behaviour a few effective stress based viscous models exists which are able to consider pore pressures. For slowly moving rock slides HM coupling could be highly relevant when low-permeability clayey-silty shear zones (fault gouges) are existing. An important parameters therefore is the hydraulic diffusivity, which is controlled by the permeability and fluid-pore compressibility of the shear zone, and by fluid viscosity. Thus time dependent pore pressure diffusion in the shear zone can either control the stability condition or the viscous behaviour (creep) of the rock slide. Numerous cases studies show that HM coupling can effect deformability, shear strength and time dependent behaviour of fractured rock masses. A process-based consideration can be important to avoid unexpected impacts on infrastructures and to understand complex rock mass as well rock slide behaviour.

  13. Rock physics models for constraining quantitative interpretation of ultrasonic data for biofilm growth and development

    NASA Astrophysics Data System (ADS)

    Alhadhrami, Fathiya Mohammed

    This study examines the use of rock physics modeling for quantitative interpretation of seismic data in the context of microbial growth and biofilm formation in unconsolidated sediment. The impetus for this research comes from geophysical experiments by Davis et al. (2010) and Kwon and Ajo-Franklin et al. (2012). These studies observed that microbial growth has a small effect on P-wave velocities (VP) but a large effect on seismic amplitudes. Davis et al. (2010) and Kwon and Ajo-Franklin et al. (2012) speculated that the amplitude variations were due to a combination of rock mechanical changes from accumulation of microbial growth related features such as biofilms. A more definite conclusion can be drawn by developing rock physics models that connect rock properties to seismic amplitudes. The primary objective of this work is to provide an explanation for high amplitude attenuation due to biofilm growth. The results suggest that biofilm formation in the Davis et al. (2010) experiment exhibit two growth styles: a loadbearing style where biofilm behaves like an additional mineral grain and a non-loadbearing mode where the biofilm grows into the pore spaces. In the loadbearing mode, the biofilms contribute to the stiffness of the sediments. We refer to this style as "filler." In the non-loadbearing mode, the biofilms contribute only to change in density of sediments without affecting their strength. We refer to this style of microbial growth as "mushroom." Both growth styles appear to be changing permeability more than the moduli or the density. As the result, while the VP velocity remains relatively unchanged, the amplitudes can change significantly depending on biofilm saturation. Interpreting seismic data from biofilm growths in term of rock physics models provide a greater insight into the sediment-fluid interaction. The models in turn can be used to understand microbial enhanced oil recovery and in assisting in solving environmental issues such as creating bio

  14. Inner crack reconstruction and mechanical analysis for rock-specimen-based phase measuring profilometry

    NASA Astrophysics Data System (ADS)

    Cao, Yiping; He, Yuhang

    2009-12-01

    A higher precise inner crack three-dimenssional(3D) reconstructed method of rock specimens is presented. Two inner shapes of the crack are measured with Phase Measuring Profilometry (PMP), and their edges are drawn out by height information instead of the traditional method based on gray information. Subsequently contour matching and height matching are carried out with algorithms due to maximum correlativity. The inner width and volume of the crack are educed according to the fissure of a rock specimen's outer surface, and the 3D profile of the crack is reconstructed with a high repetitive precision superior to 20μm. The proposed method is effective for evaluating the crack's width of rock specimens in the exploitation of petroleum and natural gas with a mechanical analysis method. The experiment shows its feasibility and practicability.

  15. Inner crack reconstruction and mechanical analysis for rock-specimen-based phase measuring profilometry

    NASA Astrophysics Data System (ADS)

    Cao, Yiping; He, Yuhang

    2010-03-01

    A higher precise inner crack three-dimenssional(3D) reconstructed method of rock specimens is presented. Two inner shapes of the crack are measured with Phase Measuring Profilometry (PMP), and their edges are drawn out by height information instead of the traditional method based on gray information. Subsequently contour matching and height matching are carried out with algorithms due to maximum correlativity. The inner width and volume of the crack are educed according to the fissure of a rock specimen's outer surface, and the 3D profile of the crack is reconstructed with a high repetitive precision superior to 20μm. The proposed method is effective for evaluating the crack's width of rock specimens in the exploitation of petroleum and natural gas with a mechanical analysis method. The experiment shows its feasibility and practicability.

  16. Model for the bit-rock interaction analysis

    NASA Astrophysics Data System (ADS)

    Palmov, Vladimir A.; Vetyukov, Yury M.

    2002-02-01

    In this report a model for drilling response of the so- called drag bits (or PDC bits) is presented. Forces acting on a single cutter are supposed to be known. Discrete and continuous cutters distribution over the bit surface are considered. Both lead to similar relations between the bit kinematics characteristics and the force factors acting on it. While the bit penetration rate into the rock is small, the force and the torque are shown to depend linearly on the ratio between the bit transitional and rotational velocities (the depth of cut per revolution). Particular cases of the bit shape are compared.

  17. Modeling the Fracturing of Rock by Fluid Injection - Comparison of Numerical and Experimental Results

    NASA Astrophysics Data System (ADS)

    Heinze, Thomas; Galvan, Boris; Miller, Stephen

    2013-04-01

    Fluid-rock interactions are mechanically fundamental to many earth processes, including fault zones and hydrothermal/volcanic systems, and to future green energy solutions such as enhanced geothermal systems and carbon capture and storage (CCS). Modeling these processes is challenging because of the strong coupling between rock fracture evolution and the consequent large changes in the hydraulic properties of the system. In this talk, we present results of a numerical model that includes poro-elastic plastic rheology (with hardening, softening, and damage), and coupled to a non-linear diffusion model for fluid pressure propagation and two-phase fluid flow. Our plane strain model is based on the poro- elastic plastic behavior of porous rock and is advanced with hardening, softening and damage using the Mohr- Coulomb failure criteria. The effective stress model of Biot (1944) is used for coupling the pore pressure and the rock behavior. Frictional hardening and cohesion softening are introduced following Vermeer and de Borst (1984) with the angle of internal friction and the cohesion as functions of the principal strain rates. The scalar damage coefficient is assumed to be a linear function of the hardening parameter. Fluid injection is modeled as a two phase mixture of water and air using the Richards equation. The theoretical model is solved using finite differences on a staggered grid. The model is benchmarked with experiments on the laboratory scale in which fluid is injected from below in a critically-stressed, dry sandstone (Stanchits et al. 2011). We simulate three experiments, a) the failure a dry specimen due to biaxial compressive loading, b) the propagation a of low pressure fluid front induced from the bottom in a critically stressed specimen, and c) the failure of a critically stressed specimen due to a high pressure fluid intrusion. Comparison of model results with the fluid injection experiments shows that the model captures most of the experimental

  18. Nondestructive Methods to Characterize Rock Mechanical Properties at Low-Temperature: Applications for Asteroid Capture Technologies

    NASA Astrophysics Data System (ADS)

    Savage, Kara A.

    Recent government initiatives and commercial activities have targeted asteroids for in situ material characterization, manipulation, and possible resource extraction. Most of these activities and missions have proposed significant robotic components, given the risks and costs associated with manned missions. To successfully execute these robotic activities, detailed mechanical characteristics of the target space bodies must be known prior to contact, in order to appropriately plan and direct the autonomous robotic protocols. Unfortunately, current estimates of asteroid mechanical properties are based on limited direct information, and significant uncertainty remains specifically concerning internal structures, strengths, and elastic properties of asteroids. One proposed method to elucidate this information is through in situ, nondestructive testing of asteroid material immediately after contact, but prior to any manipulation or resource extraction activities. While numerous nondestructive rock characterization techniques have been widely deployed for terrestrial applications, these methods must be adapted to account for unique properties of asteroid material and environmental conditions of space. For example, asteroid surface temperatures may range from -100°C to 30°C due to diurnal cycling, and these low temperatures are especially noteworthy due to their deleterious influence on non-destructive testing. As a result, this thesis investigates the effect of low temperature on the mechanical characteristics and nondestructive technique responses of rock material. Initially, a novel method to produce low temperature rock samples was developed. Dry ice and methanol cooling baths of specific formulations were used to decrease rock to temperatures ranging from -60°C to 0°C. At these temperatures, shale, chalk, and limestone rock samples were exposed to several nondestructive and conventional mechanical tests, including Schmidt hammer, ultrasonic pulse velocity, point

  19. Modeling differentiation of Karaj Dam basement igneous rocks (northern Iran)

    NASA Astrophysics Data System (ADS)

    Esmaeily, D.; M-Mashhour, R.

    2009-04-01

    The Karaj Dam basement igneous body (KDB) is located in the north of city of Karaj, 30 km from city of Tehran, which lies between 35° 50' N to 36° 05' N and between 50° 50' E to 51° 15' E. It is one of the several plutonic bodies within the E-W trending Alborz zone in northern Iran. Following the late Cretaceous orogenic movements, vast volumes of dacite, andesites and basaltic lavas with tuffaceous and other clastic sediments were deposited during Eocene time, forming Karaj Formation in central Iran and Albourz. The KDB is penetrated thorough middle and upper tuff units from Karaj Formation which is underlain by late Jurassic depositions (Shemshak Formation) and overlain by the Neogene red Conglomerates in regard to stratographic consideration. It is mainly composed of a layered series dominated by gabbro, diorite and monzonite, which is a rock sequence formed upward from the lower to upper chilled margins, respectively. The chilled margins, which have gabbroic in composition, show porphyritic texture with euhedral to subhedral plagioclase (andesine & labradorite) and pyroxene (augite) megacrysts up to 5 mm long. These rocks become coarse-grained inward and transform to equigranular texture gradually.In addition, a small fine-grained doleritic stock as well as some doleritic dykes is intrusive into the pyroclastic volcanic rocks of Karaj Formation. It is possible to observe doleritic enclaves included in the KDB, indicating that the KDB are slightly younger than the dolerites. Whole rock geochemistry and mineral chemistry of the plagioclase and pyroxene in various rock samples, suggest differentiation processes. The Mg# of the pyroxene and An% of plagioclase of the contact chilled samples can be used as an indication of the original magma and plotted between the gabbro and monzonitic samples. In addition, increasing of the Mg# within the whole rock samples from the upper of contact chilled, in comparison to the lower one, demonstrates elemental differentiation

  20. Rock-Mechanics Research. A Survey of United States Research to 1965, with a Partial Survey of Canadian Universities.

    ERIC Educational Resources Information Center

    National Academy of Sciences - National Research Council, Washington, DC.

    The results of a survey, conducted by the Committee on Rock Mechanics, to determine the status of training and research in rock mechanics in presented in this publication. In 1964 and 1965 information was gathered by questionnaires sent to industries, selected federal agencies, and universities in both the United States and Canada. Results are…

  1. Modeling magmatic accumulations in the upper crust: Metamorphic implications for the country rock

    NASA Astrophysics Data System (ADS)

    Douglas, Madison M.; Geyer, Adelina; Álvarez-Valero, Antonio M.; Martí, Joan

    2016-06-01

    Field exposures of magma chambers tend to reveal contact metamorphic aureoles in the surrounding crust, which width varies from few centimeters to kilometers. The igneous accumulation not only increases the temperature around it, but also weakens its surrounding country rock beyond the brittle-ductile transition temperature. The formation of a ductile halo around the magmatic reservoir may significantly impact into the stability and growth of the magma chamber, as well as into potential dyke injections and processes of ground deformation. In this paper, we examine how a magmatic accumulation affects the country rock through the combination of petrologic and thermal perspectives. For this, we numerically modeled (i) the conductive cooling of an instantaneously emplaced magma chamber within compositionally representative pelitic and carbonate upper crusts, and (ii) the corresponding changes in the viscosity of the host rock potentially leading to ductile regimes. We consider basaltic to rhyolitic magma chambers at different depths with oblate, prolate and spherical geometries. The resulting temperature field distribution at different time steps is integrated with crustal metamorphic effects through phase diagram modeling. Our results indicate that the geometry of the magma accumulations plays a dominant role in controlling the local metamorphic and thermal effects on the country rocks. They conclude that (i) the combination of relatively simple geothermal models with petrologic datasets can generate first order predictions for the maximum metamorphic grade and geometry of magma chamber aureoles; (ii) the possible changes in the mechanical properties of the country rock are not necessarily linked to the petrological changes in contact aureoles; and (iii) the present rheologic outcomes may be used in further studies of magma chamber stability and integrity, which may favor the understanding of the melt transfer throughout the crust.

  2. Applicability of failure criteria and empirical relations of mechanical rock properties from outcrop analogue samples for wellbore stability analyses

    NASA Astrophysics Data System (ADS)

    Reyer, D.; Philipp, S. L.

    2013-12-01

    samples from larger depths shows that it is possible to apply the obtained principal stress failure criteria on clastic and volcanic rocks, but less so for carbonates. Carbonate core samples have higher strengths and develop larger angles between fault normal and main principal stress than quarry samples. This considerably reduces the residuals between quarry failure criteria and core test results. Therefore, it is advised to use failure criteria, expressed in shear and normal stresses, for prediction of core sample failure conditions. We conclude that it is possible to apply failure criteria on samples from depth if the comparability, especially textural comparability and similar porosities, of chosen outcrop analogues samples is ensured. Applicability of empirical relations of UCS with Young's modulus and tensile strength to rocks at depths is expected. Presented results may help predict mechanical properties for in situ rocks, and thus develop suitable geomechanical models for the adaptation of the drilling strategy on rock mechanical conditions. The authors appreciate the support of 'Niedersächsisches Ministerium für Wissenschaft und Kultur' and 'Baker Hughes' within the gebo research project (http: www.gebo-nds.de).

  3. New Soft Rock Pillar Strength Formula Derived Through Parametric FEA Using a Critical State Plasticity Model

    NASA Astrophysics Data System (ADS)

    Rastiello, Giuseppe; Federico, Francesco; Screpanti, Silvio

    2015-09-01

    Many abandoned room and pillar mines have been excavated not far from the surface of large areas of important European cities. In Rome, these excavations took place at shallow depths (3-15 m below the ground surface) in weak pyroclastic soft rocks. Many of these cavities have collapsed; others appear to be in a stable condition, although an appreciable percentage of their structural components (pillars, roofs, etc.) have shown increasing signs of distress from both the morphological and mechanical points of view. In this study, the stress-strain behaviour of soft rock pillars sustaining systems of cavities under vertical loads was numerically simulated, starting from the in situ initial conditions due to excavation of the cavities. The mechanical behaviour of the constituent material of the pillar was modelled according to the Modified Cam-Clay constitutive law (elasto-plastic with strain hardening). The influence of the pillar geometry (cross-section area, shape, and height) and mechanical parameters of the soft rock on the ultimate compressive strength of the pillar as a whole was parametrically investigated first. Based on the numerical results, an original relationship for pillar strength assessment was developed. Finally, the estimated pillar strengths according to the proposed formula and well-known formulations in the literature were compared.

  4. A prediction model for uniaxial compressive strength of deteriorated pyroclastic rocks due to freeze-thaw cycle

    NASA Astrophysics Data System (ADS)

    İnce, İsmail; Fener, Mustafa

    2016-08-01

    Either directly or indirectly, building stone is exposed to diverse atmospheric interactions depending on the seasonal conditions. Due to those interactions, objects of historic and cultural heritage, as well as modern buildings, partially or completely deteriorate. Among processes involved in rock deterioration, the freeze-thaw (F-T) cycle is one of the most important. Even though pyroclastic rocks have been used as building stone worldwide due to their easy workability, they are the building stone most affected by the F-T cycle. A historical region in Central Anatolia, Turkey, Cappadoia encompasses exceptional natural wonders characterized by fairy chimneys and unique historical and cultural heritage. Human-created caves, places of worship and houses have been dug into the pyroclastic rocks, which have in turn been used in architectural construction as building stone. Using 10 pyroclastic rock samples collected from Cappadocia, we determined the rock's index-mechanical properties to develop a statistical model for estimating percentage loss of uniaxial compressive strength a critical parameter of F-T cycle's important value. We used dry density (ρd), ultrasonic velocity (Vp), point load strengths (IS(50)), and slake-durability test indexes (Id4) values of unweathered rocks in our model, which is highly reliable (R2 = 0.84) for predetermination of percentage loss of uniaxial compressive strengths of pyroclastic rocks without requiring any F-T tests.

  5. A non-Linear transport model for determining shale rock characteristics

    NASA Astrophysics Data System (ADS)

    Ali, Iftikhar; Malik, Nadeem

    2016-04-01

    Unconventional hydrocarbon reservoirs consist of tight porous rocks which are characterised by nano-scale size porous networks with ultra-low permeability [1,2]. Transport of gas through them is not well understood at the present time, and realistic transport models are needed in order to determine rock properties and for estimating future gas pressure distribution in the reservoirs. Here, we consider a recently developed non-linear gas transport equation [3], ∂p-+ U ∂p- = D ∂2p-, t > 0, (1) ∂t ∂x ∂x2 complimented with suitable initial and boundary conditions, in order to determine shale rock properties such as the permeability K, the porosity φ and the tortuosity, τ. In our new model, the apparent convection velocity, U = U(p,px), and the apparent diffusivity D = D(p), are both highly non-linear functions of the pressure. The model incorporate various flow regimes (slip, surface diffusion, transition, continuum) based upon the Knudsen number Kn, and also includes Forchchiemers turbulence correction terms. In application, the model parameters and associated compressibility factors are fully pressure dependent, giving the model more realism than previous models. See [4]. Rock properties are determined by solving an inverse problem, with model parameters adjustment to minimise the error between the model simulation and available data. It is has been found that the proposed model performs better than previous models. Results and details of the model will be presented at the conference. Corresponding author: namalik@kfupm.edu.sa and nadeem_malik@cantab.net References [1] Cui, X., Bustin, A.M. and Bustin, R., "Measurements of gas permeability and diffusivity of tight reservoir rocks: different approaches and their applications", Geofluids 9, 208-223 (2009). [2] Chiba R., Fomin S., Chugunov V., Niibori Y. and Hashida T., "Numerical Simulation of Non Fickian Diffusion and Advection in a Fractured Porous Aquifer", AIP Conference Proceedings 898, 75 (2007

  6. Rock Mechanics and Enhanced Geothermal Systems: A DOE-sponsored Workshop to Explore Research Needs

    SciTech Connect

    Francois Heuze; Peter Smeallie; Derek Elsworth; Joel L. Renner

    2003-10-01

    This workshop on rock mechanics and enhanced geothermal systems (EGS) was held in Cambridge, Mass., on June 20-21 2003, before the Soil and Rock America 2003 International Conference at MIT. Its purpose was to bring together experts in the field of rock mechanics and geothermal systems to encourage innovative thinking, explore new ideas, and identify research needs in the areas of rock mechanics and rock engineering applied to enhanced geothermal systems. The agenda is shown in Appendix A. The workshop included experts in the fields of rock mechanics and engineering, geological engineering, geophysics, drilling, the geothermal energy production from industry, universities and government agencies, and laboratories. The list of participants is shown is Appendix B. The first day consisted of formal presentations. These are summarized in Chapter 1 of the report. By the end of the first day, two broad topic areas were defined: reservoir characterization and reservoir performance. Working groups were formed for each topic. They met and reported in plenary on the second day. The working group summaries are described in Chapter 2. The final session of the workshop was devoted to reaching consensus recommendations. These recommendations are given in Chapter 3. That objective was achieved. All the working group recommendations were considered and, in order to arrive at a practical research agenda usable by the workshop sponsors, workshop recommendations were reduced to a total of seven topics. These topics were divided in three priority groups, as follows. First-priority research topics (2): {sm_bullet} Define the pre-existing and time-dependent geometry and physical characteristics of the reservoir and its fracture network. That includes the identification of hydraulically controlling fractures. {sm_bullet} Characterize the physical and chemical processes affecting the reservoir geophysical parameters and influencing the transport properties of fractures. Incorporate those

  7. Wing Rock Motion and its Flow Mechanism over a Chined-Body Configuration

    NASA Astrophysics Data System (ADS)

    Wang, Yankui; Li, Qian; Shi, Wei

    2015-11-01

    Wing rock motion is one kind of uncommanded oscillation around the body axis over the most of the aircraft at enough high angle of attack and has a strong threat to the flight safety. The purpose of this paper is to investigate the wing rock motion over a typical body-wing configuration with a chined fuselage at fixed angle of attack firstly and four kinds of wing rock motion are revealed based on the flow phenomena, namely non-oscillation, lateral deflection, limit-cycle oscillation and irregular oscillation. Simultaneously, some special relationship between the wing rock motion and the flow over the chined body configuration are discussed. In addition, the evolution of wing rock motion and its corresponding flows when the model undergoes pitching up are also given out. All the experiments have been conducted in a low-speed wind tunnel at a Reynolds number of 1.87*10E5 and angle of attack from 0deg to 65deg. National Natural Science Foundation of China(11472028) and Open fund from State Key Laboratory of Aerodynamics.

  8. Mechanical Responses to Metamorphic Fluid-Rock Reactions - Natural Examples of Weakening vs. Embrittlement

    NASA Astrophysics Data System (ADS)

    Selverstone, J.

    2006-12-01

    Metamorphic reactions can influence strain accommodation mechanisms by changing grain size and by releasing, consuming, or changing the composition of an equilibrium fluid phase. Different deformation mechanisms, in turn, can affect metamorphic reaction rates and approaches to equilibrium by changing grain size, dislocation density, the arrangement of mineral grain boundaries, and local bulk composition. Our general understanding of water weakening effects in silicate minerals might lead us to predict that dehydration reactions will contribute to enhanced crystal plasticity, and that water-consuming reactions will strengthen rocks. Natural examples of interactions between fluid-rock reactions and strain accommodation in samples from the Tauern Window, eastern Alps, illustrate cases that both support and refute these predictions. (1) Finely interlayered graphitic and nongraphitic schists show different mechanisms of strain accommodation at different stages in their history. During burial, ductile strain was localized into graphitic horizons. During decompression, however, closely spaced Mode I extension cracks and carbonic fluid inclusion (FI) planes developed throughout the graphitic layers. Nongraphitic layers lack cracks, contain aqueous FIs, and maintained strain compatibility via crystal plasticity during unroofing. During decompression, reaction between graphite and aqueous pore fluid produced increasingly carbonic fluid that inhibited dislocation climb, experienced >60% volume expansion, and promoted Mode I crack formation. In these rocks, H2O- consuming reactions thus led to embrittlement at mid-crustal depths. (2) Finely banded mafic eclogites show outcrop- and microscopic scale evidence for synchronous strain accommodation via both crystal plasticity and brittle failure at 2 GPa. These rocks also record significant heterogeneities in reaction-controlled aH2O. Layers with aH2O>0.6 initially accommodated strain by pressure solution, producing complexly zoned

  9. Model Test of Anchoring Effect on Zonal Disintegration in Deep Surrounding Rock Masses

    PubMed Central

    Chen, Xu-Guang; Zhang, Qiang-Yong; Wang, Yuan; Liu, De-Jun; Zhang, Ning

    2013-01-01

    The deep rock masses show a different mechanical behavior compared with the shallow rock masses. They are classified into alternating fractured and intact zones during the excavation, which is known as zonal disintegration. Such phenomenon is a great disaster and will induce the different excavation and anchoring methodology. In this study, a 3D geomechanics model test was conducted to research the anchoring effect of zonal disintegration. The model was constructed with anchoring in a half and nonanchoring in the other half, to compare with each other. The optical extensometer and optical sensor were adopted to measure the displacement and strain changing law in the model test. The displacement laws of the deep surrounding rocks were obtained and found to be nonmonotonic versus the distance to the periphery. Zonal disintegration occurs in the area without anchoring and did not occur in the model under anchoring condition. By contrasting the phenomenon, the anchor effect of restraining zonal disintegration was revealed. And the formation condition of zonal disintegration was decided. In the procedure of tunnel excavation, the anchor strain was found to be alternation in tension and compression. It indicates that anchor will show the nonmonotonic law during suppressing the zonal disintegration. PMID:23997683

  10. Model test of anchoring effect on zonal disintegration in deep surrounding rock masses.

    PubMed

    Chen, Xu-Guang; Zhang, Qiang-Yong; Wang, Yuan; Liu, De-Jun; Zhang, Ning

    2013-01-01

    The deep rock masses show a different mechanical behavior compared with the shallow rock masses. They are classified into alternating fractured and intact zones during the excavation, which is known as zonal disintegration. Such phenomenon is a great disaster and will induce the different excavation and anchoring methodology. In this study, a 3D geomechanics model test was conducted to research the anchoring effect of zonal disintegration. The model was constructed with anchoring in a half and nonanchoring in the other half, to compare with each other. The optical extensometer and optical sensor were adopted to measure the displacement and strain changing law in the model test. The displacement laws of the deep surrounding rocks were obtained and found to be nonmonotonic versus the distance to the periphery. Zonal disintegration occurs in the area without anchoring and did not occur in the model under anchoring condition. By contrasting the phenomenon, the anchor effect of restraining zonal disintegration was revealed. And the formation condition of zonal disintegration was decided. In the procedure of tunnel excavation, the anchor strain was found to be alternation in tension and compression. It indicates that anchor will show the nonmonotonic law during suppressing the zonal disintegration. PMID:23997683

  11. A mathematical model of microbial enhanced oil recovery (MEOR) method for mixed type rock

    SciTech Connect

    Sitnikov, A.A.; Eremin, N.A.; Ibattulin, R.R.

    1994-12-31

    This paper deals with the microbial enhanced oil recovery method. It covers: (1) Mechanism of microbial influence on the reservoir was analyzed; (2) The main groups of metabolites affected by the hydrodynamic characteristics of the reservoir were determined; (3) The criterions of use of microbial influence method on the reservoir are defined. The mathematical model of microbial influence on the reservoir was made on this basis. The injection of molasse water solution with Clostridium bacterias into the mixed type of rock was used in this model. And the results of calculations were compared with experimental data.

  12. Combining water-rock interaction experiments with reaction path and reactive transport modelling to predict reservoir rock evolution in an enhanced geothermal system

    NASA Astrophysics Data System (ADS)

    Kuesters, Tim; Mueller, Thomas; Renner, Joerg

    2016-04-01

    Reliably predicting the evolution of mechanical and chemical properties of reservoir rocks is crucial for efficient exploitation of enhanced geothermal systems (EGS). For example, dissolution and precipitation of individual rock forming minerals often result in significant volume changes, affecting the hydraulic rock properties and chemical composition of fluid and solid phases. Reactive transport models are typically used to evaluate and predict the effect of the internal feedback of these processes. However, a quantitative evaluation of chemo-mechanical interaction in polycrystalline environments is elusive due to poorly constrained kinetic data of complex mineral reactions. In addition, experimentally derived reaction rates are generally faster than reaction rates determined from natural systems, likely a consequence of the experimental design: a) determining the rate of a single process only, e.g. the dissolution of a mineral, and b) using powdered sample materials and thus providing an unrealistically high reaction surface and at the same time eliminating the restrictions on element transport faced in-situ for fairly dense rocks. In reality, multiple reactions are coupled during the alteration of a polymineralic rocks in the presence of a fluid and the rate determining process of the overall reactions is often difficult to identify. We present results of bulk rock-water interaction experiments quantifying alteration reactions between pure water and a granodiorite sample. The rock sample was chosen for its homogenous texture, small and uniform grain size (˜0.5 mm in diameter), and absence of pre-existing alteration features. The primary minerals are plagioclase (plg - 58 vol.%), quartz (qtz - 21 vol.%), K-feldspar (Kfs - 17 vol.%), biotite (bio - 3 vol.%) and white mica (wm - 1 vol.%). Three sets of batch experiments were conducted at 200 ° C to evaluate the effect of reactive surface area and different fluid path ways using (I) powders of the bulk rock with

  13. Modelling long term damage in rock slopes and impact of deglaciation

    NASA Astrophysics Data System (ADS)

    Amitrano, D.; Lacroix, P.

    2014-12-01

    Long-term observations of large active rockslides show accelerating deformation over many thousands of years since the last deglaciation. The effect of deglaciation on slope stability is however poorly understood due to (1) limited long-term observations and (2) a complex interaction between glacier retreat and hydrogeological, mechanical, and morphological processes. To assess the sensitivity of rockslide dynamics to these different processes, a model of progressive damage through intact rock mass is developed in this study, based on the finite element method. This model uses time-to-failure laws based on rock laboratory creep experiments. It is able to reproduce progressive damage localization along shear bands associated with strain rate acceleration as observed during tertiary creep. The model reproduces the different phases of deformation associated with morphologies typical of large rockslides. This model is thus suitable for simulating the dynamics of large rockslides and the transition from initiation to rapid sliding. The sensitivity of rockslide kinematics and morphology to different mechanical properties is analyzed. This analysis shows that the time evolution of the rockslide can be inferred with the knowledge of only one time parameter, independent of the knowledge of the mechanical properties of the rock mass. This parameter is here chosen as the time when the summit slope displacement has reached 10 m, a parameter that can be estimated with cosmogenic dating. The model is then used to study the effects of deglaciation on the valley flank stability and the formation of large rockslides. This study shows that the deglaciation velocity can affect the morphology of the rockslide, with the shear band of the rockslide emerging at higher elevation as the velocity decreases.We also show that the response to the deglaciation can last several thousands of years after the glacier retreat.

  14. Mechanics of graben formation in crustal rocks - A finite element analysis

    NASA Technical Reports Server (NTRS)

    Melosh, H. J.; Williams, C. A., Jr.

    1989-01-01

    The mechanics of the initial stages of graben formation are examined, showing that the configuration of a graben (a pair of antithetically dipping normal faults) is the most energetically favorable fault configuration in elastic-brittle rocks subjected to pure extension. The stress field in the vicinity of a single initial normal fault is computed with a two-dimensional FEM. It is concluded that the major factor controlling graben width is the depth of the initial fault.

  15. Modelling rock-avalanche induced impact waves: Sensitivity of the model chains to model parameters

    NASA Astrophysics Data System (ADS)

    Schaub, Yvonne; Huggel, Christian

    2014-05-01

    New lakes are forming in high-mountain areas all over the world due to glacier recession. Often they will be located below steep, destabilized flanks and are therefore exposed to impacts from rock-/ice-avalanches. Several events worldwide are known, where an outburst flood has been triggered by such an impact. In regions such as in the European Alps or in the Cordillera Blanca in Peru, where valley bottoms are densely populated, these far-travelling, high-magnitude events can result in major disasters. Usually natural hazards are assessed as single hazardous processes, for the above mentioned reasons, however, development of assessment and reproduction methods of the hazardous process chain for the purpose of hazard map generation have to be brought forward. A combination of physical process models have already been suggested and illustrated by means of lake outburst in the Cordillera Blanca, Peru, where on April 11th 2010 an ice-avalanche of approx. 300'000m3 triggered an impact wave, which overtopped the 22m freeboard of the rock-dam for 5 meters and caused and outburst flood which travelled 23 km to the city of Carhuaz. We here present a study, where we assessed the sensitivity of the model chain from ice-avalanche and impact wave to single parameters considering rock-/ice-avalanche modeling by RAMMS and impact wave modeling by IBER. Assumptions on the initial rock-/ice-avalanche volume, calibration of the friction parameters in RAMMS and assumptions on erosion considered in RAMMS were parameters tested regarding their influence on overtopping parameters that are crucial for outburst flood modeling. Further the transformation of the RAMMS-output (flow height and flow velocities on the shoreline of the lake) into an inflow-hydrograph for IBER was also considered a possible source of uncertainties. Overtopping time, volume, and wave height as much as mean and maximum discharge were considered decisive parameters for the outburst flood modeling and were therewith

  16. ACOUSTICAL IMAGING AND MECHANICAL PROPERTIES OF SOFT ROCK AND MARINE SEDIMENTS

    SciTech Connect

    Thurman E. Scott, Jr., Ph.D.; Younane Abousleiman, Ph.D.; Musharraf Zaman, Ph.D., P.E.

    2001-07-01

    Mechanically weak formations, such as chalks, high porosity sandstones, and marine sediments, pose significant problems for oil and gas operators. Problems such as compaction, subsidence, and loss of permeability can affect reservoir production operations. For example, the unexpected subsidence of the Ekofisk chalk in the North Sea required over one billion dollars to re-engineer production facilities to account for losses created during that compaction (Sulak 1991). Another problem in weak formations is that of shallow water flows (SWF). Deep water drilling operations sometimes encounter cases where the marine sediments, at shallow depths just below the seafloor, begin to uncontrollably flow up and around the drill pipe. SWF problems created a loss of $150 million for the Ursa development project in the U.S. Gulf Coast SWF (Furlow 1998a,b; 1999a,b). The goal of this project is to provide a database on both the rock mechanical properties and the geophysical properties of weak rocks and sediments. These could be used by oil and gas companies to detect, evaluate, and alleviate potential production and drilling problems. The results will be useful in, for example, pre-drill detection of events such as SWF's by allowing a correlation of seismic data (such as hazard surveys) to rock mechanical properties. The data sets could also be useful for 4-D monitoring of the compaction and subsidence of an existing reservoir and imaging the zones of damage. During the second quarter of the project the research team has: (1) completed acoustic sensor construction, (2) conducted reconnaissance tests to map the deformational behaviors of the various rocks, (3) developed a sample assembly for the measurement of dynamic elastic and poroelastic parameters during triaxial testing, and (4) conducted a detailed review of the scientific literature and compiled a bibliography of that review. During the first quarter of the project the research team acquired several rock types for testing

  17. Mechanical properties of carboniferous rocks in the Upper Silesian Coal Basin under uniaxial and triaxial compression tests

    SciTech Connect

    Bukowska, M.

    2005-04-01

    Many years' studies of geological properties of rocks from the Upper Silesian Coal Basin have resulted in acquisition of a substantial data base of mechanical parameters of rocks over the total strain range. It is found that the post-peak rock properties are closely related with the peak strength and the pre-peak properties. The relationship between the uniaxial ultimate strength, elastic modulus, and drop modulus are determined.

  18. Workshop on rock mechanics issues in repository design and performance assessment

    SciTech Connect

    1996-04-01

    The Center for Nuclear Waste Regulatory Analyses organized and hosted a workshop on ``Rock Mechanics Issues in Repository Design and Performance Assessment`` on behalf its sponsor the U.S. Nuclear Regulatory Commission (NRC). This workshop was held on September 19- 20, 1994 at the Holiday Inn Crowne Plaza, Rockville, Maryland. The objectives of the workshop were to stimulate exchange of technical information among parties actively investigating rock mechanics issues relevant to the proposed high-level waste repository at Yucca Mountain and identify/confirm rock mechanics issues important to repository design and performance assessment The workshop contained three technical sessions and two panel discussions. The participants included technical and research staffs representing the NRC and the Department of Energy and their contractors, as well as researchers from the academic, commercial, and international technical communities. These proceedings include most of the technical papers presented in the technical sessions and the transcripts for the two panel discussions. Selected papers have been indexed separately for inclusion the Energy Science and Technology Database.

  19. Modeling the growth of stylolites in sedimentary rocks

    NASA Astrophysics Data System (ADS)

    Rolland, Alexandra; Toussaint, Renaud; Baud, Patrick; Schmittbuhl, Jean; Conil, Nathalie; Koehn, Daniel; Renard, FrançOis; Gratier, Jean-Pierre

    2012-06-01

    Stylolites are ubiquitous pressure solution seams found in sedimentary rocks. Their morphology is shown to follow two self-affine regimes. Analyzing the scaling properties of their height over their average direction shows that (1) at small scale, they are self-affine surfaces with a Hurst exponent around 1, and (2) at large scale, they follow another self-affine scaling with Hurst exponent around 0.5. In the present paper, we show theoretically the influence of the main principal stress and the local geometry of the stylolitic interface on the dissolution reaction rate. We compute how it is affected by the deviation between the principal stress axis and the local interface between the rock and the soft material in the stylolite. The free energy entering in the dissolution reaction kinetics is expressed from the surface energy term and via integration from the stress perturbations due to these local misalignments. The resulting model shows the interface evolution at different stress conditions. In the stylolitic case, i.e., when the main principal stress is normal to the interface, two different stabilizing terms dominate at small and large scales which are linked respectively to the surface energy and to the elastic interactions. Integrating the presence of small-scale heterogeneities related to the rock properties of the grains in the model leads to the formulation of a Langevin equation predicting the dynamic evolution of the surface. This equation leads to saturated surfaces obeying the two observed scaling laws. Analytical and numerical analysis of this surface evolution model shows that the crossover length separating both scaling regimes depends directly on the applied far-field stress magnitude. This method gives the basis for the development of a paleostress magnitude marker. We apply the computation of this marker, i.e., the morphological analysis, on a stylolite found in the Dogger limestone layer located in the neighborhood of the ANDRA Underground

  20. Computational method for thermoviscoelasticity with application to rock mechanics. [Ph. D. Thesis

    SciTech Connect

    Lee, S.C.

    1984-01-01

    Large-scale numerical computations associated with rock mechanics problems have required efficient and economical models for predicting temperature, stress, failure, and deformed structural configuration under various loading conditions. To meet this requirement, the complex dependence of the properties of geological materials on the time and temperature is modified to yield a reduced time scale as a function of time and temperature under the thermorheologically simple material (TSM) postulate. The thermorheologically linear concept is adopted in the finite element formulation by uncoupling thermal and mechanical responses. The thermal responses, based on transient heat conduction or convective-diffusion, are formulated by using the two-point recurrence scheme and the upwinding scheme, respectively. An incremental solution procedure with the implicit time stepping scheme is proposed for the solution of the thermoviscoelastic response. The proposed thermoviscoelastic solution algorithm is based on the uniaxial creep experimental data and the corresponding temperature shift functions, and is intended to minimize computational efforts by allowing large time step size with stable solutions. A thermoelastic fracture formulation is also presented by introducing the degenerate quadratic isoparametric singular element for the thermally-induced line crack problems. The stress intensity factors are computed by use of the displacement method. Efficiency of the presented formulation and solution algorithm is initially demonstrated by comparison with other available solutions for a variety of problems. Subsequent field applications are made to simulate the post-burn and post-repose phases of an underground coal conversion (UCC) experiment and in-situ nuclear waste disposal management problems. 137 references, 48 figures, 6 tables.

  1. Evolution of stress-induced borehole breakout in inherently anisotropic rock: Insights from discrete element modeling

    NASA Astrophysics Data System (ADS)

    Duan, K.; Kwok, C. Y.

    2016-04-01

    The aim of this study is to better understand the mechanisms controlling the initiation, propagation, and ultimate pattern of borehole breakouts in shale formation when drilled parallel with and perpendicular to beddings. A two-dimensional discrete element model is constructed to explicitly represent the microstructure of inherently anisotropic rocks by inserting a series of individual smooth joints into an assembly of bonded rigid discs. Both isotropic and anisotropic hollow square-shaped samples are generated to represent the wellbores drilled perpendicular to and parallel with beddings at reduced scale. The isotropic model is validated by comparing the stress distribution around borehole wall and along X axis direction with analytical solutions. Effects of different factors including the particle size distribution, borehole diameter, far-field stress anisotropy, and rock anisotropy are systematically evaluated on the stress distribution and borehole breakout propagation. Simulation results reveal that wider particle size distribution results in the local stress perturbations which cause localization of cracks. Reduction of borehole diameter significantly alters the crack failure from tensile to shear and raises the critical pressure. Rock anisotropy plays an important role on the stress state around wellbore which lead to the formation of preferred cracks under hydrostatic stress. Far-field stress anisotropy plays a dominant role in the shape of borehole breakout when drilled perpendicular to beddings while a secondary role when drilled parallel with beddings. Results from this study can provide fundamental insights on the underlying particle-scale mechanisms for previous findings in laboratory and field on borehole stability in anisotropic rock.

  2. From minerals to rocks: Toward modeling lithologies with remote sensing

    NASA Technical Reports Server (NTRS)

    Mustard, John F.; Sunshine, Jessica M.; Pieters, Carle M.; Hoppin, Andrew; Pratt, Stephan F.

    1993-01-01

    High spectral resolution imaging spectroscopy will play an important role in future planetary missions. Sophisticated approaches will be needed to unravel subtle, super-imposed spectral features typically of natural systems, and to maximize the science return of these instruments. Carefully controlled laboratory investigations using homogeneous mineral separates have demonstrated that variations due to solid solution, changes in modal abundances, and the effects of particle size are well understood from a physical basis. In many cases, these variations can be modeled quantitatively using photometric models, mixing approaches, and deconvolution procedures. However, relative to the spectra of individual mineral components, reflectance spectra of rocks and natural surfaces exhibit a reduced spectral contrast. In addition, soils or regolith, which are likely to dominate any natural planetary surface, exhibit spectral properties that have some similarities to the parent materials, but due to weathering and alteration, differences remain that cannot yet be fully recreated in the laboratory or through mixture modeling. A significant challenge is therefore to integrate modeling approaches to derive both lithologic determinations and include the effects of alteration. We are currently conducting laboratory investigations in lithologic modeling to expand upon the basic results of previous analyses with our initial goal to more closely match physical state of natural systems. The effects of alteration are to be considered separately.

  3. Model for transient creep of southeastern New Mexico rock salt

    SciTech Connect

    Herrmann, W; Wawersik, W R; Lauson, H S

    1980-11-01

    In a previous analysis, existing experimental data pertaining to creep tests on rock salt from the Salado formation of S.E. New Mexico were fitted to an exponential transient creep law. While very early time portions of creep strain histories were not fitted very well for tests at low temperatures and stresses, initial creep rates in particular generally being underestimated, the exponential creep law has the property that the transient creep strain approaches a finite limit with time, and is therefore desirable from a creep modelling point of view. In this report, an analysis of transient creep is made. It is found that exponential transient creep can be related to steady-state creep through a universal creep curve. The resultant description is convenient for creep analyses where very early time behavior is not important.

  4. COTHERM: Modelling fluid-rock interactions in Icelandic geothermal systems

    NASA Astrophysics Data System (ADS)

    Thien, Bruno; Kosakowski, Georg; Kulik, Dmitrii

    2014-05-01

    Mineralogical alteration of reservoir rocks, driven by fluid circulation in natural or enhanced geothermal systems, is likely to influence the long-term performance of geothermal power generation. A key factor is the change of porosity due to dissolution of primary minerals and precipitation of secondary phases. Porosity changes will affect fluid circulation and solute transport, which, in turn, influence mineralogical alteration. This study is part of the Sinergia COTHERM project (COmbined hydrological, geochemical and geophysical modeling of geotTHERMal systems) that is an integrative research project aimed at improving our understanding of the sub-surface processes in magmatically-driven natural geothermal systems. We model the mineralogical and porosity evolution of Icelandic geothermal systems with 1D and 2D reactive transport models. These geothermal systems are typically high enthalphy systems where a magmatic pluton is located at a few kilometers depth. The shallow plutons increase the geothermal gradient and trigger the circulation of hydrothermal waters with a steam cap forming at shallow depth. We investigate two contrasting geothermal systems: Krafla, for which the water recharge consists of meteoritic water; and Reykjanes, for which the water recharge mainly consists of seawater. The initial rock composition is a fresh basalt. We use the GEM-Selektor geochemical modeling package [1] for calculation of kinetically controlled mineral equilibria between the rock and the ingression water. We consider basalt minerals dissolution kinetics according to Palandri & Kharaka [2]. Reactive surface areas are assumed to be geometric surface areas, and are corrected using a spherical-particle surface/mass relationship. For secondary minerals, we consider the partial equilibrium assuming that the primary mineral dissolution is slow, and the secondary mineral precipitation is fast. Comparison of our modeling results with the mineralogical assemblages observed in the

  5. ACOUSTICAL IMAGING AND MECHANICAL PROPERTIES OF SOFT ROCK AND MARINE SEDIMENTS

    SciTech Connect

    Thurman E. Scott, Jr., Ph.D.; Younane Abousleiman, Ph.D.; Musharraf Zaman, Ph.D., P.E.

    2002-11-18

    During the seven quarter of the project the research team analyzed some of the acoustic velocity data and rock deformation data. The goal is to create a series of ''deformation-velocity maps'' which can outline the types of rock deformational mechanisms which can occur at high pressures and then associate those with specific compressional or shear wave velocity signatures. During this quarter, we began to analyze both the acoustical and deformational properties of the various rock types. Some of the preliminary velocity data from the Danian chalk will be presented in this report. This rock type was selected for the initial efforts as it will be used in the tomographic imaging study outlined in Task 10. This is one of the more important rock types in the study as the Danian chalk is thought to represent an excellent analog to the Ekofisk chalk that has caused so many problems in the North Sea. Some of the preliminary acoustic velocity data obtained during this phase of the project indicates that during pore collapse and compaction of this chalk, the acoustic velocities can change by as much as 200 m/s. Theoretically, this significant velocity change should be detectable during repeated successive 3-D seismic images. In addition, research continues with an analysis of the unconsolidated sand samples at high confining pressures obtained in Task 9. The analysis of the results indicate that sands with 10% volume of fines can undergo liquefaction at lower stress conditions than sand samples which do not have fines added. This liquefaction and/or sand flow is similar to ''shallow water'' flows observed during drilling in the offshore Gulf of Mexico.

  6. Modelling kinetically controlled water-rock interactions during geothermal stimulation in typical poly-mineralic reservoir rocks from the Upper Rhine Graben, Germany

    NASA Astrophysics Data System (ADS)

    Kuesters, Tim; Mueller, Thomas; Renner, Joerg

    2015-04-01

    A quantitative understanding of geochemically controlled reaction rates and their operating reaction mechanisms is crucial for the efficient exploration and exploitation of geothermal reservoirs. For example, the volume changes associated with dissolution and precipitation reactions potentially may affect the hydraulic properties of a reservoir during the production phase. The reactivity depends on a complex interaction of various parameters such as temperature, fluid flux and chemistry, mineral composition, reactive surface areas, etc. Most shallow geothermal reservoirs are constituted by highly permeable sedimentary rocks saturated by a fluid phase. The abundance of impermeable crystalline basement rocks (magmatic and metamorphic) increases with depth. Typically, hydraulic stimulation is necessary to create fluid pathways, i.e. the permeability of the rock is increased by the generation of new fractures and the reactivation of old fractures (Enhanced Geothermal Systems, EGS). Fresh, high energy surfaces are created by this treatment, constituting potential sites for intensive water-rock interactions. An increasing number of reactive transport models using equilibrium thermodynamic data shed considerable light on water-rock interactions. However, most models simplify the involved rocks to mono-mineralic phases and/or use rate data based on powder experiments with unnatural high reactive surface areas. In this study we present a new numerical model approach to quantify the geochemical evolution and its mechanical feedback during geothermal stimulation of typical poly-mineralic reservoir rocks at elevated temperatures (150-200 °C). Rock samples representative for geothermal energy producing sites (limestone, sandstone, volcanic tuff and granite) were collected at the Upper Rhine Graben (URG) in southern Germany, i.e. one of the high potential locations for geothermal energy production in Germany. Samples have been characterized petrographically with regard to phase

  7. Experimental Investigation on the Influence of High Pressure and High Temperature on the Mechanical Properties of Deep Reservoir Rocks

    NASA Astrophysics Data System (ADS)

    Zhang, Peng; Mishra, Brijes; Heasley, Keith A.

    2015-11-01

    Deep and ultra-deep resources extraction has resulted in the challenge of drilling into high-pressure, high-temperature (HPHT) environments. Drilling challenges at such extreme conditions prompted NETL to develop a specialized ultra-deep drilling simulator (UDS) for investigating drill behavior in such conditions. Using the UDS apparatus, complex laboratory tests were performed on Carthage marble (Warsaw limestone) and Crab Orchard sandstone, which represent the rocks in the basins of the Tuscaloosa trend in southern Louisiana and the Arbuckle play in Oklahoma and North Texas. Additionally, numerical models of the UDS were developed for performing parametric analyses that would be impossible with the UDS alone. Subsequently, it was found that the input properties for these two rock types at such extreme pressure and temperature conditions were unavailable. Therefore, a suite of unconfined compressive strength, indirect tensile strength, and triaxial compression tests ( σ 1 > σ 2 = σ 3) were performed on Carthage marble and Crab Orchard sandstone for investigating their behavior in HPHT environments. The HPHT experiments were performed at confining pressures ranging from atmospheric to 200 MPa, and with temperatures ranging from 25 to 180 °C. The influences of confining pressure and temperature on the mechanical properties of two rocks were investigated.

  8. The global rock art database: developing a rock art reference model for the RADB system using the CIDOC CRM and Australian heritage examples

    NASA Astrophysics Data System (ADS)

    Haubt, R. A.

    2015-08-01

    The Rock Art Database (RADB) is a virtual organisation that aims to build a global rock art community. It brings together rock art enthusiasts and professionals from around the world in one centralized location through the deployed publicly available RADB Management System. This online platform allows users to share, manage and discuss rock art information and offers a new look at rock art data through the use of new technologies in rich media formats. Full access to the growing platform is currently only available for a selected group of users but it already links over 200 rock art projects around the globe. This paper forms a part of the larger Rock Art Database (RADB) project. It discusses the design stage of the RADB System and the development of a conceptual RADB Reference Model (RARM) that is used to inform the design of the Rock Art Database Management System. It examines the success and failure of international and national systems and uses the Australian heritage sector and Australian rock art as a test model to develop a method for the RADB System design. The system aims to help improve rock art management by introducing the CIDOC CRM in conjunction with a rock art specific domain model. It seeks to improve data compatibility and data sharing to help with the integration of a variety of resources to create the global Rock Art Database Management System.

  9. Investigation of steady and fluctuating pressures associated with the transonic buffeting and wing rock of a one-seventh scale model of the F-5A aircraft

    NASA Technical Reports Server (NTRS)

    Hwang, C.; Pi, W. S.

    1978-01-01

    A wind tunnel test of a 1/7 scale F-5A model is described. The pressure, force, and dynamic response measurements during buffet and wing rock are evaluated. Effects of Mach number, angle of attack, sideslip angle, and control surface settings were investigated. The mean and fluctuating static pressure data are presented and correlated with some corresponding flight test data of a F-5A aircraft. Details of the instrumentation and the specially designed support system which allowed the model to oscillate in roll to simulate wing rock are also described. A limit cycle mechanism causing wing rock was identified from this study, and this mechanism is presented.

  10. Rock Cutting Depth Model Based on Kinetic Energy of Abrasive Waterjet

    NASA Astrophysics Data System (ADS)

    Oh, Tae-Min; Cho, Gye-Chun

    2016-03-01

    Abrasive waterjets are widely used in the fields of civil and mechanical engineering for cutting a great variety of hard materials including rocks, metals, and other materials. Cutting depth is an important index to estimate operating time and cost, but it is very difficult to predict because there are a number of influential variables (e.g., energy, geometry, material, and nozzle system parameters). In this study, the cutting depth is correlated to the maximum kinetic energy expressed in terms of energy (i.e., water pressure, water flow rate, abrasive feed rate, and traverse speed), geometry (i.e., standoff distance), material (i.e., α and β), and nozzle system parameters (i.e., nozzle size, shape, and jet diffusion level). The maximum kinetic energy cutting depth model is verified with experimental test data that are obtained using one type of hard granite specimen for various parameters. The results show a unique curve for a specific rock type in a power function between cutting depth and maximum kinetic energy. The cutting depth model developed here can be very useful for estimating the process time when cutting rock using an abrasive waterjet.