Laboratory simulations of fluid/gas induced micro-earthquakes: application to volcano seismology

Understanding different seismic signals recorded in active volcanic regions allows geoscientists to derive insight into the processes that generate them. A key type is known as Low Frequency or Long Period (LP) event, generally understood to be generated by different fluid types resonating in cracks and faults. The physical mechanisms of these signals have been linked to either resonance/turbulence within fluids, or as a result of fluids “sloshing” due to a mixture of gas and fluid being present in the system. Less well understood, however, is the effect of the fluid type (phase) on the measured signal. To explore this, we designed an experiment in which we generated a precisely controlled liquid to gas transition in a closed system by inducing rapid decompression of fluid-filled fault zones in a sample of basalt from Mt. Etna Volcano, Italy. We find that fluid phase transition is accompanied by a marked frequency shift in the accompanying microseismic dataset that can be compared to volcano seismic data. Moreover, our induced seismic activity occurs at pressure conditions equivalent to hydrostatic depths of 200–750 m. This is consistent with recently measured dominant frequencies of LP events and with numerous models

ISRM Suggested Method for Laboratory Acoustic Emission Monitoring

Acoustic emission (AE) is defined as high-frequency elastic waves emitted from defects such as small cracks (microcracks) within a material when stressed, typically in the laboratory. AE is a similar phenomenon to microseismicity (MS), as MS is induced by fracture of rock at an engineering scale (e.g., rockbursts in mines), that is, in the field. Thus, seismic monitoring can be applied to a wide variety of rock engineering problems, and AE is a powerful method to investigate processes of rock fracture by detecting microcracks prior to macroscopic failure and by tracking crack propagation. A basic approach involves using a single channel of data acquisition, such as with a digital oscilloscope, and analyzing the number and rate of AE events. Perhaps the most valuable information from AE is the source location, which requires recording the waveform at several sensors and determining arrival times at each. Thus, investing in a multichannel data acquisition system provides the means to monitor dynamics of the fracturing process. The purpose of this suggested method is to describe the experimental setup and devices used to monitor AE in laboratory testing of rock. The instrumentation includes the AE sensor, preamplifier, frequency (noise) filter, main amplifier, AE rate counter, and A/D (analog-to-digital) recorder, to provide fundamental knowledge on material and specimen behavior in laboratory experiments. When considering in situ seismic monitoring, the reader is referred to the relevant ISRM suggested method specifically addressing that topic (Xiao et al. 2016)

Estimation of Fracture Toughness of Anisotropic Rocks by Semi-Circular Bend (SCB) Tests Under Water Vapor Pressure

In order to investigate the influence of water vapor pressure in the surrounding environment on mode I fracture toughness (KIc) of rocks, semi-circular bend (SCB) tests under various water vapor pressures were conducted. Water vapor is one of the most effective agents which promote stress corrosion of rocks. The range of water vapor pressure used was 10−2 to 103 Pa, and two anisotropic rock types, African granodiorite and Korean granite, were used in this work. The measurement of elastic wave velocity and observation of thin sections of these rocks were performed to investigate the microstructures of the rocks. It was found that the distribution of inherent microcracks and grains have a preferred orientation. Two types of specimens in different orientations, namely Type-1 and Type-3, were prepared based on the anisotropy identified by the differences in the elastic wave velocity. KIc of both rock types was dependent on the water vapor pressure in the surrounding environment and decreased with increasing water vapor pressure. It was found that the degree of the dependence is influenced by the orientation and density of inherent microcracks. The experimental results also showed that KIc depended on the material anisotropy. A fracture process was discussed on the basis of the geometry of fractures within fractured specimens visualized by the X-ray computed tomography (CT) method. It was concluded that the dominant factor causing the anisotropy of KIc is the distribution of grains rather than inherent microcracks in these rocks

Evaluation of Mode I Fracture Toughness Assisted by the Numerical Determination of K-Resistance

The fracture toughness of a rock often varies depending on the specimen shape and the loading type used to measure it. To investigate the mode I fracture toughness using semi-circular bend (SCB) specimens, we experimentally studied the fracture toughness using SCB and chevron bend (CB) specimens, the latter being one of the specimens used extensively as an International Society for Rock Mechanics (ISRM) suggested method, for comparison. The mode I fracture toughness measured using SCB specimens is lower than both the level I and level II fracture toughness values measured using CB specimens. A numerical study based on discontinuum mechanics was conducted using a two-dimensional distinct element method (DEM) for evaluating crack propagation in the SCB specimen during loading. The numerical results indicate subcritical crack growth as well as sudden crack propagation when the load reaches the maximum. A K-resistance curve is drawn using the crack extension and the load at the point of evaluation. The fracture toughness evaluated by the K-resistance curve is in agreement with the level II fracture toughness measured using CB specimens. Therefore, the SCB specimen yields an improved value for fracture toughness when the increase of K-resistance with stable crack propagation is considered

A Criterion for Brittle Failure of Rocks Using the Theory of Critical Distances

This paper presents a new analytical criterion for brittle failure of rocks and heavily overconsolidated soils. Griffith’s model of a randomly oriented defect under a biaxial stress state is used to keep the criterion simple. The Griffith’s criterion is improved because the maximum tensile strength is not evaluated at the boundary of the defect but at a certain distance from the boundary, known as the critical distance. This fracture criterion is known as the Point Method, and is part of the Theory of Critical Distances, which is utilized in fracture mechanics. The proposed failure criterion has two parameters: the inherent tensile strength, ó0, and the ratio of the half-length of the initial crack/flaw to the critical distance, a/L. These parameters are difficult to measure but they may be correlated with the uniaxial compressive and tensile strengths, óc and ót. The proposed criterion is able to reproduce the common range of strength ratios for rocks and heavily overconsolidated soils (óc/ót=3-50) and the influence of several microstructural rock properties, such as texture and porosity. Good agreement with laboratory tests reported in the literature is found for tensile and low confining stresses.The work presented was initiated during a research project on “Structural integrity assessments of notch-type defects", for the Spanish Ministry of Science and Innovation (Ref.: MAT2010-15721)

Quantifying fracture geometry with X-ray tomography: Technique of Iterative Local Thresholding (TILT) for 3D image segmentation

This paper presents a new method—the Technique of Iterative Local Thresholding (TILT)—for processing 3D X-ray computed tomography (xCT) images for visualization and quantification of rock fractures. The TILT method includes the following advancements. First, custom masks are generated by a fracture-dilation procedure, which significantly amplifies the fracture signal on the intensity histogram used for local thresholding. Second, TILT is particularly well suited for fracture characterization in granular rocks because the multi-scale Hessian fracture (MHF) filter has been incorporated to distinguish fractures from pores in the rock matrix. Third, TILT wraps the thresholding and fracture isolation steps in an optimized iterative routine for binary segmentation, minimizing human intervention and enabling automated processing of large 3D datasets. As an illustrative example, we applied TILT to 3D xCT images of reacted and unreacted fractured limestone cores. Other segmentation methods were also applied to provide insights regarding variability in image processing. The results show that TILT significantly enhanced separability of grayscale intensities, outperformed the other methods in automation, and was successful in isolating fractures from the porous rock matrix. Because the other methods are more likely to misclassify fracture edges as void and/or have limited capacity in distinguishing fractures from pores, those methods estimated larger fracture volumes (up to 80 %), surface areas (up to 60 %), and roughness (up to a factor of 2). These differences in fracture geometry would lead to significant disparities in hydraulic permeability predictions, as determined by 2D flow simulations

Common Evolution of Mechanical and Transport Properties in Thermally Cracked Westerly Granite at Elevated Hydrostatic Pressure

Increasing the damage and crack porosity in crustal rocks can result in significant changes to various key physical properties, including mechanical strength, elastic and mechanical anisotropy, and the enhancement of transport properties. Using a Non-Interactive Crack Effective Medium (NIC) theory as a fundamental tool, we show that elastic wave dispersion can be inverted to evaluate crack density as a function of temperature and is compared with optically determined crack density. Further, we show how the existence of embedded microcrack fabrics in rocks also significantly influences the fracture toughness (KIC) of rocks as measured via a suite of tensile failure experiments (chevron cracked notch Brazilian disk). Finally, we include fluid flow in our analysis via the Guéguen and Dienes crack porosity-permeability model. Using the crack density and aspect ratio recovered from the elastic-wave velocity inversion, we successfully compare permeability evolution with pressure with the laboratory measurements of permeability