When investigating the invisible propagation of fractures in cylindrical rock samples under triaxial compression, acoustic emission (AE) monitoring is commonly used to investigate the brittle or ductile failure of granite (Lei et al. 2000; Chang and Lee 2004; Li and Li 2017) and other hard rocks (Alkan et al. 2007; Yang et al. 2012). By setting a reasonable threshold for voltage, the generated AE events can be recorded as the spatial coordinates of numerous AE signals (Hardy 1972; Grosse and Ohtsu 2008). The spatial distribution of location randomness is then analyzed by describing various AE parameters statistically and conducting harmonic time-series analysis (Hsu et al. 1978; Scruby 1987). However, a remaining challenge is how to establish a quantitative relationship between the triaxial compression behavior and the AE indicator, and this challenge is due in part to the lack of analysis of the spatial correlation among scattering AE signals.

Usually, the three-dimensional (3D) visualization of randomly distributed AE signals is done using only scatter plots without spatial correlation analysis. However, the role of the spatial correlation of discrete AE signals cannot be ignored when seeking to illustrate fracture nucleation and evolution. Kaiser (1950) first proposed the memory phenomenon of pre-generated AE signals, termed the Kaiser effect, after which Goodman and Blake (1966) emphasized the spatiotemporal evolution of AE-signal distribution (Stanchits et al. 2015; Meng et al. 2016; Rodríguez et al. 2016). As well as describing the spatial distribution of AE signals statistically, the enhanced improvement of accuracy in modeling spatiotemporal evolution is influenced strongly by the spatial correlation among AE signals (Hohl et al. 2018). The key to establishing an appropriate association between the triaxial behavior and the associated AE behavior depends completely on the strategy for modeling of the spatial correlation.

Note that whatever the spatial distribution or spatial correlation of the scattering AE signals, a comprehensive and accurate description of signal cloud is very complex and can be difficult to realize if its geometry and topology are extremely sophisticated. Both fractal theory and percolation theory are important for describing such spatial randomness. Many effective efforts have been made to describe the 3D scatter distribution of AE events (Manthei 2005; Guo et al. 2015; López-Comino et al. 2017), but without analyzing the stochastic interaction and nucleation of fracture network considering the spatial correlation (Lockner 1993). Generally, the fractal description indicates the self-similarity or self-affinity of the AE signal cloud, which is based mainly on the long-range correlation among scattering AE signals (Hirata et al. 1987; Kusunose et al. 1991; Pape et al. 1999; Xie et al. 2011).

However, to describe fracture nucleation, short-range correlation seems more appropriate than long-range correlation among AE signals (Vilhelm et al. 2008; Chevy et al. 2010; Girard et al. 2012), and an extremely important reason for this is that the post-generation of the fracture network is influenced strongly by the pre-existing ones. This dependence exists mainly in a regional 3D cube, not constrained in a simply 1D single-link. Consequently, the local connection among scattering AE signals exhibits strong short-range correlation. Some achievements referring to the critical nucleation of rock failure (Chelidze 1986; Alkan 2009; Xue et al. 2018) have shown the clear effectiveness of percolation models (Broadbent and Hammersley 1957; Shante and Kirkpatrick 1971; Chelidze 1982; Bebbington et al. 1990; Hunt et al. 2014; Yuan and Bowen 2018) for describing such spatial correlation. Also, percolation models are highly advantageous for describing critical behavior, namely the surge of accumulated AE events with a stress drop in brittle failure. Considering the short-range correlation, a method known as the cube cluster modeling of scattering AE signals is proposed herein to describe the spatial correlation.

2.1 Sample preparation

The granite used in the present study was excavated by drilling at a depth of 550 m in the Beishan region of Northwestern China, which is an area that has been pre-selected for the disposal of high-level nuclear waste (Zhao et al. 2013; Wang et al. 2018). It was then processed into cylindrical samples, each with a diameter of 50 mm and a height of 100 mm.

Beishan granite comprises three main minerals, namely quartz, mica, and feldspar. Rectangular images of the granite surfaces were captured precisely using a high-definition camera, and in those images the particle boundaries were classified based on color variation and calculated using the ImageJ software. The particles were equivalent to circular areas of various sizes, and statistical diameters were determined for each of the three minerals (Fig. 1). The corresponding average diameters of the quartz, mica, and feldspar particles were calculated as 2.26, 1.91, and 2.51 mm, respectively.

Normal distributions of particle diameters for the minerals quartz, mica, and feldspar as obtained using the ImageJ software

2.2 Triaxial compression test and acoustic-emission monitoring

Triaxial compression was applied using a rock mechanics test system (815 Flex Test GT; MTS, USA), and the AE signals were monitored in real time using an AE data acquisition and digital signal processing system (PCI-2; Physical Acoustics, USA) and eight AE sensors distributed around the outside of the steel chamber (Fig. 2). Under uniaxial compression, the granite strength was determined as being 161.47 MPa, and the Poisson ratio and elastic modulus were 0.15 and 61.20 GPa, respectively. Under triaxial compression, the confining pressure was set to 10 MPa.

Schematic of the layout of the eight acoustic emission (AE) sensors for cube covering

2.3 Cluster modeling of short-range correlation

By setting a threshold of 40 dB for the eight AE sensors, the AE signals generated in the cylinder were recorded effectively. Generally, AE signals come from the opening, closing, and friction of random fractures (Tang et al. 2000; Cai et al. 2007; Xue et al. 2020). Besides the spatial coordinates, the memory effect should also include some untapped information about the fracture surface and volume. Further consideration is that before the completeness of the main fracture network, local fractures are the main influences on adjacent zones. Therefore, originating from percolation models, a covering strategy involving various sized 3D cubes could be used to represent the local damaged zone, but the greatest difficulty is how to determine the optimum cube size. The cubes contain all the near-field AE signals, so the spatial correlations of the AE scattering signals could be transformed into cube-to-cube connections and then into cluster-to-cluster correlations.

If the covering cubes are too small, with a side length that is much less than the nearest distance between two adjacent points, then all the cubes are isolated, showing a completely discrete state. By contrast, if the cubes are large enough then they contain all the AE signals. A cube containing AE signals is labeled with a randomely positive number, otherwise it is labeled with zero. Naturally, a zero cube or a cluster of such cubes represents an undamaged part of the matrix. If the cubes are made continuously larger while satisfying the connection criterion, then clusters form containing different cubes, and these various clusters are distinguished and labeled by the increasing number of cube cluster (Fig. 3), which is calculated using the Hoshen–Kopelman (HK) algorithm (Hoshen et al. 1997; Al-Futaisi and Patzek 2003).

Schematic of the classification of various cube clusters

The process of generating cube clusters depends strongly on the connection type. A 3D cube has three geometrical elements, namely points, lines, and surfaces, therefore there are three corresponding connection types, namely point-to-point contact, line-to-line contact, and surface-to-surface contact. Herein, to represent short-range correlations, we consider only surface-to-surface contact (Sugiyama et al. 2010). Different clusters form depending on the size of the cubes that are used, so the crucial issue of the optimum size must be evaluated scientifically, and the cluster corresponding to that size is termed the AE percolation cluster.

3.1 Determination of cluster criticality

Figure 4 shows the accumulated AE evolution with axial strain, as well as a comparison with three other deformable quantities, namely axial strain, hoop strain, and volumetric strain. Obviously, there is a dilatancy-induced surge in the accumulated AE signals from the dilatant stress with the maximum compressive volumetric strain to the peak stress. After the peak stress, there is a rapid drop in stress after brittle failure, which is caused mainly by the full connection of fracture network. Consequently, the evolution of the fracture network (i.e., the connection process) is an important indicator of rock failure.

Complete stress–strain curves referring to the axial, hoop, and volumetric strains as well as the evolution of accumulated AE events and the corresponding cube clusters of the largest cluster (LC) and all clusters (AC)

The suggestion here is that the cube cluster covering the AE signals represents the fracture network. Therefore, the local-to-whole connection must be evaluated precisely. We mark six points on the AE curves to monitor the connection, and in particular the sixth point with the peak stress is selected to describe the transition of the local-to-whole connection, termed the critical connection. The critical connection is therefore an extremely important indicator for evaluating the appropriateness of the covering cubes.

The key to overcoming such a challenge of searching for the critical connection is evaluating the connections accurately and selecting a proper quantitative descriptor. The HK algorithm is used to predict the connections in real time. If the cubes are too large, then the fracture network is definitely in a state of complete connection; if they are too small, then there are no local connections. There is therefore a critical transition of the local-to-whole connection corresponding to the optimum cube size. This critical connection state has a unique pore fraction that is defined as the volume ratio of the AE percolation cluster to the cylinder and is taken as the quantitative descriptor.

The pore fraction corresponds uniquely to the cube size. To determine the critical connection, Fig. 5 shows the convergence of calculating the critical cube size by the bisection method to determine critical pore fraction and the abscissa indicates the logarithmic form of cube length, which is normalized and divided by the cylinder diameter. The pore fraction is commonly used to evaluate the similarity between the largest cluster (LC) and all clusters (AC). The curves for LC and AC are similar, indicating a negligible difference, which means that it is reasonable for LC to represent AC in the percolation model. Another observation is that the geometric difference between LC and AC disappears gradually as the covering cubes become larger. Finally, the critical cube length is determined as being 1.79 mm, which is between the third state with a length of 1.71 mm and the fourth state with a length of 1.95 mm. Figure 1 shows the normal distribution of the particle diameters for each of the three minerals. The two key normal parameters of the mean \(\mu\) and standard deviation \(\sigma\) are calculated respectively as 2.26 and 1.99 mm for quartz, 1.91 and 0.91 mm for mica, and 2.51 and 2.66 mm for feldspar. The critical cube length is consistent in a confidence interval of \(\left( {\mu - \sigma ,\mu + \sigma } \right)\) with a probability of 68.27%. That confidence probability shows that the particle scale could be the best means of observing the critical transition of the local-to-whole connection of the fracture network.

Determination of critical cube length by a the bisection of pore fraction with the log of cube length divided by cylindrical diameter and b the evolution of multiscale cube clusters for LC and AC

3.2 Cluster criticality regarding anisotropy

The comfortable cube cluster depends strongly on the criticality of the local-to-whole connection. This transition is also with other critical phenomena of stress drop or surge of cumulative AE signals, so evaluating the comfortability of the cube cluster modeling is to some extent the same as estimating the criticality behavior. In percolation models, there is a classical correlation to describe the critical behavior of

where P is the connectivity probability, p and p_c indicate the causes of critical transition of connectivity and are defined as the probability density of the fracture distribution or some other variable and the critical probability, respectively, and \(\beta\) is the critical exponent, which in 3D space equals 2 (Hestir and Long 1990; Aharony and Stauffer 2014; Sakhaee-Pour and Agrawal 2018).

The 3D cube cluster has an extremely complex geometry and topography, as well as complicated connections inside the cluster. Here, the connectivity probability is estimated from the percolation fraction (PF) to illustrate the increasing proportion of LC in AC on different scales. In AC, the PF \(\varGamma_{\text{p}}\) indicates the two randomly selected cubes connected in the same cluster and it is defined by (Jarvis et al. 2017)

where s_A is the AC size and i is the cube number.

Similarly in LC, the PF \(P_{\text{L}}\) is defined by

where s_L is the LC size.

Generally, the complex structure of cube cluster can be reduced to simple geometry to investigate its anisotropy. By the transformation of cube cluster into equivalent geometric volume, the isotropy is often observed by a sphere as well as the anisotropy described by a cylinder. The radius of the equivalent sphere is determined from the cube volume as

where N is the total number of covering cubes and \(\delta\) is the cube length of side.

Figure 6 shows the nonlinear growth of the sphere’s radius with increasing LC and AC PF. The nonlinearity is described properly by means of an exponential fit. The increase in the sphere’s radius indicates the enhanced distribution of isotropy due to the regional connection of cubes in the physical cylinder. Figure 6a, b show the intersection of two fitted exponential curves for LC and AC, respectively. The point at which the curves cross matches the critical state from the bisection method in Fig. 5. Such a fixed point indicates the stability of the critical cube cluster and the criticality of the transitional behavior of isotropy. The criticality could also be verified further by two lines fitted piecewise; Fig. 7a, b show that the point at which the two lines cross for LC is similar to that for AC. More importantly, the crossing point determined by the two fitted lines corresponds to the fixed point in the exponential description. Therefore, there is an obvious criticality in the isotropy transition from weak to strong. Note that the isotropy here is correlated strongly with the dilatancy, with strong isotropy indicating enhanced dilatancy.

Equivalent sphere-based isotropic description of a LC and b AC and the comparison of two fixed points determined by the intersection of two exponential curves

Determination and comparison of two critical transition points for a LC and b AC by the intersection of two lines

Considering the cylinder shape of the physical sample, its equivalent cylinder is adopted to describe the anisotropy. The two shape parameters that must be determined are the height h and the diameter r_cyl. As an important factor, the fracture propagation induces the surge of most of the AE signals. The direction along the line linking two farthest AE points is, with a very high probability, consistent with the fracture propagation. To reveal the anisotropy, the cylinder height is defined as the biggest distance from one AE point to the farthest one, which is calculated easily by comparing the distances of two randomly selected points. Correspondingly, the bottom radius is determined by

Figure 8a shows the height evolutions for LC and AC. Obviously, the two linear trends with increasing PF are very different. For AC, the height of the equivalent cylinder is nearly constant and very close to the original height of 100 mm. However, for LC the cylinder’s height grows linearly with increasing PF, approaching that for AC.

Equivalent cylinder-based anisotropic description of a LC and AC and the coincidence of b fixed point and c critical point dependent on the exponential and linear description

Similar to the isotropy analysis, the fixed point at the intersection of two exponential curves for LC and AC (Fig. 8b) also corresponds to the turning point between the two fitting lines (Fig. 8c). The critical transition is also obvious for the anisotropy analysis. With increasing PF, the critical cube cluster is very stable in describing the isotropy or anisotropy. This means that the critical cube dependent on the full connection of the fracture network provides an invariant view of the isotropy and anisotropy.

3.3 Cluster criticality regarding pore fraction

Considering the complete damage within the cube cluster, the pore fraction is defined as the volume ratio of the cluster to the physical cylinder. The pore fraction is not equivalent to the fracture porosity, but it does include the locally damaged zones around the fractures. Beyond the damaged cluster, the other part of the granite is undamaged, so the pore fraction here indicates the proportion of damaged volume including the fracture and its adjacent zones.

Figure 9a shows the clear increase in pore fraction with increasing PF for LC and AC. In particular, the increase in pore fraction changes rapidly from slow to fast after the fixed point. The curves fit well to the scatter data while the growth is slow but exhibit serious deviation when the growth is fast. Indeed, the difference in pore fraction between LC and AC is not obvious, with a maximum difference of 0.04 at the boundary corresponding to the critical transition.

Relationship between pore fraction and percolation fraction (PF) of LC and AC by a exponential description and b linear description

Using the linear description of the correlation as shown in Fig. 9b, the transitional boundary can be verified and calculated quantitatively. Correspondingly, the critical PF is determined as 0.70 for both LC and AC. The two fractions of LC and AC are calculated respectively as 0.09 and 0.13, which are less than 0.20, indicate that most of the granite is undamaged. This means that if the fracture connections are observed on the sample scale, then too many undamaged parts will be included, thereby leading to an improper observation scale. Instead, cluster modeling is essential for determining the most appropriate scale. Totally, such transition of pore fraction indicates the criticality of full connection of fracture network.

3.4 Cluster criticality regarding specific surface area

The specific surface area (SSA) is a key parameter when analyzing the complex geometry of the cube cluster. The SSA contains the outer and inner surfaces of cube cluster because the entire cluster is inside the cylinder. Figure 10 shows the SSA decreasing with increasing PF. The dramatic drop in the total number of cubes seriously reduces the inner area due to connections. As fractures propagate, they induce the continuous generation of AE signals and correspondingly, the cube cluster enlarges with more surface areas.

Relationship between specific surface area (SSA) and PF and its descriptive correlation with fracture propagation and dilatancy

Geometrically, the sphere has the smallest SSA of all shapes, and the analysis of cluster criticality regarding isotropy showed the effectiveness of the sphere-based description of dilatancy (Fig. 6). Here, the fracture propagation is mainly accompanied by high SSA. Generally, small PF indicates that the fracture propagation is preferably observed on the small scale, whereas large PF means that the fracture dilatancy is easily noticed on the large scale. This means that using small cubes to cover the AE signals is beneficial for observing the fracture propagation. In extreme conditions, such as without any fracture dilatancy, the classical model of a single-link cluster tends to be effective in fracture propagation. However, the challenge remains of how to describe dilatancy in cluster modeling, and the key to solving this problem is to use large cubes for the covering. However, using cubes that are too large will also induce the complete attention on the fracture dilatancy ignoring the propagation behavior. Therefore, for a balance between observing propagation and dilatancy, a proper cube-based scale is needed based on the SSA.

Figure 11 shows how the SSA evolves at the six points on the AE curve in Fig. 4. Taking the third point as the transition point, the SSA decreases linearly before that point and increases linearly after that point. Also, the SSA visualized by the ellipsoid is plotted along the fitting lines. The ellipsoid becoming larger corresponds to the rapid growth of the cluster volume. The initial decrease in the SSA indicates the increasing dominance of fracture dilatancy over propagation, after which the increase shows a stronger influence of the SSA on fracture propagation than the dilatancy. Taking each cluster to be a sphere, the corresponding six SSAs are also shown in Fig. 11. Compared with the original SSAs, these new ones are all obviously lower. The ratio of the cluster SSA to that of the equivalent sphere is calculated to evaluate the deviation from isotropic dilatancy, and there is a positive linear correlation between the SSA ratio and the volumetric strain. Indeed, the fracture propagation and opening contribute to the isotropic and anisotropic dilatancy, respectively. However, the fracture propagation obviously has by far the larger influence on the rapid increase of the SSA, and so the criticality regarding the SSA mainly reflects the transition from isotropic to anisotropic dilatancy.

Relationship between SSA of cube cluster and its equivalent sphere, as well as the evolution of the SSA ratio between the cube cluster and the equivalent sphere

4.1 Relationship between pore fraction and critical cube length

In the same way as described previously, another three granite samples were tested under triaxial compression with confining pressures of 5, 20, and 30 MPa, respectively. Figure 12 shows the same dilatancy-induced surge of accumulated AE events for these three samples, but with complex cluster evolution. All the cluster geometries and topographies are completely different. By the same bisection of pore fraction and cube size, the cube length of side under 5, 20, and 30 MPa is calculated as 2.00, 2.54, and 2.30 mm, respectively. Considering the critical length of 1.79 mm at the confining pressure of 10 MPa, the critical size shows no obvious dependence on the confining pressure. However, the four cube lengths of 2.00, 1.79, 2.54, and 2.30 mm show very little difference with the average diameter of three mineral particles in Fig. 1. This suggests strongly that in cube cluster modeling, the best scale for observing the fracture connections is the particle size.

Complete stress–strain curves and cube cluster evolution of LC and AC for three samples under confining pressures of a 5 MPa, b 20 MPa, and c 30 MPa

Figure 13 shows a linear relationship between the critical cube length and the LC pore fraction, which in turn shows that cube cluster modeling of an AE signal cloud could be used to predict quantitatively the fraction of a sample that is damaged. Generally, although it is extremely important to have a reliable model for evaluating damage, most attention is paid to extrinsic influences such as the maximum principal stress, the confining pressure, and even the deviatoric stress. The prerequisite for the establishment of such reliable model is that the constraint of confining pressure could inhibit the randomness of fracture propagation. The intrinsic geometric structure is mainly correlated with the spatial distribution of the mineral particles. Without considering the variation of inner structure, the external factor of confining pressure will be ineffective in describing the correlation with intrinsic geometric variables. As mentioned above, it is difficult to establish a certain relationship between pore fraction and confining pressure, and the key to solving the problem is a geometric description of the fracture network based on a proper model. Here, the cube cluster modeling of randomly distributed AE signals provides a direct description of the fracture connections, and the highly linear relationship between the pore fraction and the critical cube length shows a better sensitivity of cube cluster on damage description than the confining pressure.

Linear correlation between pore fraction and critical cube length

4.2 Relationship between confining pressure and multifractality

Fracture networks often exhibit fractal similarity, a feature that is also used to describe the complex geometry of cube clusters. However, if the geometry becomes extremely complicated, then it cannot be described properly using a single fractal dimension. For a more detailed description, different fractal subsets must be considered for a multi-scale observation of the self-similarity, thereby necessitating a multifractal description. The geometric complexity of the cube cluster depends strongly on the stochastic distribution of the AE scattering signals. By the volume covering method (Zhou et al. 2014), cubes of various scales are used to cover the randomly distributed AE signals, and so the corresponding fractal dimension d is determined by

where N₀ is the initial cube number, and \(\delta_{0}\) is the initial cube length.

Figure 14 shows plots of the logarithm of the normalized total cube number versus that of the normalized cube length for LC and AC. Considering the AC power-law description, there are four separate linear stages with different fractal dimensions in the corresponding local ranges. The corresponding LC plot deviates considerably from the AC one when the covering cube length is less than 2.6 mm, whereas they nearly coincide when the covering cube length is greater than 2.6 mm. The point at which they bifurcate corresponds closely to the critical cube length of 1.79 mm. The gradual variation of piecewise-linear coefficients presents the multifractality of AC. However, the branching point divides the corresponding LC curve into two parts: when the cube length of side is less than the critical length, there is almost no fractal similarity; at and after the branching point, there is obvious multifractality and so the cube cluster (which depends on the critical length) still requires a multifractal description.

Plots of logarithm of normalized covering-cube number versus logarithm of normalized cube length for LC and AC to reveal the multifractality

In the multifractal description, the singularity \(\alpha\) is commonly used to characterize the multifractal spectrum \(f\left( \alpha \right)\) by

where \(\bar{p}\) is the pore fraction taken as the observation of scale invariance. Each singularity \(\alpha\) matches the total cube number \(N\left( \alpha \right)\) uniquely, so the spectrum \(f\left( \alpha \right)\) is calculated as

Chhabra and Jensen (1989) proposed a method for calculating a multifractal spectrum by introducing the discrete form of the generalized fractal dimension \(D\left( q \right)\), namely

where q is the weight of multifractal spectrum and \(k = 1,2, \cdots ,n\) is the number of subsets. By means of a Legendre transformation, the descriptive system \(\alpha \sim f\left( \alpha \right)\) is equivalent to the system \(q \sim D\left( q \right)\). The singularity then becomes

and the multifractal spectrum becomes

By constructing a series of normalized parameters \(\mu_{k}\) as the measured family of

and substituting Eq. (13) into Eqs. (11) and (12), we have

and

Figure 15 shows the evolution of the LC multifractal spectrum at the six states for the four samples under confining pressures of 5, 10, 20, and 30 MPa, respectively. To investigate the multifractal spectrum, the difference descriptors \(\Delta \alpha = \alpha_{\hbox{max} } - \alpha_{\hbox{min} }\) and \(\Delta f = f\left( {\alpha_{\hbox{min} } } \right) - f\left( {\alpha_{\hbox{max} } } \right)\) are generally used. The first descriptor, \(\Delta \alpha\), indicates the openness of the spectrum curve; the larger \(\Delta \alpha\), the greater the distribution range of the pore fraction. For the other descriptor, \(\Delta f\), usually the larger it is, the greater the difference between the maximum and minimum pore fraction. Nearly all the spectra have \(\Delta f > 0\), which means that the number of peaks is always larger than that of the accumulated valleys and that is because some constant peaks continuously distribute.

Multifractal spectrum evolution of LC under a 5 MPa, b 10 MPa, c 20 MPa, and d 30 MPa

Finally, Fig. 16 shows the piecewise linear relationship between the singularity increment \(\Delta \alpha\) and the confining pressure. Clearly, \(\Delta \alpha\) decreases with increasing confining pressure, especially at low values of the latter. Decreasing \(\Delta \alpha\) means that the pore fraction fluctuates less, which is caused mainly by the confining constraint under high pressure. In other words, the layer clusters at various positions have similar cube numbers. In general, the pore fraction given by cube cluster modeling of randomly distributed AE signals is correlated strongly with the critical cube length. However, in the multifractal description, the distribution of the pore fraction is influenced mainly by the confining pressure.

Relationship between singularity increment and confining pressure

In geotechnical engineering, AE monitoring technology is highly promising for investigating the invisible failure process in rock. However, the lack of proper modeling of the spatial correlation of randomly distributed AE signals hinders its further development to establish a predictive correlation between the AE signal cloud and the triaxial compression behavior. With the aim of providing a new way to describe the spatial correlation, the cube cluster model was proposed and analyzed herein, and the main conclusions are as follows.

Considering the criticality of connections in the fracture network, the cube cluster model of the AE scattering signals was established effectively. Intrinsically in topology, the failure of rock under triaxial compression is an evolution of fracture connectivity. The completeness of the main fracture network is defined as a critical state of connection that corresponds to the increasing-to-decreasing transition of the deviatoric stress. The criticality of the cube cluster was verified also by the other three geometric descriptors of anisotropy, pore fraction, and SSA. The criticality of the topographic connection directly induces the geometric criticality of these three descriptors.

The cube cluster model was built on the short-range correlation of the AE scattering signals. Unlike the scattering description of a randomly distributed AE signal cloud, the cube cluster model focuses on spatial correlation instead of spatial distribution, and it originates from the scientific observation of dependence of the post-generation of AE signals on pre-ones. Unlike the traditional modeling in which a single link connects any two isolated points directly, the cube cluster model is based on the local-scale observation of the interaction of adjacent signals in the same critical cube. Having built the short-range correlation in a cube, the long-range connection of the fracture network could be realized by the surface-to-surface contact.

The pore fraction considering the locally damaged zones around the fractures is influenced by the critical cube length and the confining pressure. There is a positive linear correlation between the pore fraction and the critical cube length, and based on the multifractal spectrum of the pore-fraction distribution there is a negative linear correlation between the singularity increment and the confining pressure.