Investigation of the spatial distribution pattern of 3D microcracks in single-cracked breakage

https://doi.org/10.1016/j.ijrmms.2022.105126Get rights and content

Abstract

The implementation of multi-crack recognition requires a theoretical basis from single-crack failure results. This study investigates the spatial distribution pattern of microcracks in single-crack-damaged coal specimens subjected to uniaxial compression. The distribution mechanism of the induced microcracks is analysed in terms of macrocracks and stress states. Gaussian and Weibull distributions are selected to explore the spatial distribution of microcracks, and quantile–quantile plots are used to study the probability density plots. The results show that the updated data better demonstrate the distribution characteristics of the ellipsoid shape visually, and the microcracks are observed to converge near the envelope. Linear distribution features are found in the scattered points of the quantile–quantile plots. The goodness-of-fit tests indicate that both Gaussian and Weibull distributions can be assigned to describe the microcrack distribution on the new axis. The probability density distribution results show that the shape is preferred to obey Gaussian distribution, whereas Weibull distribution exhibits a greater possibility of shape variation. The study results can enhance the existing knowledge on crack spatial distribution of acoustic emission positioning results associated with rock materials. In addition, the statistical results can interpret the cracks by recognising them with a specified distribution pattern.

Introduction

The spatial distribution of induced microcracks is important for revealing the ruptures formed inside rock samples subjected to loading.1, 2, 3 The cluster and coalescence phenomena of microcracks can lead to the formation of macrocracks and eventually failure planes, and thus contribute to increasing the initial connectivity of a fractured rock mass.4, 5, 6 Therefore, the pattern investigation of the spatial distribution characteristics of clustered microcracks is important for the understanding of rock failure.

Various methods have been applied to investigate crack distribution, such as computed tomography,7 infrared thermal imaging technology,8 fluorescent fluid tracking technology,9 optical analysis technology,10 digital image correlation technology11 and acoustic emission (AE).6,12, 13, 14, 15, 16 In particular, AE technology monitors the energy released by microdefect dislocations in the form of sound waves.17 The source location (i.e. the microcrack location) can be defined by the characteristics of the P(S) waves detected using different sensors.14 Owing to the convenience of 3D data acquisition and monitoring operability, AE technology is widely used in rock failure monitoring and damage area analysis.15,18

The study of microcracks in rock fracturing processes has shown that microcracks are randomly generated at the beginning of the failure process, with the cluster being initiated in the form of tensile microcracks.19 Moreover, microcracks in the process zone exhibit similar mechanical mechanisms of initiation, growth and coalescing.1,20, 21, 22 For instance, the development and coalescence of microcracks in granite rock types are related to the applied loading stress.22 As the load increases, the clustered microcracks become related to inelastic dissipation, which leads to macroscopic fracture and material splitting.23 Therefore, the distribution characteristics and clustered pattern of microcracks are worthy of investigation to link the potential mechanical behaviour with the visualised microcrack area in 2D (3D in space) in the rock failure process.

A statistical investigation of crack size indicated that microcracks occur in sets that exhibit well-defined statistical properties.2 In different rock types and movement modes, datasets of original faults and new fractures showed a power-law scaling,24 and the data of microcracks were evaluated to characterise the macrocrack properties.25 Furthermore, a close relation between the lengths of microcracks and macrocracks was observed as an exponential distribution26 and power-law distribution.27 The above studies on microcracks provide meaningful guidance for understanding macroscopic rock failure. However, the detailed shape pattern of clustered microcracks has scarcely been investigated. Therefore, it is important to study approaches that can facilitate a more reliable characterisation of these clustered microcracks.

To obtain potential shape characteristics from the clustered microcracks, the density of the spatial distribution of microcracks has been investigated under different loading conditions.16,28 A hydraulic fracturing test was conducted in a laboratory set up 400 m underground, and the microcracks were projected onto three vertical planes to indicate the fracture distribution. The result obtained using the kernel density estimator on each plane showed that the high-density area presents an approximately elliptical distribution in the centre of the cracking zone, whereas the low-density area is almost scattered in the remaining space.16 In addition, to visualise the AE event in the spatial distribution over time, a 3D kernel density estimation was adopted to identify the magnitude and spatial scale of the 3D microcrack zone. The high-AE-density areas in the 3D envelope were observed to move in the rock.29 The K-means algorithm was used to divide the 3D microcracks, and the centroid of the cluster far away from the end section was considered as the crack tip position.5

At a larger scale, such as microseismic monitoring, the spatial distributions of foreshocks are usually presented in both cross-sectional and subhorizontal views.30 The main cracking zones were determined by the density contour of the aftershocks.13 In addition, regardless of the front or side view, the high-density area of microseismic events in an iron mine was observed to have an ellipsoidal shape, which was confirmed to correspond to the location where the rock mass was seriously damaged.31 Different cracking areas and shapes of microcracks have been successfully obtained, and the results confirm the significance of using clustered microcrack distribution to interpret the rock failure results. However, the distribution pattern obtained from these results, not limited to areas distinguished by high density, needs to be quantitatively investigated in 3D space.

These microcracks are often regarded as induced by the failure of one macrocrack at different scales. In most cases, however, they are attributed to the joint action of various macrocracks.5,13,32 Thus, a detailed interpretation of the results of multiple macrocracks using the overall microcracks needs to be refined. Moreover, these microcracks need to be distinguished and interpreted using an appropriate distribution relation.2,33,34 Therefore, a statistical analysis of the microcracks originating from a single macrocrack is required to propose a spatial distribution pattern of clustered microcracks. Moreover, exploration of the single-crack distribution pattern is expected to play an important role in the recognition of multi-crack distribution.6 For instance, in fields that are sensitive to crack distribution, such as the nuclear waste shielding, the triggering of local crack propagation and fracturing need to be strictly controlled.35 However, it is preferred to have crack propagate for other applications such as geothermal mining or coalbed methane mining, where the crack penetration of fracturing is favoured.36,37

AE technology is widely used for detecting material fracture under multiaxial compression, which leads to complex macrocracks and currently focuses on tensile fracture.28 Therefore, the uniaxial compression test and AE monitoring of intact rock samples were performed simultaneously to obtain the single-cracked shear failure and data sample composed of microcrack coordinates. In addition, a statistical description and an inspection of potential distributions were carried out and compared. This will help to unify the microcrack distribution description of a single-cracked failure and support the understanding of multi-cracked failure.

Section snippets

Experimental procedure

In this study, uniaxial compression tests were performed at the State Key Laboratory for Geomechanics and Deep Underground Engineering (Beijing), China. Rock samples were instrumented to simultaneously measure the applied load, induced axial, and radial strains, as well as the resulting AE signals. A schematic of the experimental setup is shown in Fig. 1(a), and the AE sensor layout is shown in Fig. 1(b).

The tested coal samples were obtained directly from drilling at an underground working

Theoretical background for spatial characterisation of microcracks

In this study, a single crack resulting from the fracture was assumed to be induced by a set of stress components, and microcracks were induced along the orientation of the stress components inside the rock sample. Note that the orientation of the stress components may not exactly coincide with the orientation of the principal stresses and strains in the anisotropic materials. Therefore, in rock samples, the three stress components interacted with each other, and the induced microcracks

Results

The information of the microcracks recorded by the AE system during testing can be used to calculate the 3D location of the microcracks. The results are shown in Fig. 4.

As shown in Fig. 4, the coordinates of the microcracks can be obtained and displayed in the original Cartesian coordinate system, where the black dots represent microcracks and the grey translucent cylinders represent the tested specimens. Many studies have been conducted to optimise the AE positioning results and improve their

Analysis and discussion

Fracturing processes, including different materials and loading paths, have been extensively studied, where Weibull distribution has been found to be applicable to the evaluation of mechanical properties of rock materials, such as strength characteristics, elastic modulus, AE energy and damage parameters.43, 44, 45, 46, 47, 48 In addition, the phenomenon that cluster microcracks generally exhibit an ellipsoidal shape in a 2D plane has been previously supported using AE positioning.13,16,30,32,35

Conclusions

This paper proposes the spatial distribution pattern of 3D microcracks to provide evidence and supporting data for the application of crack recognition theory. The statistical distribution of the microcrack patterns was investigated using Gaussian and Weibull distributions. Based on the discussion presented in this paper, the following conclusions can be drawn.

  • 1.

    The stress-induced microcracks of a single-cracked failure were clustered inside the intact specimen and distributed along the

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgement

The authors would like to acknowledge the financial support from the National Natural Science Foundation of China, China (Grant nos. 52074294, 51574246) and Fundamental Research Funds for the Central Universities, China (Grant no. 2011QZ01).

References (49)

  • T. Wong et al.

    Microcrack statistics, Weibull distribution and micromechanical modeling of compressive failure in rock

    Mech Mater

    (2006)
  • H. Zhou et al.

    Acoustic emission based mechanical behaviors of Beishan granite under conventional triaxial compression and hydro-mechanical coupling tests

    Int J Rock Mech Min

    (2019)
  • M. Johnson

    Waveform based clustering and classification of AE transients in composite laminates using principal component analysis

    NDT E Int

    (2002)
  • R. Anay et al.

    Identification of damage mechanisms in cement paste based on acoustic emission

    Construct Build Mater

    (2018)
  • A. Behnia et al.

    Advanced damage detection technique by integration of unsupervised clustering into acoustic emission

    Eng Fract Mech

    (2019)
  • R. Vidya Sagar et al.

    Statistical analysis of acoustic emissions generated during unconfined uniaxial compression of cementitious materials

    Construct Build Mater

    (2019)
  • L. Wong et al.

    A method for multiscale interpretation of fracture processes in carrara marble specimen containing a single flaw under uniaxial compression

    J Geophys Res Solid Earth

    (2018)
  • X. Zhou et al.

    Effects of microfracture on wave propagation through rock mass

    Int J GeoMech

    (2017)
  • C. Wang et al.

    Three-dimensional crack recognition by unsupervised machine learning

    Rock Mech Rock Eng

    (2021)
  • P.M. Benson et al.

    Imaging slow failure in triaxially deformed Etna basalt using 3D acoustic-emission location and X-ray computed tomography

    Geophys Res Lett

    (2007)
  • M. Nasseri et al.

    Fracture toughness measurements and acoustic emission activity in brittle rocks

    Pure Appl Geophys

    (2006)
  • X. Zhou et al.

    Fracturing behavior study of three-flawed specimens by uniaxial compression and 3D digital image correlation: sensitivity to brittleness

    Rock Mech Rock Eng

    (2019)
  • S. Zhang et al.

    Acoustic emission associated with self-sustaining failure in low-porosity sandstone under uniaxial compression

    Rock Mech Rock Eng

    (2018)
  • Y. Yabe et al.

    Nucleation process of an M2 earthquake in a deep gold mine in South Africa inferred from on-fault foreshock activity

    J Geophys Res Solid Earth

    (2015)
  • Cited by (15)

    • Experimental investigation on the spatio-temporal-energy evolution pattern of limestone fracture using acoustic emission monitoring

      2022, Journal of Applied Geophysics
      Citation Excerpt :

      The edge sensor is placed 20 mm from each end of the specimen to reduce the effect of the end effect. Apply a proper amount of couplant between the sensor and the surface of the rock sample to ensure sufficient contact, and fix it to set the threshold of the AE monitoring system to 45 dB to reduce the influence of environmental noise (Rodríguez et al., 2016; Wang et al., 2022). The sampling frequency is 1 MHz, and the preamplifier gain is Set to 40 dB.

    View all citing articles on Scopus
    View full text