You are here
Home > BLOG > Hydrulic Fracture > Three-dimensional characterization of micro-fractures in shale reservoir rocks

Three-dimensional characterization of micro-fractures in shale reservoir rocks

Fractures are crucial for unconventional hydrocarbon exploitation, but it is difficult to accurately observe the 3D spatial distribution characteristics of fractures. Microtomography (micro-CT) technology makes it possible to observe the 3D structures of fractures at micro-scale.

Fig. 6. Statistics on the angles between fracture plane and three coordinate axes (a) Sample YP10, (b) Sample YP30 (islands removed), (c) Sample YP40.

Fig. 6. Statistics on the angles between fracture plane and three coordinate axes (a) Sample YP10, (b) Sample YP30 (islands removed), (c) Sample YP40.

The histograms of the isotropy index of Samples YP10, YP30 and YP40 are shown in Fig. 7. The isotropy index of the most clusters in both Samples YP30 and YP40 is less than 0.2, characteristic of the shapes of cracks. Comparatively, Sample YP10 has more cracks which are not ideal fracture morphology, i.e., unlike Griffith cracks. Sample YP40 has the probability of 80% for the isotropy index between 0.0 and 0.15, which can be interpreted that the cracks in Sample YP40 are closer to ideal fracture morphology than in other samples.

Fig. 7. Statistics of isotropy index of core samples. (a) Sample YP10, (b) Sample YP30 (islands removed), (c) Sample YP40.

Fig. 7. Statistics of isotropy index of core samples. (a) Sample YP10, (b) Sample YP30 (islands removed), (c) Sample YP40.

For better understanding of the shape of cracks, we use Mayavi2 software package to display the orientation matrix of clusters by ellipsoidal glyphs in Sample YP10 (Fig. 8). Fig. 8(a) displays 110 largest clusters (with≥10000 voxels) and Fig. 8(b) displays 164 small clusters (with less than 1000 and more than 700 voxels). To obtain the best visualization of the shapes of the cracks, the magnitude factor of ellipsoids used is 3 for large cracks (Fig. 8(a)) and 125 for small cracks (Fig. 8(b)), respectively.

It is clear that the relatively large clusters are very thin, showing typical crack morphology, while the small clusters have blade-like, rod-like and spherical morphological shapes. The visualization well explains the results in Fig. 7(a), in which a large number of small clusters with high isotropy index in the sample. Thus it causes that a high percentage of clusters do not have typical crack morphology.

Fig. 8. Visualization of clustered equivalent ellipsoid of Sample YP10. (a) 110 largest clusters, the magnitude factor of ellipsoids is 3; (b) 164 smallest clusters, the magnitude factor of ellipsoids is 125.

Fig. 8. Visualization of clustered equivalent ellipsoid of Sample YP10. (a) 110 largest clusters, the magnitude factor of ellipsoids is 3; (b) 164 smallest clusters, the magnitude factor of ellipsoids is 125.

The Zingg Diagram (Zingg, 1935; Voigt and Twala, 2012) is a versatile shape classification scheme that plots intermediate/long axis versus short/intermediate axis to identify if the grains are spheroids, discoids, rods or blades. We use the modified Zingg Diagram to further describe the crack shapes of the shale samples. The classified characters of clusters (cracks) are shown with (τ13)1/2 as the horizontal-axis versus (τ23)1/2 as the vertical-axis. Fig. 9 shows the statistical results of Samples YP10, YP30 and YP40.

All the points are located above the diagonal because of τ123. It can be seen that the points are mainly distributed in the range adjacent to the vertical axis and particularly large cracks concentrate in the range of (τ13)1/2<0.2 or even 0.1, indicating that most of the large cracks are close to ideal fracture morphology. Those clusters with the isotropy index greater than 0.2 may be caused by the bending of the cracks. The (τ23)1/2 values of the cracks are relatively distributed over a wide range, implying that the width of cracks change greatly.

Fig. 9. Statistics of fracture morphological characteristics. (a) Sample YP10, (b) Sample YP30 (islands removed), (c) Sample YP40.

Fig. 9. Statistics of fracture morphological characteristics. (a) Sample YP10, (b) Sample YP30 (islands removed), (c) Sample YP40.

Finally, the fitting of the fractal dimension of three samples are given below. Fig. 10 shows the graph of LN(N(δ)) against LN(δ) with box length δ from 2, 4, 8, 16, 32, 64 to 128 voxels. The plots show a very good linear relationship between LN(δ) and LN(N(δ)), and the absolute values of the slope is the fractal dimension. The correlation coefficients R are all above 0.98, which confirms that the fitted results are reliable and the characteristics of cracks do show fractal character. Fig. 10(b) compares the fittings of the fractal dimension before and after removing islands, which shows that the small pores in the sample have little effect on the fractal dimension, but only causes lower correlation coefficient.

Fig. 10. Fitting results of fracture fractal dimensions. (a) Sample YP10, (b) Sample YP30 (before and after islands being removed), (c) Sample YP40.

Fig. 10. Fitting results of fracture fractal dimensions. (a) Sample YP10, (b) Sample YP30 (before and after islands being removed), (c) Sample YP40.

Table 2 lists the fractal dimensions of all 6 analyzed samples. For these 5 samples of isothermal experiments, the fractal dimensions trend to decrease with increasing fluid pressure. This can be interpreted as some of the holes and micro-fractures in shale are closed with increasing pressure. The original Sample YS (the last row of Table 2) is not subject to any temperature or pressure, and its fractal dimension is different from that the other 5 samples.

Although all these samples are taken from the same specimen, the heterogeneity of the natural rocks may cause structural differences between the 6 samples. Thus, the results are insufficient to derive the relationship between the fractal dimension of the cracks and the pressure. To obtain the reliable law of the fractal dimensions of cracks with temperature or pressure, it is required the observation of real-time scanning of the evolution of cracks in the identical sample under high-pressure and high-temperature.

Table 2. Comparison in sample fractal dimensions.


(1)   We come to the conclusion that to effectively characterize fractures in 3D space, there are three basic requirements: 1) high-resolution 3D images; 2) individual fractures can be separated if they are intersected and form a network; 3) The scheme of fracture characterization.

(2)   Micro-tomography provides high-resolution 3D images of samples at micro-scale. We demonstrate that in this study, 6 samples with natural and artificial fractures can be effectively separated manually in spite of heavy workload. Our in-house code CTSTA is proven to be a very effective tool in characterizing fractures.

(3)   With all these prerequisites met, we fulfill a detailed characterization of micro-fractures in shale reservoir rocks from the Permian Lucaogou Formation, Junggar Basin. Our work provides detailed information of fractures including porosity, specific surface area, percolation, the shape, orientation and dimensions of individual fractures.

(4)   In addition, the fractal dimensions of fractures are extracted. The fractal dimension is the most commonly used scaling parameter that describes the characteristics of cracks across scales. Our results show that fractures in the shale samples analyzed are characteristic of fractal distribution. Furthermore, five out of the six samples that experience isothermal variable-pressure action show that the fractal dimension decreases with increasing pressure.


This work was financially supported by the National Basic Research Program of China (973 Project) “Formation mechanism and enrichment law of China’s terrestrial tight oil (shale oil)” (2014CB239004). We would like to acknowledge the support of Professor Zou Caineng, the chief scientist of the project.

Annotation symbol

Symbol Three-dimensional characterization of micro-fractures in shale reservoir rocks


Arena et al., 2014   A. Arena, C.D. Delle, J. Sarout   A new computational approach to cracks quantification from 2D image analysis: application to micro-cracks description in rocks   Comput. Geosci., 66 (1) (2014), pp. 106-120

Ambrose et al., 2010   R.J. Ambrose, R.C. Hartman, M.D. Campos, I.Y. Akkutlu, C. Sondergeld   New Pore-scale Considerations for Shale Gas in Place Calculations   Society of Petroleum Engineers, Tulsa (2010), pp. 219-229

Bai et al., 2013   B. Bai, R.G. Zhu, S.T. Wu, W.J. Yang, G. Jeff, X.X. Zhang, L. Su   Multi-scale method of Nano(Micro)-CT study on microscopic pore structure of tight sandstone of Yanchang Formation, Ordos Basin   Petrol. Explor. Dev., 40 (3) (2013), pp. 329-333

Bažant, 1997  Z.P. Bažant   Scaling of quasibrittle fracture: asymptotic analysis   Int. J. Fract., 83 (83) (1997), pp. 19-40

Borodich, 1999   F.M. Borodich   Fractals and fractal scaling in fracture mechanics   Int. J. Fract., 95 (1∼4) (1999), pp. 239-259

Bour and Davy, 1997   O. Bour, P. Davy   Connectivity of random fault networks following a power law fault length distribution   Water Resour. Res., 33 (7) (1997), pp. 1567-1583

Cao et al., 2016   T.T. Cao, Z.G. Song, S.B. Wang, J. Xia, Q.G. Cao   Micro Pore characteristics and their controlling factors of Permian shale reservoir in southern Anhui   J. Earth Sci. Environ., 38 (5) (2016), pp. 668-684   (in Chinese)

Clarkson et al., 2012   C.R. Clarkson, J. Wood, S. Burgis, S. Aquino, M. Freeman   Nanopore-structure analysis and permeability predictions for a tight gas siltstone reservoir by use of low-pressure adsorption and mercury-intrusion techniques   SPE Reservoir Eval. Eng., 15 (6) (2012), pp. 648-661

Dai et al., 2016   Q.Q. Dai, Q. Suo, C. Zhang, C.J. Lu, Y.Z. Zhang, S.J. Lu, Y. Zhao   Pore structure characteristics of tight-oil sandstone reservoir based on a new parameter measured by NMR experiment: a case study of seventh Member in Yanchang Formation, Ordos Basin   Acta Pet. Sin., 37 (7) (2016), pp. 887-897   (in Chinese)

Davy et al., 2013   P. Davy, R.L. Goc, C. Darcel   A model of fracture nucleation, growth and arrest, and consequences for fracture density and scaling   J. Geophys. Res.: Solid Earth, 118 (4) (2013), pp. 1393-1407

Griffith, 1921   A.A. Griffith   The phenomena of rupture and flow in solids   Phil. Trans. Roy. Soc. Lond., 221 (1921), pp. 163-198

Guo et al., 2014   Y.T. Guo, C.H. Yang, C.G. Jia, J.B. Xu, L. Wang, D. Li   Research on hydraulic fracturing physical simulation of shale and fracture characterization methods   Chin. J. Rock Mech. Eng., 33 (01) (2014), pp. 52-59   (in Chinese)

Jia et al., 2012   C.Z. Jia, M. Zheng, Y.F. Zhang   Unconventional hydrocarbon resources in China and the prospect of exploration and development   Petrol. Explor. Dev., 39 (2) (2012), pp. 129-136   (in Chinese)

Jia et al., 2014   C.Z. Jia, M. Zheng, Y.F. Zhang   Four important theoretical issues of unconventional petroleum geology   Acta Pet. Sin., 35 (1) (2014), pp. 1-10   (in Chinese)

Law, 2002   B.E. Law   Basin-centered gas systems   AAPG (Am. Assoc. Pet. Geol.) Bull., 86 (11) (2002), pp. 1891-1919

Lee et al., 2013   B.Y. Lee, Y.Y. Kim, S. Yi, J. Kim   Automated image processing technique for detecting and analysing concrete surface cracks   Struct. Infrastruct. Eng., 9 (6) (2013), pp. 1-11

Li et al., 2016   W.G. Li, B. Zhong, H.Z. Yang, X.F. Yang, Z.M. Hu, M. Chen   A new method for gas diffusivity evaluation in matrix rocks of shale reservoir   Acta Pet. Sin., 37 (1) (2016), pp. 88-96   (in Chinese)

Liu et al., 2016   J. Liu, G.G. Pereira, Q.B. Liu, K. Regenauer-Lieb   Computational challenges in the analyses of petrophysics using microtomography and upscaling: a review    Comput. Geosci., 89 (2016), pp. 107-117

Liu et al., 2009   J. Liu, K. Regenauer-Lieb, C. Hines, K. Liu, O. Gaede, A. Squelch   Improved estimates of percolation and anisotropic permeability from 3-D X-ray microtomography using stochastic analyses and visualization   G-cubed, 10 (5) (2009), pp. 1-13

Liu et al., 2013   J. Liu, K. Regenauer-Lieb, C. Hines, S. Zhang, P. Bourke, F. Fusseis, D.A. Yuen   Applications of Microtomography to Multiscale System Dynamics: Visualisation, Characterisation and High-performance Computation   Springer, Berlin (2013), pp. 653-674

Ma et al., 2014   Y. Ma, N.N. Zhong, X.Y. Huang, Z.P. Guo, L.P. Yao  The application of focused ion beam scanning electron microscope (FIB-SEM) to the nanometer-sized pores in shales   J. Chin. Electron Microsc. Soc., 33 (3) (2014), pp. 251-256   (in Chinese)

Qiu et al., 1999   Z.S. Qiu, R. Ding, L.X. Yu   The determing method researching of specific surface area for mud shale   Drill. Fluid Complet. Fluid, 16 (1) (1999), pp. 9-11   (in Chinese)

Russell et al., 1980   D.A. Russell, J.D. Hanson, E. Ott   Dimension of strange attractors   Phys. Rev. Lett., 45 (14) (1980), pp. 1175-1178

Schmoker, 2002   J.W. Schmoker   Resource-assessment perspectives for unconventional gas systems   AAPG (Am. Assoc. Pet. Geol.) Bull., 86 (11) (2002), pp. 1993-1999

Song, 2011   J. Song   The Research on Micro-fracture Characteristics Based on Digital Analysis Method   Chengdu University of Technology, Chengdu (2011)

Voigt and Twala, 2012   A.E. Voigt, C. Twala   Novel size and shape measurements applied to jig plant performance analysis   J. S. Afr. Inst. Min. Metall, 112 (3) (2012), pp. 171-177

Wang et al., 2012   P.R. Wang, H.W. Huang, Y.D. Xue   Automatic recognition of cracks in tunnel lining based on characteristics of local grids in images   Chin. J. Rock Mech. Eng., 31 (5) (2012), pp. 991-999   (in Chinese)

Wang, 2016   X.Z. Wang   Advances in unconventional gas exploration and development of Yanchang petroleum group   Acta Pet. Sin., 37 (1) (2016), pp. 137-144   (in Chinese)

Wilcox, 1990   R.D. Wilcox   Surface area approach key to borehole stability   Oil Gas J., 88 (9) (1990), p. 9

Yang et al., 2013   F. Yang, Z.F. Ning, C.P. Hu, B. Wang, K. Peng, H.Q. Liu   Characterization microscopic pore structures in shale reservoirs   Acta Pet. Sin., 34 (2) (2013), pp. 301-311   (in Chinese)

Yuan et al., 2016   D.D. Yuan, S.F. Lu, W.F. Chen, H. Xiao, Y.H. Wu   Shale microscopic pore structure characterization in well Pengye1 of southeast Chongqing   Special Oil Gas Reservoirs, 23 (1) (2016), pp. 49-53   (in Chinese)

Zhang et al., 2017   Q. Zhang, Z.L. Pang, R.H. Liu, et al.   Pore structure characteristics of Taiyuan formation (C3t) in Eastern uplift of Liaohe depression SW Pet. Univ. (Sci. Technol. Ed.), 39 (2) (2017), pp. 35-42   (in Chinese)

Zhao et al., 2016   J.Z. Zhao, Q. Cao, Y.B. Bai, C. Er, J. Li, W.T. Wu   Petroleum accumulation from continuous to discontinuous: concept, classification and distribution   Acta Pet. Sin., 37 (2) (2016), pp. 145-159   (in Chinese)

Zhou et al., 2004   Y.J. Zhou, J.X. Zhang, Y.S. Du, Q. Wu   Fractal dimension of core fractures in Chengbei 30 buried hill reservoir and its relation to fracture density   Oil Gas Recovery Technol., 11 (3) (2004), pp. 19-21   (in Chinese)

Zingg, 1935   T. Zingg   Beitrag zur schotteranalyse   Schweiz. Mineral. Petrogr. Mittl., 15 (1935), pp. 39-140Zou et al., 2015   C.N. Zou, G.M. Zhai, G.Y. Zhang   Formation, distribution, potential and prediction of global conventional and unconventional hydrocarbon resources   Petrol. Explor. Dev., 42 (1) (2015), pp. 14-28   (in Chinese)


* Corresponding author. School of Earth Sciences and Engineering, Sun Yat-sen

University, Guangdong, 510275, China.

E-mail addresses: [email protected] (C. Qi), [email protected] (J. Liu).

2096-2495/© 2018 Chinese Petroleum Society. This is an open access article under the CC BY-NC-ND license (

Emanuel Martin
Emanuel Martin is a Petroleum Engineer graduate from the Faculty of Engineering and a musician educate in the Arts Faculty at National University of Cuyo. In an independent way he’s researching about shale gas & tight oil and building this website to spread the scientist knowledge of the shale industry.

Leave a Reply

16 + 19 =