Anisotropic Mechanical Behaviors of Shale Rock and Their Relation to Hydraulic Fracturing in a Shale Reservoir: A Review

: Shale gas is an important supplement to the supply of natural gas resources and plays an important role on the world’s energy stage. The e ﬃ cient implementation of hydraulic fracturing is the key issue in the exploration and exploitation of shale gas. The existence of bedding structure results in a distinct anisotropy of shale rock formation. The anisotropic behaviors of shale rock have important impacts on wellbore stability, hydraulic fracture propagation, and the formation of complex fracture networks. This paper brie ﬂ y reviews previous work on the anisotropic mechanical properties of shale rock and their relation to hydraulic fracturing in shale reservoirs. In this paper, the research status of work addressing the lithological characteristics of shale rock is summarized ﬁ rst, particularly work considering the mineral constituent, which determines its physical and mechanical behavior in essence. Then the anisotropic physical and mechanical properties of shale specimens, including ultrasonic anisotropy, mechanical behavior under uniaxial and triaxial compression tests, and tensile property under the Brazilian test, are summarized, and the state of the literature on fracture toughness anisotropy is discussed. The concerns of anisotropic mechanical behavior under laboratory tests are emphasized in this paper, particularly the evaluation of shale bri tt leness based on mechanical characteristics, which is discussed in detail. Finally, further concerns such as the e ﬀ ects of bedding plane on hydraulic fracturing failure strength, crack propagation, and failure pa tt ern are also drawn out. This review study will provide a be tt er understanding of current research ﬁ ndings on the anisotropic mechanical properties of shale rock, which can provide insight into the shale anisotropy related to the fracture propagation of hydraulic fracturing in shale reservoirs.


Introduction
As a typical sedimentary rock, shale is widely distributed in the earth's crust.Shale gas has now become a strategically and globally significant unconventional resource [1,2].Hydraulic fracturing is the core technology of shale gas development and is generally employed in the exploitation of shale gas [3][4][5].The economic and efficient way for hydraulic fracturing to be performed is to create large-scale complex fracture networks, the in situ stress coefficient, strength and deformation properties of the rock, rock anisotropy and heterogeneity, distribution of natural fractures, and brittleness of the rock.Therefore, the mechanical properties of shale are an important part of shale gas reservoir evaluation [8][9][10].The mechanical properties of the rock matrix and bedding development, as well as the hardness and brittleness of the rock, can lead the borehole wall of shale gas to collapse, causing leakage, and other instability problems [11,12].Understanding the anisotropic behaviors of shale rock has important impacts for shale energy exploration, wellbore stability, the interpretation of micro-seismic monitoring, etc. [13].Shale has many fine pores and micro-cracks due to its fine mineral particles and structure.The developmental characteristics of these micro-pores and micro-cracks will have certain impacts on the mechanical characteristics of shale at the macro level.The analysis of the internal microstructure of shale specimens can provide some auxiliary reference value for the subsequent test analysis.Therefore, it is of practical significance to investigate the anisotropic mechanical properties of shale, for the purpose of better application to the hydraulic fracturing of the shale reservoir.
As a highly differentiated stratified rock, shale has strong anisotropy, and is generally considered a transversely isotropic material [14][15][16][17].The anisotropy of shale rock can be mainly attributed to certain factors, such as mineral orientation [18][19][20], microcracks [21], stress state [22,23], and kerogen content [24].Overall, it is important to conduct comprehensive research on the shale's characteristics while considering its heterogeneity and anisotropy [25].Scholars have summarized the effects of various minerals on the mechanical characteristics of shale [26,27].These studies show that clay and a TOC content of 30-40 wt% are critical for micro-structural deformation [28], ensuring the various flow properties [29] and elastic properties [30], and revealing that mineralogy controls the mechanical response of the rock matrix.Therefore, to evaluate a shale reservoir as to whether it is a good candidate for fracture stimulation, rock mineralogy analysis is necessary, and a result indicating high levels of quartz and a low clay content is quite positive [31].
Shale anisotropy plays a significant role in engineering activities, and weak beddings have important impacts for the initiation, propagation, and formation of the hydraulic fracture networks of shale [32].In several different studies, a distinct anisotropy of shale rock has been proven to have had a distinct effect on mechanical behaviors under different stress conditions [25,27,[32][33][34][35][36][37], as well as predominate in the initiation and propagation of hydraulic fractures [38][39][40].Many previous experimental investigations have studied the mechanical behaviors of different clayey rocks.In the prior experimental investigations, the shale was processed in different bedding inclinations.Many laboratory-based mechanical tests on shale specimens from outcrops and reservoirs have been conducted, including compressive tests [32,[41][42][43], tensile tests [44], three-point bending tests [45], direct shear tests [36,46], etc.These studies provided adequate data and grounds for comprehending the anisotropic mechanical properties of the shale specimens.
Brittleness plays a critical role in hydraulic fracturing design [47].However, brittleness cannot be measured directly from seismic data, or well-log data, either [48].Using the mineralogy and mechanical parameters, various brittleness indices have been developed to evaluate the fracability of shale formations.A brittleness index based on mechanical parameters may provide a more accurate result, as these parameters can reflect the influence of stress on the brittleness [2], since the brittleness changes as the shale-formation depth increases [49].An index based on mineralogy has been widely employed in the oil and gas industry because the input parameters are easy to obtain; however, due to the high cost of well exploration, the required data are difficult to obtain [50].The reliability of the brittleness evaluation of shale should be improved.However, many of the brittleness indices function mostly according to the assumptions of isotropy theories [51], and it is well known that the anisotropy of shale rock rises to a high degree, thus, the anisotropic properties should be considered in the brittleness evaluation.
Moreover, the previous tests on shale mechanics [52][53][54] show that the anisotropy of the fracture toughness of shale is high, and these results show scatter between different shales.During the shale fracturing, the beddings are usually weaker than the rock matrix, and this can easily lead to crack propagation, i.e., when the fracture meets a bedding plane it can either extend along it or penetrate across it [36].It is established that when a hydraulic fracture extends in the shale rock, it will be positively influenced by the bedding planes.The complex interaction between hydraulic fracture and bedding planes could be affected by several factors, i.e., the bedding strength, inclination, fluid injection rate, and so on.Compared to other discontinuities, bedding planes with large scales can be found more frequently in shale, particularly in continental stratum [55].In general, the bedding planes have distinct effects on several aspects, including weakening the shear and tensile strength, disturbing fracture propagation, and complicating the interaction between the existing and induced fractures [56][57][58][59][60].These results revealed that the stimulated fracture networks and the well extraction process are significantly affected by these discontinuities.Therefore, it is necessary to study the relationships between the bedding plane and the hydraulic fracture characteristics [61].
This paper briefly reviews previous work on the anisotropic mechanical properties of shale rock and their relation to hydraulic fracturing in shale reservoirs.The lithological characteristics of shale rock, particularly the mineral constituent, determine its physical and mechanical behavior in essence.In this paper, the status of the research addressing the lithological characteristics of shale rock is summarized first, and then the anisotropic physical and mechanical properties of shale specimens, including ultrasonic anisotropy, mechanical behavior under compression tests, and tensile property under the Brazilian test, as well as fracture toughness anisotropy in the literature, are summarized.The topics associated with anisotropic mechanical behavior under laboratory tests are addressed in this paper; in particular, the evaluation of shale brittleness based on its mechanical characteristics is discussed in detail.Finally, further concerns, such as the effects of the bedding on hydraulic fracturing failure strength, crack propagation, and failure pattern, are also drawn out.This review study will provide a better understanding of the current research findings on the anisotropic mechanical properties of shale rock, which might provide insight into the aspects of shale anisotropy related to hydraulic fracture propagation in shale reservoirs.

The Lithological Characteristics of Shale Rock
Rock mineralogy is essential for the evaluation of shale reservoirs [62,63], and can greatly determine the brittleness [47].The anisotropy is induced by the packing density, which reveals that, in situ, the anisotropy of shale is from the deposition of the clay particles rather than from the intrinsic anisotropy of the mineral.
The mineral components of shale vary with different depositional and burial environments and aspects of its mineralogy [64].The basic mineral components of shale rock are carbonates, quartz, clays, and feldspar.In many research works, the shales can be named according to the mineral enrichment, such as clay-rich shale [65][66][67], organic-rich shale [68][69][70], and carbonate-rich shale [71,72].The relevant mineralogy plays an essential role in hydraulic fracturing implementation.Different proportions of quartz, carbonates, and clays could result in different mechanical properties in shale rocks.Shales with abundant quartz and a low clay content usually are of high brittleness, and shales that have a high clay content and low quartz content are low in brittleness [47,73].
As for the clay-rich shale, the common intrinsic clay minerals in petroleum reservoirs are kaolinite, chlorite, illite, smectite, and mixed-layer clays [74].For instance, Pierre Shale [65], is a clay-rich rock, with its mineralogy dominated by mixed-layer illite-smectite, discrete illite, and quartz.Its microfabric shows a distinct amount of fine-grained rigid minerals floating in a clay matrix.This kind of shale is layered in places, with some high aspect-ratio rigid minerals, organic matter, and clays aligned sub-parallel to the bedding.Clay-rich shales have better reservoir properties, and they show higher porosity, pore volume, and specific surface.Clay-poor shales have higher concentrations of brittle minerals, resulting in limited shale gas storage, which makes the shale inadequate for shale gas exploration [66].Organic-rich shales are usually considered to be good for petroleum and gas extraction [75].The micro-pores are often developed at the nano-level in the organicrich shale, leading to a prominent capillary effect and a high spontaneous imbibition effect [69].As for carbonate-rich shale, such as Mancos Shale [71], as shown in the XRF maps shown in Figure 1, the light band from the petrographic data indicates the quartz-calciterich laminae.It can be seen that Si is distributed in the matrix, while Ca, Al, Fe, K, and Mg are distributed between the optically lighter and darker layers.It also can be found that Al, K, Fe, and Mg are intensively distributed in the more fine-grained (optically darker) layers, and Ca is intensively distributed in the more coarse-grained (optically lighter) layers.The micro-texture of the quartz-calcite-rich laminae consists of isolated organic material and coarser-grained quartz, dolomite, and feldspar minerals, cemented by calcite.High quartz contents and low clay contents contribute to high levels of brittleness and help hydraulic fracture propagation.Li et al. [76] conducted research on black siliceous shale specimens from the southeastern edge of the Sichuan Basin of the Upper Yangtze Plate.The X-ray diffraction (XRD) results show that the brittle minerals of the shale specimen are 60-80 wt%, and that this gradually increases as the burial depth increases (Figure 2a). Figure 2b also shows that the clay content is low, and that the shale in the upper has clay-rich minerals, while the shale in the lower has brittleness-rich minerals like quartz [50].

The Anisotropic Mechanical Behavior of Shale Specimens under Different Loading Conditions
The anisotropy should be carefully considered when attempting to predict the mechanical behaviors of shale [77].Many research efforts have been conducted to study the mechanical behaviors of different rock materials containing high clay contents.As for those experiments, the shale specimen was processed in various loading directions relative to the bedding planes.In those studies, it was found that the anisotropic inclination of a shale specimen under a loading test is generally defined by two aspects: one is the angle of the loading line to the bedding plane [77][78][79], and the other one is the angle of the loading line relative to the normal direction of the bedding plane [11,25,80].These two angles are complementary.To better compare and describe the test results of shale specimens in different studies, the anisotropic inclinations must be unified.As depicted in Figure 3, the angle θ of the loading line to the normal direction of the bedding plane is employed in this paper [2,11,81].As illustrated in Figure 3a, the cylinder specimens with different bedding inclinations were loaded under the compression test, while, as seen in Figure 3b, the disc specimens were loaded under the Brazilian split test.Therefore, the anisotropic inclination θ is also the angle of the bedding plane relative to the horizontal line.The test data from different studies described in the following text are all standardized to this anisotropic inclination θ.

Ultrasonic Anisotropy of Shale Specimens
Ultrasonic techniques are maneuverable and nondestructive, and therefore useful for applications in laboratory conditions, and are promising for geophysical applications, including the interpretation of the dynamic mechanical behaviors of shale [77].The P-wave travels faster than the S-wave.The anisotropy affecting wave velocity in the shale subsurface is considered to be a central difficulty, one resulting in significant issues for understanding dynamic elasticity in geophysical interpretation.Elastic constants of geologic materials, such as shale reservoirs, can be determined by measuring wave velocities [82].Lo et al. [83] conducted a series of tests on the wave velocities associated with Chicopee shale by using ultrasonic techniques.Their experimental results show that the mineral orientation and the micro pores or cracks are the main factors resulting in elastic anisotropy.Additionally, the elastic anisotropy decreases with increasing confining pressure.Allan et al. [84] investigated the anisotropy of wave velocity in an organic-rich shale by proposing a multiscale methodology as well as various analyses based on XRD tests, BSE-SEM tests, and X-ray microtomography, which were applied to determine and understand the sources of wave velocity anisotropy.
Based on some ultrasonic test data from shale specimens [77,85,86], Figure 4 depicts the relationships between ultrasonic velocities and bedding angles.It can be clearly seen that the ultrasonic velocities are obviously affected by the bedding planes, and that the velocities of P-and S-waves parallel to the bedding plane are greater than those of the waves orthogonal to the bedding plane.The energy decreases due to the wave scattering on different beddings, leading to increases in the time duration of the wave traveling through the specimens.The above conclusions are similar to the results from Dewhurst et al. [87] and Zhubayev et al. [88].The experimental results from Mokhtari et al. [44] show that the P-wave velocity of a shale specimen with a bedding inclination of 0° is 30% lower than a specimen with a bedding inclination of 90°.Based on the propagation and polarization directions relative to the bedding plane, six ultrasonic wave velocities can be measured in shale specimens [89].Figure 5 illustrates the six velocities, as follows [90]: one compressional wave (Vpv) and two shear waves (Vsv1 = Vsv2), which propagate normally to the bedding plane in a vertical plug (as shown in Figure 5a, θ = 0°); and one compressional wave (Vph) and two shear waves (Vsh > Vsv), which propagate parallel to the bedding plane in a horizontal plug (as shown in Figure 5b, θ = 90°).These six ultrasonic velocity parameters could then be employed to calculate the elastic parameters, which are used for describing the anisotropic elastic properties of a vertically transverse isotropic (VTI) material like shale rock [91,92].In the study conducted by Yin [93], the wave's attenuation by bedding planes in the transversely isotropic rock is also described.As depicted in Figure 6, when the ultrasonic wave travels normal to the bedding, the transmission and reflection occur synchronously at the bedding plane (Figure 6a), and the wave's reflection results in energy loss and reduces the wave's velocity; with more bedding planes orthogonal to the propagation direction, greater energy losses occur, leading to a lower velocity of the ultrasonic wave.At a bedding inclination of 0° (Figure 6a), the bedding planes are orthogonal to the propagation direction of the P-wave, resulting in strong reflection at the bedding planes and a lower wave-velocity.With the increases of bedding inclination, the beddings are slowly brought parallel to the propagation direction, and the energy attenuation becomes less, leading to an increase in wave velocity.When at a bedding inclination of 90° (Figure 6b), the bedding planes are parallel to the propagation direction of the P-wave, and there is almost no reflection and the energy attenuation is minimum; thus, the wave velocity is at its maximum.

Uniaxial and Triaxial Compression Test
Shale specimens with various bedding inclinations are frequently studied to investigate their anisotropic behaviors.Many researchers have conducted experiments on various shale specimens under compression conditions [36,41,42,46,78,[94][95][96][97][98][99][100][101][102][103][104].The above studies reported the mechanical properties of anisotropic rocks; in particular, in Ramamurthy's research [105], the anisotropy of rocks is generally classified into three modes: "U" type, undulatory type, and "shoulder" type.These are are shown in Figure 7. Strength under triaxial compression conditions is the key characteristic of rock material.Some tests are conducted under uniaxial compression [2,11,25,36,[77][78][79][80]85,86,95,106,107], and some tests are conducted under triaxial compression conditions [2,36,78,85,108].The various shale specimens taken from different sampling sites in these studies display different magnitudes of strength, but their values all show a similar law of variation concerning the bedding inclination.Based on the test data from the abovementioned studies, the uniaxial compressive strength (UCS) and triaxial compressive strength (TCS) variations with respect to bedding inclination are displayed in Figures 8  and 9, respectively.From Figures 8 and 9, it can be concluded that the peak strength instances all display a "U" type or a similar "shoulder" type variation with increased bedding inclination under different confining pressures.Particularly, the maximum peak strength is observed at θ = 0° or 90°, and the minimum values are observed at 45° < θ < 75°.To better predict the strength of anisotropic rocks, many strength criteria have been proposed.The earliest attempt was the theory of a single weakness plane [109], and this was followed by various subsequent studies [110][111][112][113][114][115][116][117]; in one, an empirical strength criterion was proposed for predicting the non-linear strength behavior of intact anisotropic rocks [111,112]: where σ1 and σ3 are the major and minor principal stresses, respectively, and σcj is the UCS at the particular θ.The values for αj and Bj are defined herein to consider the anisotropy of the strength, and the functions are as follows: where σc90 is the UCS at θ = 90°, and α90 and B90 are the values of αj and Bj at θ = 90°.The values for αj and Bj at each anisotropy orientation can be calculated by substituting the triaxial compression test data into Equation (1).Additionally, another empirical criterion under triaxial conditions can be expressed as follows [118]: where σci is the UCS of intact rock; A and B are constant parameters; and r is the strength reduction factor, with a value equal to 0 for intact rock and 1 for highly jointed rock masses.To better apply the above failure criterion to transversely isotropic rocks, Saeidi et al. [114] modified the failure criterion (Equation ( 4)) to the following: where σcθ is the UCS of anisotropic rock at θ, and α is the strength reduction parameter when taking into account the anisotropy, which is considered here in order to extend the generalization of Equation ( 4) to anisotropic rocks.
In addition, several mechanical models based on a damage theory have been proposed for anisotropic rocks by a series of researchers [119][120][121][122][123].These models show a good capability for describing the main mechanical behaviors of anisotropic rocks, but still need to be verified by more experimental data for different types of anisotropic rocks.To evaluate the effect of the confining pressure on anisotropy, the strength anisotropy degree Ra, represented in terms of σpmax/σpmin, can be calculated in each confining pressure [113,124].The strength anisotropy degree Ra indicates the maximum degree of strength deviation under different bedding inclinations.Figure 10 shows the variations of Ra with increases in confining pressure, the data for which are calculated based on the studies in Figures 8  and 9.It can be seen that the strength anisotropy degree Ra of a shale specimen decreases as the confining pressure increases.In other words, the confining pressure can reduce the strength anisotropy of rock material.The reason is that the confining pressure is hydrostatic pressure, and when it is applied to the specimen, the beddings are compacted, and thus the bedding strength is enhanced, and the bedding effect is weakened.In addition, Figure 11 shows the elastic modulus with bedding inclination increases under uniaxial compression of shale specimens, based on data drawn from some previous studies [11,25,36,77,80,86,106,107].These studies all reported an increasing variation with the bedding inclination increases overall, even though there were some local fluctuations.The explanation for this can be developed as follows: (1) the directional alignment of the clay minerals and micro-cracks in the bedding plane is the inherent reason for the strong elastic anisotropy of the shale specimen [125], and the content of the clay and the development of micro-cracks are different in the specimens under different bedding inclinations; this results in the anisotropy of the elastic modulus [86]; (2) as typical sedimentary rock, the compaction degree of the bedding in shale rock is low in the diagenetic process.Therefore, when the principal stress is parallel to the bedding (θ = 90°), the bedding is under a tensile state and the rock matrix is under compression; thus, the elastic modulus and deformation modulus are large and, this is in accord with the rock matrix in this case.When the principal stress is perpendicular to the bedding (θ = 0°), the bedding is in a state of compaction, and the axial strain is larger, resulting in a smaller elastic modulus and deformation modulus.As for the failure behavior of shale rock under the conventional triaxial compression test, the fracture propagation during the failure process is generally affected by the bedding inclination and the confining pressure.The fracture angle strongly depends on the bedding inclination and the confining pressure.In the summary of Niandou et al. [78], for a bedding inclination of 0° < θ < 30°, tensile fracture along the loading direction occurs under low confining pressure and shear fracture across the bedding plane occurs under high confining pressure.For a bedding inclination of 30° < θ < 75°, the failure generally occurs by the bedding sliding; however, shear fracture in the matrix can occur and may cross the bedding plane when under a high confining pressure.For a bedding inclination of 75° < θ < 90°, when under low confining pressure, the failure occurs by the bedding plane splitting, and failure occurs by a shearing across the bedding plane when under a high confining pressure.In the study of Yang et al. [2], the researchers conducted conventional triaxial compression tests on Longmaxi Formation shale specimens under different bedding inclinations.According to the test results, they classified the shale specimens into four failure modes, as illustrated in Figure 12 and Table 1.Four failure modes are reported: tensile fracture through the bedding plane (T-T), tensile fracture along the bedding plane (T-A), shearing through the bedding plane (S-T), and shearing along the bedding plane (S-A).T-T generally occurs under uniaxial compression (or very low confining pressure) with low and medium bedding inclination (0°-45°).S-A mainly takes place in specimens with medium or high bedding inclinations (45°-75°).With increases in the confining pressure, S-T occurs in specimens with low or medium bedding inclinations (0°-45°), instead of T-T, and even in some high bedding inclinations, the S-T dominates the failure.As for the bedding inclination of 90°, T-A occurs when under uniaxial compression, as the bedding is activated by tension along the loading direction.With the increases in confining increases, S-T takes place in the rock matrix.It can be concluded that the bedding planes at medium and high inclinations have a distinct effect on the failure fracture, resulting in shear slipping along the bedding plane; with the increases in confining pressure, high confining pressure dominates the failure mode and shear fracture in the matrix (across the bedding plane) generally occurs, which is similar to the mode of isotropic rock material.These findings are also in accord with the studies of Liu et al. [80] and Zhao et al. [126].Table 1.Failure mode summary for shale specimens, from the study of Yang et al. [2].

Tensile Properties under Brazilian Split Test
The tensile strength of shale plays an important role in the initiation and propagation of hydraulic fractures.As a method classically used to evaluate the tensile strength, the Brazilian test with different bedding inclinations is widely employed for determining structural anisotropy.The analytical solution developed to measure the tensile strength of anisotropic rocks was developed by Amadei et al. [127] according to the theoretical relation between the stress and strain proposed by Lekhnitskii et al. [128].The effects of bedding inclination and thickness on tensile strength and the fracture pattern of anisotropic rocks under diametrical loading conditions have been investigated by conducting various Brazilian tests [15,[129][130][131][132][133][134].
Vervoort et al. [133] summarized the failure anisotropy of nine different anisotropic rocks under Brazilian test conditions into four trends, as follows (as shown in Figure 13).Trend 1: Constant value over the entire range of anisotropy angles (Figure 13a).Trend 2: Constant value between 0°and 45°, followed by a linear decrease (Figure 13b).
Trend 3: Decrease in the strength over the entire interval, but a rather systematic decrease, approximating a linear variation (Figure 13c).
The above four variation trends describe different types of anisotropic rocks with varying bedding planes.As for the shale rock specimens from different studies, the BTS variation with respect to bedding inclination is similar.Figure 14 summarizes some BTS data of shale specimens from different studies.From their test results, depicted in Figure 14, it can be seen that the BTS variation of shale specimens is more like the above trend 2 and trend 4. The studies of Jin et al. [25], Yang et al. [81], and Cho et al. [95] show that the BTS value is constant between 0°and 60° and then followed by a linear decrease, which is in accord with trend 2. The other studies, such as Yang et al. [135], Hou et al. [33,136], and Du et al. [137], show a trend similar to trend 4; their BTS decreases first to a low value (the anisotropy angle is around 60°) and then keeps on, unchanged.The BTS from the study of He and Afolagboye [138] shows the variation of trend 3, in which BTS decreases over the entire interval.From the above summary, it can be seen that, although the BTS variation of shale specimens from different sampling sites shows a slight difference, in general, their BTS values do not demonstrate much change from 0° to 30°, and then decrease to a low level around 45°-60°, and finally keep within a stable range around 75°-90°.To further explain the failure behavior of the shale disc specimen with respect to the bedding plane, the failure pattern in the Brazilian test should be discussed.In the study of Hou et al. [136], shale disc specimens are loaded under Brazilian tests with different bedding inclinations.They classify the Brazilian split failure process of shale specimens into three stages: compaction, elastic deformation, and destruction stages.Their test results show that the anisotropic characteristics of tensile strength and failure pattern are highly distinct at the peak point of stress.The bedding plane has a strong effect at the bedding inclination of 60° and a weak effect at 0°.They also reported the transformation angles of the failure mechanisms, which are the bedding inclinations of 60° and 45°.At 60°, a tension fracture along the bedding turns to cross the bedding plane, and a crossbedding plane fracture turns to shear slip at 45°.These findings are similar to those in the studies of Cho et al. [95] on Boryeong shale, Vervoort et al. [133] on Freiberger Gneiss, Du et al. [137] on Shaanxi shale, and Yin and Yang [104] on layered sandstone.Cho et al. [95] investigated the BTS of Boryeong shale under different loading angles.The study's results show that the maximum tensile strength occurs near 15°, and the minimum value occurs when the bedding line is parallel to the axial direction (90°).Additionally, the shale begins to fracture along the bedding when θ > 30°.When 0° ≤ θ ≤ 15°, the crack propagates along the loading direction.
From the above studies and discussion on the anisotropic characteristics of BTS and failure patterns, it can be concluded that, with the bedding closely aligned to the loading direction, the bedding's effect on the BTS turns out to be significant.In other words, there seems to be a critical loading angle; the BTS could not be affected before the critical angle, and decreases rapidly when it exceeds it.In the research of Yin and Yang [104], many disc specimens of layered sandstone were tested and a critical bedding inclination θ of 71.6° was derived by theoretical analysis.In their studies, to better explain the stress state at the center, the normal and shear stress coefficients kN and kS are defined.As shown in Figure 15, when 0° < θ < 71.6°, the compressive stress and shear stress are applied on the bedding plane, and the disc specimen can evince tension in the rock matrix and shearing at the bedding; and when 71.6° < θ < 90°, the bedding plane is under tensile and shear states simultaneously, and when the shear stress is at a high level, the disc specimen can be shear or tensile along the bedding plane.The bedding may be under compression, shearing, or tension states with respect to the variation of bedding inclination.When the bedding inclination is low, the bedding is dominated by compression, does not easy fail.When at a high bedding inclination, the bedding is under shearing and tension simultaneously, and easily experiences shear slip and tensile split.Figure 16 presents some failure examples for shale disk specimens under the Brazilian test, with different bedding inclinations from different studies shown.From these fractured disk specimens with different loading directions, it can be seen that the fracture pattern of the anisotropic rock is much different from those of isotropic rock.Although the fractures generally propagate along the loading direction, the bedding directions have distinct effects on the fracture propagation to a certain extent.In some specimens with high bedding inclinations of 60° and 75°, such as the Mancos and Boryeong Shale with 60° (Figure 16a from Simpson et al. [139] and Figure 16b from Cho et al. [95]), Chongqing shale with 75° (Figure 16c from Hou et al. [136] and Figure 16d from Wang et al. [32]), the fractures extend along or divert to the bedding planes.The principal reason for this pattern is that the bedding plane is very closely aligned to the loading direction when under a high bedding inclination, and the shear stress or tensile stress is dominant in the bedding plane and this easily leads to shear or tensile fracture along it.This is in accord with the study and explanation of Yin and Yang [104] in Figure 15.By conducting a series of Brazilian tests on layered sandstone, Tavallali and Vervoort [131] considered the specimens after failure and classified the observed fractures into three different types (see Figure 17): (1) layer activation (LA), in this pattern, the fractures extend along the bedding planes; (2) central fractures (CF), in this pattern, fractures are roughly located in the central part of the specimen and parallel to the loading direction; (3) noncentral fractures, in these cases, the fractures are outside the central part.As shown in Figure 17, the mode of LA mostly occurs at a high loading direction of 60°-90°, and the mode of CF mostly occurs at a low loading direction of 0°-30°, while the combination of LA and CF appears at a loading direction of 45°-75°.The failure modes at θ = 0° and 90° look similar; these are all central fractures along the loading direction.However, their fracture mechanism is different: the central fracture in specimen θ = 0° extends across the beddings, which shows the rock matrix fractured by tension; and the central fracture in specimen θ = 90° extends at the beddings, which shows the bedding fractured by tension and is more in accord with the layer activation.This is also can be explained by Figure 15: the BTS of θ = 0° is rather similar to the BTS of the rock matrix, and the BTS of θ = 90° corresponds least to the BTS of the bedding plane [131].

Fracture Toughness Anisotropy of Shale Specimens
The fracture toughness of mode I (KIc), is an important parameter that controls the hydraulic fracture propagation [140,141].To determine the KIc of brittle rock material, many experimental methods and configurations have been designed; some common methods and their configurations are listed in Table 2.As described in Table 2, the configurations of CCCD, SCB, and SENB are specimens with straight-through notches, and the others are specimens with chevron notches.To make sure that fracture behavior can be characterized appropriately by the linear elastic fracture mechanics (LEFM), large specimens are selected, such as deep beam of NBD and CNDB.

Methods and Specimen Configuration Graphic Illustration References
Centrally cracked circular disc (CCCD) Awaji and Sato [142], Atkinson et al. [143] Semicircular bend specimens (SCB) Chong and Kuruppu [144], Lim et al. [145], Dai et al. [146], Funatsu et al. [147], Ren et al. [148] Single-edged notched beam specimens (SENB) ASTM [149] Straight edge cracked round bar bending specimens (SECRBB) Bush [150] Single edge-notched deep beam specimens (SENDB) Luo et al. [37,151] Chevron-notched beam specimens (CNB) Wu [152] Chevron-notched deep beam specimens (CNDB) Ren et al. [153] Cracked chevron-notched Brazilian disc specimens (CCNBD) Sheity et al. [154], Fowell [155], Dai et al. [156] Chevron-notched semicircular bend specimens (CNSCB) Kuruppu [157], Dai et al. [158] Short-rod specimens (SR) Barker [159], Ouchterlony [160] Traditional Brazilian disk and flatted disk specimens without pre-existing flaws Wang and Xing [161] To determine the fracture toughness anisotropy of the shale specimens, three main crack orientations can be designed [162], as shown in Figure 18: arrester, divider, and short-transverse.By conducting loading tests on these three kinds of specimens using the different methods listed in Table 2, the fracture toughness tested in different crack orientations shows distinct differences.The fracture toughness of specimens with cracks oriented parallel to the bedding plane (short-transverse, Figure 18c) is smaller than in the other two crack orientations [54,151,153,162,163].This result is easy to understand, and the reason for it is that the specimens with a pre-existing crack parallel to the bedding planes are easy to initiate, due to the weak tensile strength of the bedding plane [164].From Figure 18, we can see that in the case of Arrester (Figure 18a), the crack extends through the isotropic plane (bedding plane).In the case of Divider (Figure 18b), the crack also propagates in the rock matrix but does not need to extend through the bedding plane.Thus, comparing the cases of Arrester and Divider, the crack propagation is slightly harder in Arrester, as the crack needs to expend more energy to penetrate through the isotropic plane.This is clearly illustrated in Table 3, which summarizes experimental results from various studies conducted using different testing methods.From Table 3, it can be found that almost all the fracture toughness values of KIc, A are greater than KIc, D, while, meanwhile, their ratio is nearly always around 1, which indicates that the crack propagation in Arrester and Divider is dominated by the rock matrix.On the other hand, it can be clearly seen that in the case of Short-Transverse (Figure 18c) the crack initiation and propagation are mainly needed to open the beddings, and this expends less energy than used in the other two cases due to the weak tensile strength of bedding plane.Thus, as displayed in Table 3, the value KIc, ST is much lower than the KIc, A and KIc, D. As for other anisotropic directions of crack propagation, Lei et al. [169] have conducted investigations regarding the fracture behavior of shale samples with different bedding strengths.Their beam specimens are processed with different bedding inclinations according to the straight-through notch and tested under three-point bending loading.The micro-crack propagation process was observed by the SEM equipment with a loading system.The results revealed that the outcrop shale of the Longmaxi Formation from Pengshui County, in the southeast of Chongqing, China, shows strong anisotropy and that this has a great effect on hydraulic fracture propagation.In Figure 19a, it can be observed that the micro-cracks initiate at the notch tip, and then propagate along the loading direction.Figure 19b shows the simulated failure patterns of three-point bending specimens under different bedding properties.The tensile strength of the smooth joint decreases from mode Ⅰ to mode Ⅳ.The results revealed that the fracture characteristic of shale is affected by bedding inclination and bedding strength.It is also can be determined that the cracks exhibit strong tortuosity under a low tensile strength of the bedding plane.On the other hand, as shown in Figure 20a, the study results of Lei et al. [169] reveal that, when under the same bedding strength, the fracture toughness decreases with the increases of bedding inclination, with a maximum at 0° and a minimum at 90°.Particularly, the fracture toughness decreases gently between 0° and 45°, and then decreases rapidly after 45°.It indicates that when the bedding inclination approaches the loading direction, the bedding's effect on fracture toughness turns to be more obvious.It can be clearly seen from Figure 19b that the crack propagation eventually occurs at the bedding planes more obviously when the bedding inclination increases after 45°.In particular, as displayed in Figure 19b, mode Ⅳ has the lowest bedding strength, resulting in the most obvious crack deflection towards the bedding plane.Luo et al. [37] and Shi et al. [170] have investigated fracture toughness anisotropy under different bedding inclinations.The shale specimens used by Luo et al. [37] and Shi et al. [170] are all sampled from the Longmaxi Formation.The former study processed NDB shale specimens, with a lengthto-width ratio (L/W) of 2.0 and an edge crack with inclined angle  = 0° and bedding inclination θ = 0°, 30°, 60°, and 90°.The latter study employed the CCNBD specimens with chevron notches with angles of 0°, 30°, 45°, 60°, and 90°. Figure 20b depicts the KIc variations concerning the shale bedding inclination from the above two studies.The results show that the bedding inclination has a distinct effect on the fracture toughness, and the tested values all decrease with an increase of bedding inclination.
In the above studies, the test specimens with bedding inclinations of 0° and 90° are configured in accord with the cases of Arrester (shown in Figure 18a) and Short-Transverse (shown in Figure 18c), respectively.It can be determined from Figure 20 that the maximum value is at 0° and the minimum at 90°, which is in accord with the findings from Table 3.The crack's approach angle relative to the bedding plane and bedding strength are the key influence factors in hydraulic fracture propagation.The high approach angle and low bedding strength contribute to the fracture propagation along bedding planes.The Short-Transverse specimen has the lowest fracture toughness, due to the weaker bedding plane, and has a looser grain arrangement and is less resistant to crack propagation [171].The microcracks could extend and coalesce, resulting in the failure of grain clusters [172].

Quantitative Evaluation of Shale Brittleness Based on Mechanical Properties
Brittleness is commonly applied to characterize the failure behavior of the rock material [173].For hydraulic fracturing, the high brittleness of shale may contribute to the deeper and more complex fracture propagation.Therefore, the evaluation of brittleness is important in engineering.Based on various previous studies, brittleness can be defined as the lack of ductility [174,175], as a failure at or only slightly exceeding the yield stress [176], as the cohesion destruction [177], as a failure with small or no plastic deformation [178], and as a self-sustaining failure process [179].
As for hydraulic fracturing in shale, the brittleness index (BI) is usually employed to quantify the brittleness degree of rock, which is useful for evaluating the ability of fracture propagation (fracability) and fracture network generation.To accurately measure the brittleness of shale rock, many factors should be considered.To date, many BIs have been proposed to quantify its extent [73,[179][180][181][182][183][184][185][186][187].Table 4 lists some common BIs, mainly based on mineral composition, and strength parameters.
Among the test methods in Table 4, mineral composition is easy to obtain and can be precisely determined by conducting laboratory analyses such as XRD testing.Methods based on strength parameters, stress-strain characteristics, and energy balance analysis can be determined by stress-strain curves.As the stress-strain curves can be obtained easily by conducting the triaxial compression test, this measurement is commonly used to determine strength parameters.The brittleness is a mechanical response to a specific stress state.For a material, the composition is constant, but the brittleness can vary if different stress states are applied.Therefore, the specific stress state must be considered when its brittleness is evaluated.Thus, BI definitions based on stress conditions, such as strength parameters, stress-strain characteristics, and energy balance analysis, are more significant [2].As for the brittleness indices BI5~BI9, which are based on strength parameters, parameters such as UCS and BTS are easily obtained by a laboratory test.All of these BIs display a positive relationship [179] with confining pressure.On the other hand, rocks that show different stress-strain curve shapes may be of the same BI values (BI11~BI13).Furthermore, the combination of strength and strain performances (BI10 and BI13) is more precise in predicting rock brittleness [173].Additionally, some BIs based on stress-strain (BI11~BI13) and energy balance (BI17) only consider the effect of the pre-peak state and ignore the stress drop rate; however, some rock materials show a distinct ductility when high confining pressures are applied [185,195].
The drilling of horizontal wells for hydraulic fracturing in shale usually needs to cross bedding planes in different orientations; thus, the shale anisotropy should not be ignored in a brittleness evaluation [50].According to the various studies on anisotropic mechanical properties of shale specimens from Sections 3.2 and 3.3, the UCS shows a "U" type variation and BTS shows a negative relationship with bedding inclination.Thus, as illustrated in Figure 21, from the equations of the BI5 and BI6 in Table 4, which are calculated by UCS and BTS, it can be concluded that the variation of the two BIs concerning bedding inclination will display a similar "U" type.Before the inflection point, the UCS and BTS decrease, and the UCS decreases faster than the BTS; therefore, the BI5 and BI6 decrease gently.When the bedding inclination increases to a high level, such as 90°, the BTS decreases to a very low level and UCS is at a high level; thus, the BI5 and BI6 increase to high magnitudes rapidly.It should be noted that the BI5 and BI6 variations before the inflection point closely depend on the margins of UCS and BTS reductions; as shown in Figures 8 and 14, the UCS and BTS reduction from various shale specimens are different.Therefore, in Figure 21 the BI5 and BI6 variation before the inflection point may show a decrease or a gentle increase, or, alternatively, it can be stable.In the study on shale brittleness evaluation from Yang et al. [2], the researchers obtained data from the triaxial compression test.They proposed two new brittleness indices, BI1 * and BI2 * , to evaluate the shale brittleness variation with respect to bedding inclinations.By comparing the two new brittleness indices with the BI from Kivi et al. [49], which is referred to herein as BI25, the BI variation concerning bedding inclinations is depicted in Figure 22.The new brittleness indices BI * 1 and BI * 2 from Yang et al. [2] and BI25 from Kivi et al. [49] are expressed as follows: * =  /( +  ) = 0.5( / +  /( +  )) where σf is the strength and σr is the residual strength, respectively, and M is the modulus after the peak.Uei is the ideal elastic energy, as defined in Yang's study.[2].The other parameters are noted in Table 4.
From Figure 22a,c, we can see that the BI1* and BI25 display very similar variations, roughly increasing with respect to the bedding inclination; in particular, they increase rapidly when the bedding inclination is at 90°.This indicates that the shale specimen shows higher brittleness when the bedding inclination is at a high level.This experimental finding is in accordance with the conclusion stated in Figure 21.In the study of Qian et al. [50], a novel BI is proposed for analyzing the anisotropic characteristics of brittleness [196,197].Their theoretical research suggests that Young's modulus in the vertical is lower than that in the horizontal, and Poisson's ratio in the vertical can be either higher or lower than that in the horizontal.Thus, the BIs, which are determined by parameters of Young's modulus and Poisson's ratio, can vary significantly when drilling in shale formations with strong anisotropy.
The brittleness, under uniaxial compressive conditions, which results in tensile failure is higher than the brittleness when confining pressure is applied.Therefore, it can be concluded that the shale specimens that failed by tension show higher brittleness than those that failed by shearing.The reason may be attributed to the shearing slip along the fracture plane having led to a more plastic deformation.Specimens with fractures extending along the bedding show higher brittleness than those samples with fractures that extend across the bedding.[2].Therefore, according to the conclusions from Figure 12 and Table 1, the order of brittleness can be as follows: T-A > T-T > S-A > S-T.

The Anisotropic Effect of Bedding Plane on Hydraulic Fracture Initiation and Propagation of Shale
The weak bedding has a key effect on hydraulic fracturing characteristics [61].Many researchers have attempted to reveal the relationship between the breakdown pressure and the bedding inclination.A general conclusion has been reported that the breakdown pressure is considered to be decreased with increases of bedding inclination [40,198].By comparing different experiments results, it is found that the breakdown pressure and BTS show similar variations, in which the maximum values and minimum values occur at the bedding inclinations of 0° and 90°, respectively [40].However, the breakdown pressure also can be increased first and then decreased relative to the bedding inclination [61].Figure 23 depicts the relationship of breakdown pressure and bedding inclination.As can be seen in Figure 23, the breakdown pressure first increases and then decreases, as the bedding inclination increases.The maximum and minimum breakdown pressures occur at 45° and 90°, respectively.This variation agrees well with the tendency determined in the research of Chong et al. [199], which synchronously considered the effect of confining pressure and under a condition of σv = 20 MPa, σc = 10 MPa, and Qinj = 6 mL/min.In the research of Lin et al. [40] and Zhang et al. [198], the breakdown pressure under triaxial hydraulic fracturing (σv = 25 MPa, σc = 20 MPa, Qinj = 12 mL/min) shows a decreasing tendency as the bedding inclination increases.The above two different conclusions reveal that the anisotropy of shale breakdown strength depends on the stress conditions, injection rates, and fluid viscosity [198,200].However, the breakdown pressure seems to fluctuate within a bedding inclination range from 30° to 60°, and the breakdown pressure, mostly at the bedding inclination of 90°, is minimum.The reason for this variation may be attributed to the effects of tensile strength and the critical stress intensity factor (critical SIF).A smaller bedding inclination contributes to more stress components on the bedding planes, enhancing bedding compression, and thus, higher breakdown pressure is needed to initiate a fracture.Additionally, according to some relative theoretical analysis [201,202], it has been revealed that when θ < 45°, the critical SIF increases as the bedding inclination increases, and when θ > 45°, the critical SIF turns to become smaller as the bedding inclination increases.Higher fluid pressure is needed for the fracture propagation when under higher critical SIF, and then breakdown pressure turns to be at its maximum value at θ = 45°.There are significant interactions between a hydraulic fracture (HF) and joints, bedding planes (BP), and faults [38].The development of BP is one of the key factors affecting hydraulic fracturing in shale formation [39].The angle of the HF as it approaches the BP is generally considered to be an important parameter in determining the interaction between them [203][204][205][206].At a low angle (<30°), the HF propagates along the BP, and crosses the BP at a high angle (>60°), while at moderate angles (30°~60°) the fracture arrests [38,39].These weak BP would be activated during fracturing if the friction coefficient or cohesion is at a low level, which contributes to creating a complex fracture network [203,[205][206][207][208][209].Wang et al. [210] conducted experiments on shale specimens to investigate the fracturing mechanism and the effect of BP.The shale was marine black carbonaceous with obvious bedding structures, which lead to strong anisotropy.They summarized the fracturing failure patterns after hydraulic fracturing.As shown in Figure 24, the BP induce the HF propagating orientation to change.The HF is of a complex pattern controlled by bedding inclination and stress condition.From the test result, the failure patterns can be summarized as follows: (a) central planar fracture, which mostly occurs when θ = 0° and 90°; (b) deflected fracture, which mostly occurs when θ ≤ 45°; and (c) layer-activated fracture, which mostly occurs when θ > 45°.It should be noted that the central planar fracture in the specimen of θ = 90° is also a layer-activated fracture, which is induced by the vertical bedding planes.In general, two propagation modes can be classified when an HF encounters BP, namely crossing and opening.For further classification, as shown in Figure 25, an HF encountering BP can cause four cases: penetration, diversion, offset, and termination [206,211,212].When under a true triaxial stress condition, the offset mode is observed during the intersection with BP [213].Penetration (Figure 25a) refers to an HF that crosses the BP but does not change its propagation path.Diversion (Figure 25b) refers to a vertical HF that is deflected into the BP and is divided into two branches.Offset (Figure 25c) refers to an HF re-initiating and leaving a step-over at the BP.To explain these propagation modes, one possible reason is that the shear stress of the interface easily results in HFs crossing the barrier or extending along the BP [203].Additionally, the strength of the BP may not be the same at different distances, leading to differences in the degree to which the bedding opens [214], and in that case some HFs may cross the BP and some terminate at the BP.The HF propagation is slightly affected by the BP and the natural fracture (NF) when under large vertical and horizontal stress anisotropy conditions [212].For the HF interacting with multi-BP in laminated shale formation, complex fracture networks are of great practical value.The activation of BP is the optimum condition for forming a complex fracture network in the vertical plane [214].A possible fracture mode is summarized and shown in Figure 26.The depicted four fracture modes are very similar to the fracture types generalized in the study of Yin et al. [215], in which the complex fracture network is formed by primary cross-bedding fractures and secondary interbedding fractures.From the laboratory test results (depicted in Figure 27), it can be determined that interbedding fractures and cross-bedding fractures are generated.
The basic HF patterns could be summarized as follows [214] (Figure 28): (1) Mode I (Figure 28a), an HF initiates and propagates perpendicular to the BP; (2) Mode II (Figure 28b), an HF initiates and propagates parallel to the BP; (3) Mode III (Figure 28c), an HF initiates and propagates normal to the BP and a complex fracture network is induced by the weak BP; (4) Mode IV (Figure 28d), an HF initiates and propagates parallel to the BP, and diverts into another propagation path when it meets a cemented fracture; and (5) Mode V (Figure 28e), an HF initiates and propagates from a few NFs, and is diverted into another propagation path by BP.To form the above patterns, the BP strength is the key: In the cases of Mode I and Mode III, the cement strength of the BP should be high, but in Mode III some weak points should be presented.In the cases of Mode II and Mode IV, the cement strength of the BP should be low, but some strong cement BP should be presented in Mode IV.In the case of Mode V, weak NFs occur near the initiation point and the cement of the BP is low.The above summaries are also consistent with the simulation results of Yin et al. [215], which reveal the important effect of the mechanical properties of BP on the propagation of HF.Additionally, the BP failure was closely related to its cohesion and friction coefficient.If the BP had a low friction coefficient or low cohesion, shear failure would obviously increase.If the cementing strength of the BP at different distances is not the same, the degree of the bedding opening would differ [206,216].

Conclusions
The efficient implementation of hydraulic fracturing is the key issue in the exploration of shale gas.The existence of bedding structure results in distinct anisotropy of the shale rock formation.The anisotropy has an important influence on wellbore stability and hydraulic fracturing implements.This paper first briefly reviews the previous research works on the lithological characteristics and anisotropic mechanical behavior of shale specimens, and then a brief discussion on the relationship of shale anisotropic mechanical properties to hydraulic fracturing evaluation is conducted.The concerns of anisotropic mechanical behaviors under laboratory tests are addressed; in particular, the evaluation of shale brittleness based on mechanical characteristics is discussed in detail.Further concerns, such as the bedding's effect on the hydraulic fracturing failure strength, crack propagation, and fracture pattern are also drawn.Some main conclusions are drawn, as follows: (1) The rock mineralogy is critical for the intrinsic formation of anisotropy in shale rock.The proportions of quartz-carbonates-clays could result in very different rock mechanical properties.Clay-rich shale mostly contains highly developed weak planes.Shale samples containing high levels of quartz and low levels of clay have a relatively high Young's modulus and low Poisson's ratio, contributing to a high brittleness and helping hydraulic fracture propagation.
(2) The anisotropy of the wave velocity is thought to result in some geophysical interpretative problems for shale.The effects of mineral grain inclination and micro-cracks are the main factors attributed to an anisotropy.The wave velocities parallel to the bedding are higher than those normal to the bedding, as the ultrasonic wave energy is attenuated by the bedding's scattering.The attenuation anisotropy is affected by the stress, the bedding inclination, and the aspect ratios of micro-cracks.
(3) The mechanical properties of shale specimens under compression show distinct anisotropy.Most determinations of the compression strengths of shale specimens display a "U" type or a similar "shoulder" type variation trend with increasing bedding angle.Particularly, the maximum peak strength is observed at a θ of 0° or 90°, and the minimum values are observed at 45° < θ < 75°.The elastic modulus shows an increasing variation with the bedding inclination increases overall; this is in accord with the variation of ultrasonic velocities.The confining stress reduces the anisotropy of rock material, as the confining pressure is hydrostatic pressure, and when it is applied to the specimen, the beddings are compacted, resulting in the reduction of the primary defects effect.The fracture propagation during the failure process is generally affected by the bedding inclination and confining pressure.For bedding inclination of 0° < θ < 30°, tensile fracture along the loading direction occurs under low confining pressure, and shear fracture across the bedding plane occurs under high confining pressure; for a bedding inclination of 30° < θ < 75°, the fracture generally occurs by the sliding of the bedding, although shear fracture in the matrix can occur, and may cross the bedding; for 75° < θ < 90°, the failure occurs by bedding splitting when under low confining stress, and by bedding shearing under high confining stress.
(4) The tensile strength has a distinct effect on fracture initiation and propagation.The BTS variation in shale specimens from different sampling sites shows a slight difference; in general, their BTS values do not have much change from bedding inclination of 0° to 30°, and then decrease to a low level around 45°-60°, and finally keep within a stable range of around 75°-90°.The center of the disc specimen is in a compressive state when under low bedding inclination and experiencing shearing and tension when under a high bedding inclination.The layer activation occurs under a high bedding inclination, and the central fracture occurs under a low bedding inclination.
(5) The value KIc, ST is much lower than the KIc, A and KIc, D; additionally, KIc, A is slightly greater than KIc, D. The bedding inclination has a distinct effect on the fracture toughness, the values for which all decrease as the bedding inclination increases.The high approach angle and low bedding strength contribute to the fracture propagation along the bedding planes.(6) In shale hydraulic fracturing, the brittleness index (BI) is usually employed to quantify the brittleness of the shale rock formation, for the purpose of evaluating the possibilities of fracture propagation (fracability) and fracture network generation.The BI definitions considering mechanical response are more significant.The horizontal well for hydraulic fracturing usually crosses the bedding plane in different orientations; thus, shale anisotropy should not be ignored in brittleness evaluation.The BIs, which are determined by parameters of Young's modulus and Poisson's ratio, can vary significantly when drilling in shale formations with strong anisotropy.(7) The bedding planes (BP) have an important effect on the fracturing characteristics.The breakdown pressure shows a similar variation with the BTS of shale specimens, in which the maximum and minimum values occur at the bedding inclinations of 0° and 90°, respectively.A smaller bedding inclination contributes to more stress components on the bedding planes, enhancing compression on the beddings and resulting in higher breakdown pressure.The bedding shear stress has an important effect on whether the HFs cross the rock formation or extend along the BP, as well as the fact that the strength of the BP may not be the same at different distances, leading to differences in the bedding opening, in which case some HFs may cross the BP and some terminate at the BP.Moreover, the activation of the BP is the optimum condition for forming fracture networks, as opposed to a simple fracture.

Figure 2 .
Figure 2. Mineral components of shale at different formation depths for the (a) shale of the Wufeng-Longmaxi Formation in well XK-3 [76]; (b) the mineralogy of target depths shows the volumetric fractions of constituents [50].

Figure 3 .
Figure 3. Preparation of shale specimens in different bedding inclinations.

Figure 6 .
Figure 6.A schematic diagram of longitudinal wave propagation in transversely isotropic rock.(Adapted from Yin and Yang [93]).

Figure 7 .
Figure 7. Three types of strength variation relative to the bedding inclination of the anisotropic rock.

Figure 10 .
Figure 10.Strength anisotropy degree variation with respect to the bedding inclination of shale specimens.

Figure 14 .
Figure 14.Variation of average failure load, relative to the failure load for loading perpendicular to weak planes, for shale specimens[25,33,81,95,135,136,137,138].

Figure 15 .
Figure 15.The variation of the stress state at the center with respect to bedding inclination [104].

Figure 20 .
Figure 20.Variation of KIc with respect to the bedding inclination of the shale specimen [37,170].

Figure 26 .
Figure 26.Some basic modes of fracture propagation in a vertical plane, as determined by experiments [214].

Table 2 .
Summary of some common test methods and configurations for KIc determination of a rock specimen.

Table 3 .
Summary of the shale rock KIc values obtained from different studies.

Table 4 .
Summary of some universal BI definitions.