Mathematical Problems in Engineering
Volume 2019 (2019), Article ID 6208207, 10 pages
https://doi.org/10.1155/2019/6208207

Research Article

Research on the Development Mechanism of Postpeak Cracks in Sandstone under Different Confining Pressures

Jihua Zhang,¹ Lianguo Wang,² Yun Dong,¹ Yadong Chen,¹ Yang Jiang,¹ Huasheng Sun,¹ Jiarui Chen,¹ and Yuqing Fan³

¹Faculty of Architecture and Civil Engineering, Huaiyin Institute of Technology, Huai’an, Jiangsu 223001, China
²State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining & Technology, Xuzhou, Jiangsu 221116, China
³School of Mining Engineering, Guizhou University of Engineering Science, Bijie, Guizhou 551700, China

Correspondence should be addressed to Yun Dong; hadyun@163.com

Received 9 January 2019; Accepted 17 April 2019; Published 2 May 2019

Academic Editor: Surajit Kumar Paul

Copyright © 2019 Jihua Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Water inrush happens occasionally during deeper roadways excavation. It is mainly due to the lack of understanding in the formation and development mechanism of cracks and its spatial distribution pattern under ground stress. In view of this, this paper used different stress levels to represent the fracture state of different parts of the surrounding rock of the deep roadway; CT detection technique is used to scan the fractured sandstone in the postpeak state; and the CT images under different confining pressures are thus obtained. The geometrical parameters such as area, length, and width of the crack are used to describe the distribution patterns based on CT images processing technique and statistical principle. These patterns are then analyzed under varying postpeak stress levels and confining pressures. The result shows that, as the area, length, and width of the cracks get larger, number of cracks increases with decreasing stress level; at different stress levels, sandstone crack area probability density, crack length probability density, and crack width probability density form exceptional, linear, and Gaussian distribution, respectively. The amount of confining pressure affects the size of cracks and the extent of expansion. This means that the higher the confining pressure is, the easier the internal crack will be penetrated and expanded and the bigger the cracks are, and the number of cracks gets lesser. Such research results can be used to describe the propagation and evolution law of cracks under different stress states of postpeak rock, which also provide an important basis for further analysis of its permeability and the stability of roadway surrounding rock.

1. Introduction

Fractures in rock are major channels for groundwater seepage. After excavation of roadway, adjustment of stress balancing would cause changes of fracture field, which is the source of common disasters of tunnel exploration such as fracture water, aquifer, collapse column, and fault [1–3]. Objects of surrounding rock of roadway timbering are ruptured rock mass in the back of peak, and underground water enters into roadway through underground watercourse, causing reduction of bearing capacity of timbering, failures of anchoring, and unstable roadway. In addition, under the effects of crustal stress, surrounding rock of roadway has different permeability in different space regions, which is closely related to number of internal fractures, space distribution, and distribution characteristics. Therefore, study on distribution rules of cracks in the back of rock peak under different stress and their quantitative expression methods has important guiding significance and application value to evaluate forms and permeability of internal cracks of rock in the back of peak and reveal disaster-causing mechanism of transfusion [4–6].

For formation and expansion mechanism of internal cracks of rock, scholars had many works. At the very beginning, they studied expansion of crack by establishing theoretical models [7–9] and then introduced rock mechanics test, numerical simulation, and CT detection. Zhao Yanlin [10–12], Wong [13–15], Nazife Erarslan [16], and Yang Shengqi et al. [17–19] studied expansion rules of rock or rock precrack with different parameters (crack density, crack porosity, crack dip angle, crack length, and crack gap) through experiment and their influence on rock mechanical property. Through DDM, discrete element, RFPA, and microscopic damage numerical method, Camonesa [20], Liang Zhengzhao [21], Yang Shengqi [22], Li Lianchong [23], Wang [24], and Lu Yinlong et al. [25, 26] studied expansion rules of rock crack and fracture, as well as destruction process of rock test pieces. Moreover, Ren Jianxi [27], Li Yanchun [28], Yin Xiaotao [29, 30], Wang Chuanyang [31], and Yang Yongming et al. [32] adopted CT scanning technology to analyze production and expansion of internal crack of rock under loading process, and their influence on rock mechanical property. The research results had important significance to reveal expansion mechanism and rules of rock crack under effects of external applied load. However, due to anisotropy of rock material, distribution changes of internal crack under the effects of external applied load would be more complex, and meanwhile, there are no quantitative description methods for quantity, form and space distribution of crack [33–36], thus making people have difficulties to grasp distribution characters of cracks under effects of crustal stress and their influence factors. To study number and geometrical morphology characteristics of crack in the back of peak under different crustal stress, the essay adopted CT detection technology to scan ruptured sandstone under different crustal stress on the basis of common triaxial test to obtain CT image of crack distribution, used image processing techniques and principle of statistics to extract geometrical morphology of crack distribution in the back of sandstone, and employed geometrical parameters including width, length, and area of crack to describe crack distribution characteristics and reveal influence of confining pressure size and stress level in the back of peak on geometrical characteristics of crack in the back of peak.

2. CT Scanning Test of Crack

After roadway is excavated, due to readjustment and distribution of surrounding rock stress, stress state in the inside of surrounding rock of roadway would be different, thus causing damage degree of internal surrounding rock in different positions has big differences. Generally speaking, the closer to the superficial part of surrounding rock, the bigger the damage degree will be; thus different stress levels of rock in the back of peak can be used to stand for fracture state in different positions around surrounding rock, which is as shown in Figure 1.

Figure 1: Relationship between loose state of surrounding rock and the stress-strain curve of rock.

The material that the test used is sandstone in certain coal mine of Huaibei mining group. Test pieces were processed into cylinders with a diameter of 50mm and a height of 100mm. 60 samples were processed, and 27 of them were selected for triaxial damage test and CT scanning test under different confining pressure and stress levels in the back of peak. To avoid test errors, unevenness degree of the two planes of processed test pieces could not surpass 0.005mm, irregularity degree of each plane cannot exceed 0.02mm, and angle between the vertical line of each plane with the central line cannot surpass 0.25°. Test samples are as shown in Figure 2.

Figure 2: Rock samples after processing.

Based on common triaxial compression test, three groups of confining pressure stresses were designed: 2MPa, 2.5 MPa, and 3MPa, corresponding to peak point, intermediate point of strain softening, and residual strength point of stress level in the back of peak, respectively (A, B, and C in Figure 1). For the convenience of operation, postpeak stress levels were 90%, 65%, and 40% of peak strength, respectively, and corresponding mean intensity was shown in Table 1.

Table 1: CT scanning experiment scheme for post-peak sandstone.

Confining pressure/MPa
Post-peak stress state /MPa

Peak intensity (A)
Strain softening (B)
Residual strength (C)

2
94
68
42

2.5
96
70
43

3
98
71
44

ACTIS300-320/225 industry CT (as shown in Figure 3(a)) was adopted to scan from of ruptured space of test pieces (see in the testing program in Table 1) under different stress levels to obtain CT images of crack distribution under different stress states (as shown in Figures 3(b) and 3(c)). ACTIS300-320/225 industry CT was mainly applied in nondestructive testing and flaw detection of industrial products, and its detection principles are as shown in Figure 4. It can closely and accurately represent internal three-dimensional structure of objects and quantificationally provide internal physical and mechanical properties such as position and size, density changes, and level of defected positions, shape, and accurate size of difference type structure, as well as internal impurities and distribution of objects.

(a) Industrial CT

(b) Cross section

(c) Longitudinal section

Figure 3: ACTIS300-320/225 industrial CT and the scanned pictures of rock specimen.

Figure 4: Detection principle of industrial CT scanning system.

For rock postpeak test pieces scanned by CT, cut out one pieces from up to bottom with an interval of 1mm and cut out 100 levels in total. Since the CT images are 1024×1024 grey-scale maps, the grey scale of changes within the scope of 0~255. To extract information of rock postpeak crack, CT images should have follow-up treatment. At present, analysis of CT images of rock is mainly to observe whether it has defects like holes and crack in the inside through images so as to roughly judge their positions, but it fails to make use of micro physical information included by the images. If contrast ratio of CT images is low, boundaries between crack and hole would be not obvious at all, and together with noises, sensory information is hard to effectively reveal macro and micro information hidden in images. Most analysis of CT data described evolution rules of crack in the inside of rock through mean value CT data or variance changes and often neglected unevenness of rock. Therefore, to highlight useful image information of rock such as holes and crack and enhance contrast ratio of images, CT images should have data treatment. Adopting Matlab software which used self-programmed procedures, the essay obtained the two-dimensional diagram including crack postpeak geometrical morphology through histogram enhancement of CT image, median filter noise reduction, and OTSU threshold segmentation and eliminated remaining isolated black and white noisy points through opening and closing operation of mathematical morphology so as to obtain binary images which can reflect postpeak crack form and distribution of sandstone, which is as shown in Figure 5. At last, it had edge detection through Canny operator to obtain microstructure outline images of postpeak sandstone with white standing for boundary between holes and crack with rock particles, just as shown in Figure 6.

Figure 5: Binary images of cracks distribution of postpeak sandstone.

Figure 6: Edge detection results of cracks distribution of postpeak sandstone.

3. Characterization Methods of Rock Crack

Mechanical behaviors of rock material has nonlinearity; thus to study postpeak distributions rules of crack, it is necessary to introduce a series of parameters to describe irregular crack forms, including distance, shape factor, length, width, and area, to have quantitative description and statistics of the above processed image crack forms.

(1) Distance. To avoid errors generated by conventional measurement of images, introduced geodesic distance to realize length measurement of irregular geometrical shape, and the definition is as follows [35]:

In the formula, stands for supremum, is the geodesic distance between a and b in the image, standing for that at least there is a shortest length of connection arc between and . Differences of geodesic distance and Euclidean distance are shown in Figure 7.

Figure 7: Measurement effect of geodesic distance and Euclidean distance.

(2) Shape Factor. Rock postpeak crack includes various irregular physical characteristics. To depict the geometrical shapes reasonably, it is necessary to introduce shape factors, which can be expressed by the following formula [29]:

In the formula, is the perimeter of irregular geometrical shape of the study area; is the acreage of the study area. If the study area is round, value is 1; if shape of X significantly deviates from roundness, value will be bigger. Perimeter of the study area is obtained by the boundary extraction algorithm of the mathematical morphology; that is to say, U(A) =(A-(AΘB)), indicating results of image subtracting the part which is eaten off by B, which is as shown in Figure 8; the acreage of the study area can be obtained by making the element number in the crack area multiplying acreage of the resolution element, which is related to precision of measuring instrument.

Figure 8: Boundary extraction algorithm of binary morphology.

(3) Crack Width. When estimating crack width in the study area, crack in curve form can be divided into several rectangular L sections to study, just as shown in Figure 9. In the figure, the solid line range is the designated area when there is no crack, while the imaginary line area is the crack range and and are length and width when there is no crack, respectively, and and are length and width, respectively, after equivalent processing.

(a) No crack

(b) After crack

Figure 9: Calculation of crack width.

By use of resolution ratio of CT machine and CT numbers when scanning the rock crack area, width of crack can be calculated [36]:

In the formula, is the acreage of the study area and ; is the average value of CT number when there is no crack in the study area; is the average value of CT number after there is crack in the study area; and -1000 is CT number of air.

4. Test Results and Analysis

Adopting the above characterization methods for crack form used self-compiled program to extract and gather all crack geometrical morphology information in CT images of processed rock postpeak and analyzed influence of confining pressure size and postpeak stress level on geometrical characteristics of above crack distribution.

(1) Distribution Rules of Crack Area. Divide crack areas obtained by scanning test pieces with CT into 8 areas from small to big and calculate probability of number of crack in each area. Distribution of crack area under different stress level is as shown in Figure 10. In the figure, peak strength, strain softening, and residual strength stand for , , and in the above test, respectively.

(a) 2MPa

(b) 2.5MPa

(c) 3MPa

Figure 10: Distribution probability of cracks area under different stress level for postpeak rock.

From Figure 10, under the same confining pressure, different postpeak stress level and distribution of crack area have big differences, and distribution of postpeak sandstone crack area under different confining stress conditions is also different. In regard of the three kinds of postpeak stress levels, namely, peak strength (A), strain softening (B), and residual strength (C), when confining pressure is 2MPa crack with an area less than 3mm2 makes up 70%, 64%, and 57%; and crack with an area bigger than 5mm2 is 5%, 10%, and 16%, respectively. When confining pressure is 2.5MPa, crack with an area less than 4mm2 makes up 67%, 61%, and 56%; and crack with an area bigger than 6mm2 makes up 4%, 7%, and 9%; when confining pressure is 3MPa, crack with an area smaller than 5mm2 makes up 77%, 73%, and 53%; and crack with an area bigger than 7mm2 makes up 4%, 9%, and 13%. Thus, we can know that, with increase of confining pressure, postpeak crack area of sandstone will gradually shift to the area with big crack, while with decrease of postpeak stress level of sandstone, small crack would continuously slip, expand and run through, and gradually shift to big crack, causing crack area gradually increase. By fitting test data, it can found that distribution of crack area nearly presents the distribution of exponential decrease and can be expressed by the following formula:

In the formula, , , and are undetermined parameters, and values of various undetermined parameters under different stress levels are as shown in Table 2.

Table 2: Fitting parameters of cracks-area distribution under different stress level for post-peak sandstone.

Parameter
Confining pressure
Confining pressure
Confining pressure

2MPa
2.5MPa
3MPa

A
B
C
A
B
C
A
B
C

0.012
0.047
0.075
0.007
0.021
0.031
0.002
0.039
0.03

0.614
0.677
0.755
0.555
0.507
0.459
0.473
0.52
0.375

1.668
1.32
1.025
2.108
2.373
2.686
2.925
2.277
3.607

(2) Distribution Rules of Crack Length. Under different postpeak stress levels, distribution of crack length is as shown in Figure 11, and it is observed that, under different postpeak stress levels, with decrease of stress level, crack length will gradually increase and its distribution probability density presents distribution of linear type decline. In regard of three kinds of postpeak stress levels, namely, peak strength (A), strain softening (B), and residual strength (C), when confining pressure is 2MPa, crack with an area smaller than 30mm makes up 67%, 63%, and 41%, respectively; and crack with an area bigger than 50mm takes up 8%, 9%, and 26%, respectively. When confining pressure is 2.5MPa, crack with an area smaller than 40mm makes up 80%, 60%, and 59%, respectively, and crack with an area bigger than 60mm makes up 10%, 14%, and 25%, respectively; when confining pressure is 3MPa, crack with an area smaller than 60mm makes up 89%, 71%, and 63%, respectively, and crack with an area bigger than 80mm makes up 3%, 11%, and 15%, respectively. Thus, we can see that postpeak stress level has a great influence on extension and distribution of crack length, and size of confining pressure affects slope or decides speed and extension direction of transition from small crack to big crack.

(a) 2MPa

(b) 2.5MPa

(c) 3MPa

Figure 11: Distribution probability of cracks length under different stress level for postpeak rock.

(3) Distribution Rules of Crack Width. Under different postpeak stress levels, distribution of crack width is as shown in Figure 12. Through fitting of test data, it can be seen that distribution probability of crack width nearly presents Gaussian distribution, which can be expressed with the following formula:

(a) 2MPa

(b) 2.5MPa

(c) 3MPa

Figure 12: Distribution probability of cracks width under different stress level for postpeak rock.

In the formula, , , , and are undermined, and values of various undetermined parameters under different postpeak stress levels are as shown in Table 3.

Table 3: Fitting parameters of cracks width distribution under different stress level for post-peak sandstone.

Parameter
Confining pressure
Confining pressure
Confining pressure

2MPa
2.5MPa
3MPa

A
B
C
A
B
C
A
B
C

y₁
0.055
0.066
0.052
0.052
0.066
0.086
0.031
0.091
0.061

A₁
0.045
0.046
0.053
0.022
0.024
0.021
0.043
0.009
0.043

0.115
0.130
0.150
0.127
0.142
0.161
0.167
0.183
0.215

0.012
0.012
0.012
0.009
0.011
0.006
0.019
0.004
0.015

From Figure 12 and Table 3, we can see, under the same confining pressure, with decrease of postpeak stress level, peak point of fitting curve would gradually increase; that is to say, crack width gradually increases; for example, when confining pressure is 3MPa, value of A, B, and C is 0.167, 0.183, and 0.215, respectively. Under the same postpeak stress level, with increase of confining pressure, peak point will increase correspondingly, for example, when they are in residual strength C, under the three kinds of confining pressure conditions, of peak point is 0.15, 0.161, and 0.215, respectively, and influence of postpeak stress level on postpeak crack width is bigger than that of confining pressure. In brief, with increase of confining pressure, narrow crack would gradually reduce and transform into wide crack.

5. Conclusions

By CT scanning sandstone test pieces under different postpeak stress levels, distribution rules of sandstone postpeak crack can be obtained:

(1) Under different confining pressure, sandstone postpeak crack area probability density, crack length probability, and crack width probability density present exponential distribution, linear distribution, and Gaussian distribution.

(2) With increase of confining pressure, crack area probability density, crack length probability, and crack width probability density gradually shift into big area, but probability density value would be smaller, indicating that confining pressure affects crack size and extension direction; that is to say, the bigger the confining pressure, the easier big crack will come into being, but the number of crack would reduce.

(3) With decrease of postpeak stress level, area, length, and width of crack would continuously increase and number of crack would also increase; small crack would slip, extend, run through, and gradually transform into big crack.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Projects nos. 51708245, 51808246, 51804129, and 51404245), the Scientific and Technological Guidance Project of Jiangsu Construction System (Projects nos. 2018ZD268 and 2017ZD246), the Science and Technology Project of the Ministry of Housing and Urban-Rural Construction of China (Project no. 2017K6004), the Natural Science Foundation for Colleges and Universities in Jiangsu Province (Projects nos. 17KJB130003, 17KJA560001, and 17KJB620002), the Project with Science and Technology of Huai’an City (Projects nos. HAG201606 and HAB201708), the Natural Science Foundation of Jiangsu Province (Projects nos. BK20160426 and BK20170457), the Department of Education Foundation of Guizhou Province (Project no. QJHKYZ273), and the Cooperation Project with Science and Technology of Guizhou Province (Project no. QKHLHZ7047).

References

N. Zhang, X. L. Xu, and G. Li, “Fissure-evolving laws of surrounding rock mass of roadway and control of seepage disasters,” Chinese Journal of Rock Mechanics and Engineering, vol. 28, no. 2, pp. 330–335, 2009.

H. Wu, N. Zhang, W. Wang, Y. Zhao, and P. Cao, “Characteristics of deformation and stress distribution of small coal pillars under leading abutment pressure,” International Journal of Mining Science and Technology, vol. 25, no. 6, pp. 921–926, 2015.

G. C. Li, Z. H. Jiang, C. X. Lv, C. Huang, G. Chen, and M. Y. Li, “Instability mechanism and control technology of soft rock roadway affected by mining and high confined water,” International Journal of Mining Science and Technology, vol. 25, no. 4, pp. 573–580, 2015.

P. Cao, H. Lin, and T. Y. Liu, “Evolution procedure of multiple rock cracks under seepage pressure,” Mathematical Problems in Engineering, vol. 1, pp. 211–244, 2013.

Y. L. Zhao, P. Cao, W. H. Ma, J. Z. Tang, W. Wang, and W. J. Wang, “Coupling model of seepage-damage-fracture in fractured rock masses and its application,” Zhongnan Daxue Xuebao, vol. 48, no. 3, pp. 794–803, 2017.

S. Krisnanto, H. Rahardjo, D. G. Fredlund, and E. C. Leong, “Mapping of cracked soils and lateral water flow characteristics through a network of cracks,” Engineering Geology, vol. 172, pp. 12–25, 2014.

A. E. Green and I. N. Sneddon, “The distribution of stress in the neighbourhood of a flat elliptical crack in an elastic solid,” Mathematical Proceedings of the Cambridge Philosophical Society, vol. 46, no. 1, pp. 159–163, 1950.

R. P. Nordgren, “Propagation of a vertical hydraulic fracture,” SPE Journal, vol. 12, no. 04, pp. 306–314, 2013.

S.-C. Lin and J. F. Abel, “Variational approach for a new direct-integration form of the virtual crack extension method,” International Journal of Fracture, vol. 38, no. 3, pp. 217–235, 1988.

Y.-L. Zhao, W. Wan, W.-J. Wang, M. Wang, and Q.-Y. Peng, “Fracture experiments on ordered multi-crack body in rock-like materials under uniaxial compression and numerical simulation of wing cracks,” Chinese Journal of Geotechnical Engineering, vol. 35, no. 11, pp. 2097–2109, 2013.

Y.-L. Zhao, W. Wan, W.-J. Wang, F.-J. Zhao, and P. Cao, “Compressive-shear rheological fracture of rock-like cracks and subcritical crack propagation test and fracture mechanism,” Chinese Journal of Geotechnical Engineering, vol. 34, no. 6, pp. 1050–1059, 2012.

Y. L. Zhao, P. Cao, W. Wang, W. J. Wang, and Q. Y. Peng, “Rock cracks subcritical propagation test and compression-shear rheological fracture model,” Journal of Central South University, vol. 45, no. 1, pp. 276–286, 2014.

R. H. C. Wong, K. T. Chau, C. A. Tang, and P. Lin, “Analysis of crack coalescence in rock-like materials containing three flaws—part I: experimental approach,” International Journal of Rock Mechanics and Mining Sciences, vol. 38, no. 7, pp. 909–924, 2001.

R. H. C. Wong and S. W. Wang, “Experiment and numerical study on the effect of material property, normal stress and the position of joint on the progressive failure under direct shear,” in NARMS-TAC, Mining and Roadwaying Innovation and Opportunity, pp. 1009–1016, Toronto, 2002.

Y. S. H. Guo, R. H. C. Wong, K. T. Chau, W. S. Zhu, and S. C. Li, “Crack growth mechanisms from 3-D surface flaw with varied dipping angle under uniaxial compression,” Key Engineering Materials, vol. 353-358, no. 4, pp. 2353–2356, 2007.

N. Erarslan, “Microstructural investigation of subcritical crack propagation and Fracture Process Zone (FPZ) by the reduction of rock fracture toughness under cyclic loading,” Engineering Geology, vol. 208, pp. 181–190, 2016.

S. Q. Yang, D. S. Yang, H. W. Jing, Y. H. Li, and S. Y. Wang, “An experimental study of the fracture coalescence behaviour of brittle sandstone specimens containing three fissures,” Rock Mechanics and Rock Engineering, vol. 45, no. 4, pp. 563–582, 2012.

S.-Q. Yang, H.-W. Jing, and T. Xu, “Mechanical behavior and failure analysis of brittle sandstone specimens containing combined flaws under uniaxial compression,” Journal of Central South University, vol. 21, no. 5, pp. 2059–2073, 2014.

S.-Q. Yang, Y.-H. Huang, W.-L. Tian, and J.-B. Zhu, “An experimental investigation on strength, deformation and crack evolution behavior of sandstone containing two oval flaws under uniaxial compression,” Engineering Geology, vol. 217, pp. 35–48, 2017.

L. A. M. Camones, E. do Amaral Vargas, R. P. de Figueiredo, and R. Q. Velloso, “Application of the discrete element method for modeling of rock crack propagation and coalescence in the step-path failure mechanism,” Engineering Geology, vol. 153, pp. 80–94, 2013.

Z. Z. Liang, L. C. Li, S. B. Tang, and YB. Zhang, “3D numerical simulation of growth of surface crack of rock specimens,” Chinese Journal of Geotechnical Engineering, vol. 33, no. 10, pp. 1615–1622, 2011 (Chinese).

S.-Q. Yang, W.-L. Tian, and Y.-H. Huang, “Failure mechanical behavior of pre-holed granite specimens after elevated temperature treatment by particle flow code,” Geothermics, vol. 72, pp. 124–137, 2018.

L. C. Li, C. A. Tang, G. Li, S. Y. Wang, Z. Z. Liang, and Y. B. Zhang, “Numerical simulation of 3D hydraulic fracturing based on an improved flow-stress-damage model and a parallel FEM technique,” Rock Mechanics and Rock Engineering, vol. 45, no. 5, pp. 801–818, 2012.

S. Y. Wang, S. W. Sloan, D. C. Sheng, and C. A. Tang, “3D numerical analysis of crack propagation of heterogeneous notched rock under uniaxial tension,” Tectonophysics, vol. 677-678, pp. 45–67, 2016.

Y. Lu and L. Wang, “Numerical simulation of mining-induced fracture evolution and water flow in coal seam floor above a confined aquifer,” Computers & Geosciences, vol. 67, pp. 157–171, 2015.

Y. L. Lu, D. Elsworth, and L. G. Wang, “Microcrack-based coupled damage and flow modeling of fracturing evolution in permeable brittle rocks,” Computers & Geosciences, vol. 49, pp. 226–244, 2013.

J. X. Ren, X. G. Feng, and H. Liu, “CT testing on meso-damage propagation mechanism of cracked-sandstone sample under triaxial compression loading,” Journal of Xi'an University of Science and Technology, vol. 29, no. 3, pp. 300–304, 2009.

T. C. Li and H. B. Lv, “CT real-time scanning tests on single crack propagation under triaxial compression,” Chinese Journal of Rock Mechanics and Engineering, vol. 29, no. 2, pp. 289–296, 2010.

X.-T. Yin, F.-N. Dang, W.-H. Ding, and H.-Q. Chen, “Morphologic measurement of crack in ct images of rock and soil,” Chinese Journal of Rock Mechanics and Engineering, vol. 25, no. 3, pp. 539–544, 2006.

X. T. Yin, S. L. Wang, F. N. Dang, W. H. Ding, and H. Q. Chen, “Study on fractal characteristics of sandstone damage-fracture under ct test condition,” Chinese Journal of Rock Mechanics and Engineering, no. S1, pp. 2721–2726, 2008.

C. Y. Wang, C. H. Yang, H. Shuai, and H. J. Mao, “CT test for evolution of mudstone fractures under compressive load,” Rock and Soil Mechanics, vol. 36, no. 6, pp. 1591–1597, 2015.

Y. Yang, Y. Ju, F. Li, F. Gao, and H. Sun, “The fractal characteristics and energy mechanism of crack propagation in tight reservoir sandstone subjected to triaxial stresses,” Journal of Natural Gas Science and Engineering, vol. 32, pp. 415–422, 2016.

Y.-M. Yang, Y. Ju, and L.-T. Mao, “Growth distribution laws and characterization methods of cracks of compact sandstone subjected to triaxial stress,” Chinese Journal of Geotechnical Engineering, vol. 36, no. 5, pp. 864–872, 2014.

L. Sui, Y. Ju, Y. Yang, Y. Yang, and A. Li, “A quantification method for shale fracability based on analytic hierarchy process,” Energy, vol. 115, pp. 637–645, 2016.

E.-K. Kim, H. Choi, K. Park, and W.-H. Kang, “Deterministic and probabilistic investigation on multiple crack interactions in a semi-infinite domain,” Mathematical Problems in Engineering, vol. 2015, Article ID 305397, 9 pages, 2015.

N. R. Backeberg, F. Iacoviello, M. Rittner et al., “Quantifying the anisotropy and tortuosity of permeable pathways in clay-rich mudstones using models based on X-ray tomography,” Scientific Reports, vol. 7, no. 1, p. 14838, 2017.