“…Especially, deformation and strength anisotropy are detected in both ultrasonic pulse velocity measurements [8,9,16], Brazilian tensile [11, 18, 19], uniaxial compression [11, 17], triaxial compression [13, 20], and triaxial creep tests [9,10]. In addition, researchers reported that the heterogeneous and anisotropic nature of shale has an impact on break down pressure, fracture initiation, and fracture containment during hydraulic fracturing processes [21][22][23]. Therefore, a reliable numerical model for the mechanical characterization of shale needs to take into account these properties.As shale at the macroscale and mesoscale, layered rocks are usually modeled by the numerical methods classified into continuum-based and discontinuum-based approaches [24,25].…”