Identification of rock and fracture kinematics in high alpine rockwalls under the influence of elevation

In alpine environments, tectonic processes, past glaciation and weathering processes fracture rock and prepare or trigger rockfalls, which are important processes of rock slope evolution and natural hazards. In this study, I quantify thermally and ice-induced rock and fracture kinematics and place these in the context of their role in producing rockfall and climate change. I conducted laboratory measurements on intact rock samples and installed temperature loggers and crackmeters at four rockwalls reaching from 2585 to 2935 m in elevation in the Hungerli Valley, Swiss Alps. My laboratory data show that thermal expansion followed three phases of rock kinematics, which resulted in a hysteresis effect. In the field, control crackmeters on intact rock reflected these temperature phases, and based on thermal expansion coefficients of these observed phases, I modelled thermal stress. Model results show that thermal stress magnitudes were predominantly below rock strengths. Crackmeters across fractures revealed fracture opening during cooling and reverse closing behaviour during warming on daily timescales. Elevation-dependent snow cover controlled the number of daily temperature changes and thermal stresses affecting both intact and fractured rock, while the magnitude is controlled by topographic factors influencing insolation. On a seasonal scale, slow ice-segregation-induced fracture opening can occur within lithology-dependent temperature regimes called frost cracking windows. Shear plane dipping controlled whether fractures opened or closed irreversibly with time due to thermally induced block crawling on an annual scale. Climate change will shorten snow duration and increase temperature extremes and will, therefore, affect the number and the magnitude of thermal changes and associated stresses. Earlier snowmelt in combination with temperature increase will shift the ice-induced kinematic processes to higher elevations. In conclusion, climate change will affect and change rock and fracture kinematics and, therefore, change rockfall patterns in alpine environments. Future work should quantify rockfall patterns and link these patterns to climatic drivers.

Characterization of Hoek–Brown constant mi of quasi-isotropic intact rock using rigidity index approach

An accurate determination of Hoek–Brown constant mi is of great significance in the estimation of the failure criteria of brittle rock materials. So far, different approaches such as rigidity index method (R-index), uniaxial compressive strength-based method, and tensile strength-based method, and the combination of these two methods (combination based method) have been proposed to calculate the value of mi. This paper aims to thoroughly review the previously existing methods to calculate the value of mi and make comparison between the obtain results to propose the new material constants that provide the best fit with the experimental data. In order to fulfill this goal, a large number of data for different quasi-isotropic intact rock types from the literature were collected and statistically analyzed. Additionally, based on rock types, new material constants are introduced for igneous, sedimentary, and metamorphic rocks. The obtained results proves that for different rock groups (igneous, sedimentary, and metamorphic rocks), R-index method provides the best fit with the experimental data among the others, and it is also independent of rock type. Interestingly enough, there is significant differences in the predicted mi values using different methods, which is more probably due to the quantity and quality of data used in the statistical analysis.

ADAPTATION OF CRUSHED ROCK ANALYSIS TO INTACT ROCK ANALYSIS FOR IMPROVING WATER SATURATION ASSESSMENT AND FAST PRESSURE DECAY PERMEABILITY QUANTIFICATION

Legacy crushed rock analysis, as applied to unconventional formations, has shown great success in evaluating total porosity and water saturation over the previous three decades. The procedure of crushing rock into small particles improves the efficiency of fluid recovery and grain volume measurements in a laboratory environment. However, a caveat to crushed rock analysis is that water and volatile hydrocarbon evaporate from the rock during the preparatory crushing process, causing significant uncertainty in water saturation assessment. A modified crushed rock analysis incorporates nuclear magnetic resonance (NMR) measurements before and after the crushing process to quantify the volume of fluid loss. The advancements improve the overall total saturation quantification. However, challenges remain in the quantification of partitioned water and hydrocarbon loss currently derived from NMR spectrum along with its uncertainty. Furthermore, pressure decay permeability from crushed rock analysis has been reported to have two to three orders of magnitude difference between different labs. The calculated pressure decay permeability of the same rock could even vary several orders of magnitude difference with different crushed size, which questions the quality of the crushed pressure decay permeability. In this paper, we introduce an intact rock analysis workflow on unconventional cores for improved assessment of water saturation and enhanced quantification of fast pressure decay matrix permeability from intact rock. The workflow starts with acquisition of NMR T2 and bulk density measurements on the as-received state intact rock. Instead of crushing the rock, the intact rock is directly transferred to a retort chamber and heated to 300 °C for thermal extraction. The volumes of thermally-recovered fluids are quantified through an image-based process. The grain volume measurement and a second NMR T2 measurement are performed on post retort intact rock. The pressure decay curve during grain volume measurement is then used for calculating pressure decay matrix permeability. Total porosity is calculated using bulk volume and grain volume of the rock. Water saturation is quantified using total volume of recovered water. In addition, the twin as-received state rocks are processed through the crushed rock analysis workflow for an apple-to-apple comparison. Meanwhile, pressure decay permeability is cross-validated against the steady state permeability of the same sample. The introduced workflow has been successfully tested on different formations, including Bakken, Bone Spring, Eagle Ford, Cotton Valley, and Niobrara. The results show that total porosities calculated from intact rock analysis are consistent with total porosities from crushed rock analysis, while water saturations from the new workflow are average 8%SU (0.2–0.7%PU of bulk volume water) higher than those from the prior crushed rock workflow. The study also indicated that for some formations (e.g., Bone Spring) the fluid loss during crushing process is dominated by water, however, for some other formations (e.g., Bakken), hydrocarbon loss is significant. Pressure decay permeability quantified using intact rock analysis is also confirmed within an order of magnitude of steady state matrix permeability.

The effect of saturation on the physical and mechanical behavior of some rock samples

Different major factors control the strength of solid rocks. Moisture content is one of the most important factors, which can change the physical and mechanical behavior of intact rock as well as rock mass. Several early studies have shown that rock is weaker if tested wet rather than dry. In this paper, the density, P-wave velocity, uniaxial compressive strength, Brazilian tensile strength, and modulus of elasticity of seven different intact rock samples were measured under both dry and saturated conditions. The porosity of the samples was reported as well. Based on the obtained results, some correlations were proposed for estimating the saturated physical and mechanical properties from dry ones. The proposed correlations include different rock types and are more general than the previously reported ones. Comparing the obtained results showed that the mechanical and physical properties of weaker samples are more sensitive to the saturation process. Moreover, among the different mechanical parameters, Brazilian tensile strength exhibited more sensitivity to saturation. Comparing the results with the calculated porosities revealed that porosity is one of the key factors in the effect of saturation on physical and mechanical parameters. It seems that in the more porous rock samples, greater changes in the different measured parameters occur after saturation.

The Impact of Fracture Persistence and Intact Rock Bridge Failure on the In Situ Block Area Distribution

The in situ block size distribution is an essential characteristic of fractured rock masses and impacts the assessment of rockfall hazards and other fields of rock mechanics. The block size distribution can be estimated rather easily for fully persistent fractures, but it is a challenge to determine this parameter when non-persistent fractures in a rock mass should be considered. In many approaches, the block size distribution is estimated by assuming that the fractures are fully persistent, resulting in an underestimation of the block sizes for many fracture geometries. In addition, the block size distribution is influenced by intact rock bridge failure, especially in rock masses with non-persistent fractures, either in a short-term perspective during a slope failure event when the rock mass increasingly disintegrates or in a long-term view when the rock mass progressively weakens. The quantification of intact rock bridge failure in a rock mass is highly complex, comprising fracture coalescence and crack growth driven by time-dependent changes of the in situ stresses due to thermal, freezing-thawing, and pore water pressure fluctuations. This contribution presents stochastic analyses of the two-dimensional in situ block area distribution and the mean block area of non-persistent fracture networks. The applied 2D discrete fracture network approach takes into account the potential failure of intact rock bridges based on a pre-defined threshold length and relies on input parameters that can be easily measured in the field by classical discontinuity mapping methods (e.g., scanline mapping). In addition, on the basis of these discrete fracture network analyses, an empirical relationship was determined between (i) the mean block area for persistent fractures, (ii) the mean block area for non-persistent fractures, and (iii) the mean interconnectivity factor. The further adaptation of this 2D approach to 3D block geometries is discussed on the basis of general considerations. The calculations carried out in this contribution highlight the large impact of non-persistent fractures and intact rock bridge failure for rock mass characterization, e.g., rockfall assessment.