Numerical simulation and fractal analysis of mesoscopic scale failure in shale using digital images

• Published in: Journal of petroleum science & engineering Vol. 145; pp. 592 - 599
• Main Authors: Wu, Zhonghu, Zuo, Yujun, Wang, Shanyong, Yi, Tongsheng, Chen, Shijiang, Yu, Qinglei, Li, Wei, Sunwen, Jibin, Xu, Yunfei, Wang, Ruyue, Li, Bingfeng
• Format: Journal Article
• Language: English
• Published: Elsevier B.V 01.09.2016
• Subjects: Calcite; Compressive strength; Computer simulation; Digital imaging; Failure mode; Failure modes; Fractal analysis; Fractal dimension; Mathematical models; Meso-scale structure; Mesoscale phenomena; Numerical simulation; Shale; Meso-scale structure; Failure mode; Numerical simulation; Digital imaging; Fractal dimension; Shale
• ISSN: 0920-4105; 1873-4715
• DOI: 10.1016/j.petrol.2016.06.036

Shale rock formations often contain calcite, which is a mineral that plays an important role in the failure of shale and in fracturing the reservoir volume. This study uses digital images to characterise the mesoscopic-scale non-uniformity of shale resulting from the shape, size and distribution of calcite in the shale. The geometry of the calcite is mapped to a finite element mesh based on the meso-scale structure of the rock. The numerical failure analysis model RFPA2D-DIP, which can accurately represent the meso-scale structure of shale, is used to simulate uniaxial compression. The influence of the meso-scale structure on the compressive strength and the failure process under various loading conditions are studied. The simulation results show that the meso-scale structure has an important influence on the compressive strength and the ultimate failure mode of samples loaded in various directions α. Furthermore, the results show that the compressive strength varies significantly with the direction of the load, indicating anisotropy. Depending on the direction, the failure mode can take the form of an inverted V (α=0°, 30°, 45°), an inverted Z (α=15°), a line (α=75°), or a V (α=60°, 90°). The spatial distribution of the acoustic emissions (AE) we observed is similar to those observed in macroscopic failure modes. The fractal dimension (D) of the shale fractures under various loading conditions is calculated using the box-counting method. The fractal dimension is an effective indicator of the failure mode. The value of D corresponding to the various failure modes is as follows: for the inverted Z the maximum value is 1.688016, for the line the minimum value is 1.481904, and for both the inverted V and the V, the value is between 1.688016 and 1.481904. Therefore, the failure mode is more complex for larger values of D. •Digital images to characterise the mesoscopic-scale non-uniformity of shale.•The numerical failure analysis model RFPA2D-DIP to simulate uniaxial compression.•Failure modes of shale under different loading directions.•Fractal dimension of different failure modes of shale.