Tensile and Shear Mechanical Characteristics of Longmaxi Shale Laminae Dependent on the Mineral Composition and Morphology

: Laminae are well developed in shale and generally inﬂuence fracture propagation during hydraulic fracturing. Hence, comprehensively understanding the tension and shear behaviors of shale laminae is crucial. There have been limited systematic studies thus far on the tensile and shear strength as well as fracture morphology of shale laminae. In this study, the Lower Silurian Longmaxi Shale (China) was investigated via Brazilian tensile and angle-varied plate shear tests. Five lamina types were tested, i.e., calcite (Cal), pyrite (Py), organic-enriched (Oc), the interface between Cal and Oc (Cal-Oc), and the interface between Py and Oc (Py-Oc) laminae. Results showed that the tensile strength was in the range 0.43–8.22 MPa, mainly in the order of Cal > Py > Cal-Oc > Py-Oc > Oc. The modes of fracture morphology were highly related to the occurrence, continuity, and mineralogy ﬁllings of laminae. Shear strength parameters were within the range 22.50–29.64 MPa for cohesion and 37.29–43.60 ◦ for internal friction angle. Fracture surface roughness was strongly related to its cohesion. Calcite laminae considerably inﬂuenced the tensile fracturing of shale, suggesting that the geometry and properties of calcite lamina should receive more attention during the design of shale gas exploration.


Introduction
Shale is generally characterized by layered structures and consequent intrinsic anisotropy owing to distinct sequential sedimentation [1]. The frequently developed sandy laminae in shale have a significant influence on the accumulation [2], storage, and permeability [3] of shale gas as well as on the mechanical properties of shales [4,5]. Outcrop observation and core analysis have indicated that black shale deposited in marine environments was abundant with calcite-/pyrite-rich laminae and lenses as well, such as the Longmaxi Shale [6] and the Devonian Horn River Shale [7]. Calcite and pyrite within shale are generally regarded as brittle minerals and classified into the calcite group and quartz group respectively, which can enhance the brittleness of shale and make the shale more likely to generate and maintain fractures when stimulated [8,9].
Quantitative moduli mapping of organic and inorganic components in shale at a nanometer scale demonstrated that the Young's modulus of different components is distinguishable and varies significantly, ranging from 15 MPa for organic matter to nearly 100 MPa for quartz, except for pyrite being the stiffest in the range 250-312 MPa [10,11]. The heterogeneity of mineralogical components is a

Materials and Methods
The laminated marine shale was sampled from the typical outcrop of the Longmaxi Shale Formation, Sichuan Basin, Shizhu County, Chongqing city, southwestern China ( Figure 1). This zone has been targeted as the main shale play for marine shale gas exploitation [25,26]. First, cylindrical cores with a diameter of 50 mm and cuboid cores with a length of 30 mm were processed from larger shale cubes in the direction parallel to the laminae. The disc-shaped specimens with a diameter of 50 mm and a thickness of 25 mm as well as the cubes with a length of 30 mm were then processed. All the end surfaces were ground until smooth. The center regions of discs and the medial cross-sections of cubes were set with one of five types of laminae or interfaces, i.e., calcite lamina (Cal), pyrite lamina (Py), organic-enriched lamina (Oc), the interface of Cal and Oc lamina (Cal-Oc), and the interface of Py and Oc lamina (Py-Oc) (see Figures 2 and 3). A total of 25 disc specimens (Φ 50 mm × 25 mm) for the Brazilian tensile tests and 25 cubic samples (30 mm × 30 mm × 30 mm) for the angle-varied plate shear tests were prepared, with five samples for each lamina or interface, respectively. The microcomputer-controlled electronic universal testing machine MTS-SANS CMT integrated with three sets of shear clamps (see Figure 4) was employed to conduct the Brazilian tensile and angle-varied plate shear tests. Disc samples were loaded at a constant velocity of 0.06 mm/min until failure, while cubic samples were loaded at a constant loading rate of 1 kN/s, with the inclination angle of shear clamps varying from 20°, 30°, and 40°. Images of the specimens with various laminae before and after Brazilian tensile tests were taken using cameras for a fracture morphology analysis. Meanwhile, the fracture surface roughness of broken cubic specimens was measured using the combination of a high-resolution stereo microscope M205A and a Leica Map Start software module.   The microcomputer-controlled electronic universal testing machine MTS-SANS CMT integrated with three sets of shear clamps (see Figure 4) was employed to conduct the Brazilian tensile and angle-varied plate shear tests. Disc samples were loaded at a constant velocity of 0.06 mm/min until failure, while cubic samples were loaded at a constant loading rate of 1 kN/s, with the inclination angle of shear clamps varying from 20°, 30°, and 40°. Images of the specimens with various laminae before and after Brazilian tensile tests were taken using cameras for a fracture morphology analysis. Meanwhile, the fracture surface roughness of broken cubic specimens was measured using the combination of a high-resolution stereo microscope M205A and a Leica Map Start software module.

Tensile Strength
As recommended by the International Society for Rock Mechanics, the tensile strength was calculated as follows [27]: where σt is the tensile strength (MPa), P is the load at failure (N), D is the diameter of the test specimen (mm), and t is the thickness of the specimen measured at the center (mm). Equation (1) was derived for isotropic materials assuming that the specimen failed with a vertical fracture along the loading diameter. Because the typical vertical splitting failure of a disc cannot be fully achieved for transversely isotropic rock material, in this study, the calculated "tensile strength" is set into quotation marks, referring to the maximum tensile stress at the center of the disc specimen once the specimen has lost its total sustaining capacity [28]. Consequently, the formula was adopted only for comparison among specimens with various laminae; further, the results need not characterize the true tensile strength of the samples [29].
The original test numbers BFW01-05, BHW01-05, BYW01-05, BFY01-05, and BHY01-05 named on samples corresponds to the five sample types with laminae of Cal, Py, Oc, Cal-Oc, and Py-Oc, respectively. Therefore, we substituted the original numbers with Cal01B-05B, Py01B-05B, Oc01B-05B, Cal-Oc01B-05B, and Py-Oc01B-05B to make the consequent analysis clearer and easier to compare. Tensile strengths obtained from the Brazilian tensile tests were in the range 0.43-8.22 MPa for shale with different laminae, as listed in Table 1. Since the occurrence of laminae has unavoidable differences with that of design, the tensile strength varied considerably even if the specimen were The microcomputer-controlled electronic universal testing machine MTS-SANS CMT integrated with three sets of shear clamps (see Figure 4) was employed to conduct the Brazilian tensile and angle-varied plate shear tests. Disc samples were loaded at a constant velocity of 0.06 mm/min until failure, while cubic samples were loaded at a constant loading rate of 1 kN/s, with the inclination angle of shear clamps varying from 20 • , 30 • , and 40 • . Images of the specimens with various laminae before and after Brazilian tensile tests were taken using cameras for a fracture morphology analysis. Meanwhile, the fracture surface roughness of broken cubic specimens was measured using the combination of a high-resolution stereo microscope M205A and a Leica Map Start software module.
within those of the same central laminae type. In general, the shale specimens with calcite/pyrite laminae (Cal/Py) had the highest tensile strength, approximately 6~8 MPa, followed by specimens with the interface of calcite/pyrite laminae and organic-enriched laminae (Cal-Oc/Py-Oc) (2~6 MPa), and the lowest tensile strength (0.43-1.66 MPa) was measured in shale with organic-enriched laminae (Oc). Core observations and mechanical tests demonstrate that the pyrite and calcite laminae were often synchronously deposited in natural bedding-parallel fractures. Consequently, the specimens of Cal, Py, Cal-Oc, and Py-Oc mostly had both calcite and pyrite laminae near the center of the discs, and fractures occasionally initiated from or propagated along with their interfaces (Cal-Py) rather than along the designed central laminae, leading to a relatively large fluctuation in the tensile strength of the same lamiae types. From the scatter plot of the tensile strength of the Longmaxi Shale with various laminae (see Figure 5), five prominent levels of tensile strength could be easily classified: I (7-8 MPa), II (6-7 MPa), III (5-6 MPa), IV (2)(3)(4)(5), and V (0-2 MPa), corresponding to specimen with Cal, Py, Cal-Oc, Py-Oc, and Oc, respectively. The results suggest that the strength heterogeneity of shale laminae is apparent, which should be considered for numerical simulation or physical model test design. According to the classification and realignment of the specimens, we recalculated the average tensile strength of each central mineral interlayer type, which were 7.39 MPa, 6.35 MPa, 5.47 MPa, 3.69 MPa, and 0.96 MPa for shale laminae Cal, Py, Cal-Oc, Py-Oc, and Oc, respectively (Table 1).

Tensile Strength
As recommended by the International Society for Rock Mechanics, the tensile strength was calculated as follows [27]: where σ t is the tensile strength (MPa), P is the load at failure (N), D is the diameter of the test specimen (mm), and t is the thickness of the specimen measured at the center (mm). Equation (1) was derived for isotropic materials assuming that the specimen failed with a vertical fracture along the loading diameter. Because the typical vertical splitting failure of a disc cannot be fully achieved for transversely isotropic rock material, in this study, the calculated "tensile strength" is set into quotation marks, referring to the maximum tensile stress at the center of the disc specimen once the specimen has lost its total sustaining capacity [28]. Consequently, the formula was adopted only for comparison among specimens with various laminae; further, the results need not characterize the true tensile strength of the samples [29].
The original test numbers BFW01-05, BHW01-05, BYW01-05, BFY01-05, and BHY01-05 named on samples corresponds to the five sample types with laminae of Cal, Py, Oc, Cal-Oc, and Py-Oc, respectively. Therefore, we substituted the original numbers with Cal01B-05B, Py01B-05B, Oc01B-05B, Cal-Oc01B-05B, and Py-Oc01B-05B to make the consequent analysis clearer and easier to compare. Tensile strengths obtained from the Brazilian tensile tests were in the range 0.43-8.22 MPa for shale with different laminae, as listed in Table 1. Since the occurrence of laminae has unavoidable differences with that of design, the tensile strength varied considerably even if the specimen were within those of the same central laminae type. In general, the shale specimens with calcite/pyrite laminae (Cal/Py) had the highest tensile strength, approximately 6~8 MPa, followed by specimens with the interface of calcite/pyrite laminae and organic-enriched laminae (Cal-Oc/Py-Oc) (2~6 MPa), and the lowest tensile strength (0.43-1.66 MPa) was measured in shale with organic-enriched laminae (Oc). Core observations and mechanical tests demonstrate that the pyrite and calcite laminae were often synchronously deposited in natural bedding-parallel fractures. Consequently, the specimens of Cal, Py, Cal-Oc, and Py-Oc mostly had both calcite and pyrite laminae near the center of the discs, and fractures occasionally initiated from or propagated along with their interfaces (Cal-Py) rather than along the designed central laminae, leading to a relatively large fluctuation in the tensile strength of the same lamina types. From the scatter plot of the tensile strength of the Longmaxi Shale with various laminae (see Figure 5), five prominent levels of tensile strength could be easily classified: I (7-8 MPa), II (6-7 MPa), III (5-6 MPa), IV (2-5 MPa), and V (0-2 MPa), corresponding to specimen with Cal, Py, Cal-Oc, Py-Oc, and Oc, respectively. The results suggest that the strength heterogeneity of shale laminae is apparent, which should be considered for numerical simulation or physical model test design. According to the classification and realignment of the specimens, we recalculated the average tensile strength of each central mineral interlayer type, which were 7.39 MPa, 6.35 MPa, 5.47 MPa, 3.69 MPa, and 0.96 MPa for shale laminae Cal, Py, Cal-Oc, Py-Oc, and Oc, respectively (Table 1).
The shear strength of shale laminae was calculated as follows: where σ and τ are the normal stress and shear stress (MPa), respectively, P is the maximum failure load of the specimen (N), α is the inclination angle ( • ), and A is the shear area of the specimen (mm 2 ). The experimental and averaged normal stress and shear stress at various inclination angle obtained during the angle-varied plate shear test is shown in Figure 6. The scatter points represent the stress values of tested samples and the lines are fitted linearly from corresponding points. The fitting formula is characterized as τ = σtanφ + c, of which tanφ represents the slope of the fitted line and c is the intercept.
To characterize the shear strength characteristics without considering normal stress, the cohesion and internal friction angle of the shale cubes were analyzed according to the linear regression of the normal stress-shear stress curve. As shown in Figure 7 and Table 2, the shear strength parameters of shale laminae have strong heterogeneity due to the difference in mineral components, with a range of 22.50-29.64 MPa for cohesion and 37.29-43.60 • for the internal friction angle. The Cal shale laminae had the highest internal friction angle (43.60 • ) but the lowest cohesion (22.50 MPa). Nevertheless, the specimen with the highest cohesion did not present the smallest internal friction, suggesting that the two parameters were not in a linear negative relation as the mineral components varied.

Tensile Fracture Morphology
Test results suggested that the samples with different lamina had different failure modes. We analyzed the potential corresponding relationship between tensile strength and failure morphology, and the eventual facture morphologies of representative specimens were extracted and realigned according to the tensile strength level, as illustrated in Figure 8.
As mentioned above, there were five levels I-V of tensile strength which could be classified, corresponding to the specimens with Cal, Py, Cal-Oc, Py-Oc, and Oc, respectively. Different lamina had different fracture surfaces. The highest level I failed with tortuous fractures along fluctuations of Cal-Py, while level II tended to fail along the tortuous Cal-Py containing ~1/2 calcite and ~1/2 pyrite. Level III ruptured with straight fractures nearly along Cal-Py or turning from Cal-Py into Py-Oc

Tensile Fracture Morphology
Test results suggested that the samples with different lamina had different failure modes. We analyzed the potential corresponding relationship between tensile strength and failure morphology, Energies 2020, 13, 2977 9 of 18 and the eventual fracture morphologies of representative specimens were extracted and realigned according to the tensile strength level, as illustrated in Figure 8.

Shear Fracture Surface Morphology
The fracture surface of specimens after shear tests was observed and measured using a highresolution stereo microscope M205A, trying to interpret the internal mechanism of the shear strength parameter difference due to various laminae. Even though the fracture surfaces were irregular and rugged, high-resolution Montage multi-focus images can be captured by combining the LAS Montage Multi-focus software module with a microscope, with the depth of focus (namely height) recorded synchronously for each pixel in the images.
The maximum field of the microscope was approximately 17 mm × 13 mm, while the specimen surface was about 30 mm × 30 mm. Thus, we took approximately six Montage images for each specimen surface and these were merged into one final picture. The representative pictures of the fracture surfaces with various laminae are exhibited in Figure 9. Fracture surfaces of the organicenriched laminae seemed to be flattest, while the other four types appeared to fluctuate. As mentioned above, there were five levels I-V of tensile strength which could be classified, corresponding to the specimens with Cal, Py, Cal-Oc, Py-Oc, and Oc, respectively. Different laminae had different fracture surfaces. The highest level I failed with tortuous fractures along fluctuations of Cal-Py, while level II tended to fail along the tortuous Cal-Py containing~1/2 calcite and~1/2 pyrite. Level III ruptured with straight fractures nearly along Cal-Py or turning from Cal-Py into Py-Oc containing low-density pyrite, while level IV corresponded with fractures turning from straight Cal-Py into Oc or just along Py-Oc containing low-density pyrite. The lowest level V was characterized by straight tension fractures and local shear fractures within organic-enriched laminae (Oc).
Thus, the shale specimens have tensile fractures along with tortuous interfaces and generally consume more energy and have higher tensile strength than those along the linear interfaces. Meanwhile, the larger the contact area of fractures with high content calcite or pyrite lamina, the higher the tensile strength. The eventual failure morphology of shale with various tensile strength and laminae are summarized in Table 3. The tensile strength heterogeneity of the Longmaxi Shale with calcite/pyrite laminae (3.69-7.39 MPa) was similar to the variation of Brazilian tensile strength (3.30-7.63 MPa) of the Marcellus Shale with different anisotropy angles [30], suggesting that the strength heterogeneity was as prominent as the anisotropy of shale. Therefore, the influence of both the mineral difference and bedding orientation should be considered for analyzing the tension behavior of shales.

Shear Fracture Surface Morphology
The fracture surface of specimens after shear tests was observed and measured using a high-resolution stereo microscope M205A, trying to interpret the internal mechanism of the shear strength parameter difference due to various laminae. Even though the fracture surfaces were irregular and rugged, high-resolution Montage multi-focus images can be captured by combining the LAS Montage Multi-focus software module with a microscope, with the depth of focus (namely height) recorded synchronously for each pixel in the images.
The maximum field of the microscope was approximately 17 mm × 13 mm, while the specimen surface was about 30 mm × 30 mm. Thus, we took approximately six Montage images for each specimen surface and these were merged into one final picture. The representative pictures of the fracture surfaces with various laminae are exhibited in Figure 9. Fracture surfaces of the organic-enriched laminae seemed to be flattest, while the other four types appeared to fluctuate. Based on the depth of focus recorded for each pixel in the multi-focus images, the statistical analysis of fracture surface roughness was conducted with the cooperation of the Leica Map Start and 3D Viewer software modules. First, the three-dimensional (3D) view of the fracture surface was generated from a two-dimensional (2D) Montage image to reflect the stereo roughness of the shear plane, as shown in Figure 10a,b. The Montages image was then imported into the Leica Map Start to generate the 2D and 3D depth maps with a depth of focus for each pixel in the image characterized by Based on the depth of focus recorded for each pixel in the multi-focus images, the statistical analysis of fracture surface roughness was conducted with the cooperation of the Leica Map Start and 3D Viewer software modules. First, the three-dimensional (3D) view of the fracture surface was generated from a two-dimensional (2D) Montage image to reflect the stereo roughness of the shear plane, as shown in Figure 10a,b. The Montages image was then imported into the Leica Map Start to generate the 2D and 3D depth maps with a depth of focus for each pixel in the image characterized by gradient color (Figure 10c,d). The shades of red represent large values, while those of blue, small values. Based on the depth of focus recorded for each pixel in the multi-focus images, the statistical analysis of fracture surface roughness was conducted with the cooperation of the Leica Map Start and 3D Viewer software modules. First, the three-dimensional (3D) view of the fracture surface was generated from a two-dimensional (2D) Montage image to reflect the stereo roughness of the shear plane, as shown in Figure 10a,b. The Montages image was then imported into the Leica Map Start to generate the 2D and 3D depth maps with a depth of focus for each pixel in the image characterized by gradient color (Figsure 10c,d). The shades of red represent large values, while those of blue, small values. The fracture surfaces of two representative samples for each lamina type are illustrated in Figure 11, among which the height scales of the corresponding gradient color for shale laminae Py-Oc and Oc were the largest and the smallest, respectively. As the root-mean-square height, namely the standard deviation of the height distribution, is among the most commonly used parameters for characterizing surface roughness [31,32], herein we adopted the root-mean-square height to quantify the surface roughness. The larger value indicates a higher surface roughness.
The fracture surfaces of two representative samples for each lamina type are illustrated in Figure 11, among which the height scales of the corresponding gradient color for shale laminae Py-Oc and Oc were the largest and the smallest, respectively. As the root-mean-square height, namely the standard deviation of the height distribution, is among the most commonly used parameters for characterizing surface roughness [31,32], herein we adopted the root-mean-square height to quantify the surface roughness. The larger value indicates a higher surface roughness. The root-mean-square height formula is included in the ISO 25178 standard as follows: where A is the projected area of the fracture surface and z is the height perpendicular to the surface, namely the depth of focus, which is the function of horizontal and vertical coordinates x and y.
Since the final fracture surface of each specimen after the angle-varied plate shear test was merged from several sub-images, the root-mean-square height of each specimen was calculated based on the corresponding roughness of sub-images. Hence, first, the root-mean-square of each sub-image was calculated with the depth of focus of numerous pixels (with a pixel resolution of 2560 × 1920) on the image. The root-mean-square height of each specimen was then obtained by averaging the value of its sub-images. Finally, fracture surface roughness was characterized quantitatively with the average root-mean-square height of five parallel specimens of each type, as summarized in Table 4. The results indicated that the fracture surfaces of shale lamina Py-Oc were the roughest, while the fracture plane of shale lamina Oc was the smoothest. We infer that the largest roughness of shear planes along the interfaces of pyrite and organic-enriched laminae could be attributed to the largest strength difference between pyrite and organic matters. As we compared the fracture surface roughness parameter with shear strength parameters, good consistency was revealed between the variation of cohesion and the fracture roughness due to the mineral composition of the laminae (see Figure 12). In other words, the larger the cohesion, the rougher the fracture surface. This finding could shed light on the optimization of shear fracture roughness according to the features of shear strength parameters. The root-mean-square height formula is included in the ISO 25178 standard as follows: where A is the projected area of the fracture surface and z is the height perpendicular to the surface, namely the depth of focus, which is the function of horizontal and vertical coordinates x and y.
Since the final fracture surface of each specimen after the angle-varied plate shear test was merged from several sub-images, the root-mean-square height of each specimen was calculated based on the corresponding roughness of sub-images. Hence, first, the root-mean-square of each sub-image was calculated with the depth of focus of numerous pixels (with a pixel resolution of 2560 × 1920) on the image. The root-mean-square height of each specimen was then obtained by averaging the value of its sub-images. Finally, fracture surface roughness was characterized quantitatively with the average root-mean-square height of five parallel specimens of each type, as summarized in Table 4. The results indicated that the fracture surfaces of shale lamina Py-Oc were the roughest, while the fracture plane of shale lamina Oc was the smoothest. We infer that the largest roughness of shear planes along the interfaces of pyrite and organic-enriched laminae could be attributed to the largest strength difference between pyrite and organic matters. As we compared the fracture surface roughness parameter with shear strength parameters, good consistency was revealed between the variation of cohesion and the fracture roughness due to the mineral composition of the laminae (see Figure 12). In other words, the larger the cohesion, the rougher the fracture surface. This finding could shed light on the optimization of shear fracture roughness according to the features of shear strength parameters.

Discussion
This study found that pyrite laminae were commonly associated with calcite and that their association had a great influence on the fracture modes of shale. As the fracture position of the sample in shear tests was controllable and the test results revealed the strength of each lamina type, herein we analyze only the tensile test results. Figure 13 shows the association development of pyrite lamina and calcite lamina (upper part) and their broken features after tests (bottom part). As shown in the figure, there were three association types of pyrite and calcite laminae. First, pyrite lamina was embedded in the calcite lamina; second, pyrite developed tightly adjacent to the calcite lamina; and third, the first and second types developed together and coexisted.

Discussion
This study found that pyrite laminae were commonly associated with calcite and that their association had a great influence on the fracture modes of shale. As the fracture position of the sample in shear tests was controllable and the test results revealed the strength of each lamina type, herein we analyze only the tensile test results. Figure 13 shows the association development of pyrite lamina and calcite lamina (upper part) and their broken features after tests (bottom part). As shown in the figure, there were three association types of pyrite and calcite laminae. First, pyrite lamina was embedded in the calcite lamina; second, pyrite developed tightly adjacent to the calcite lamina; and third, the first and second types developed together and coexisted.

Discussion
This study found that pyrite laminae were commonly associated with calcite and that their association had a great influence on the fracture modes of shale. As the fracture position of the sample in shear tests was controllable and the test results revealed the strength of each lamina type, herein we analyze only the tensile test results. Figure 13 shows the association development of pyrite lamina and calcite lamina (upper part) and their broken features after tests (bottom part). As shown in the figure, there were three association types of pyrite and calcite laminae. First, pyrite lamina was embedded in the calcite lamina; second, pyrite developed tightly adjacent to the calcite lamina; and third, the first and second types developed together and coexisted.   Figure 13a shows the typical discs with central calcite laminae (Cal), with the tensile fracture having mainly occurred in the calcite lamina crossing the pyrite lamina that it was embedded in. These specimens were relatively brittle such that four in five specimens broke with half-discs falling down to the bearing plate, likely due to the relatively large energy release within a short time period. Figure 13b demonstrates that the fracture was captured by the interface between pyrite and calcite laminae, although the fracturing began in the pyrite lamina in the upper part. Compared with Cal discs, the broken behaviors of Py discs were relatively gentle, with only one specimen having failed with a half-disc thrown to the ground. While in Figure 13c, the fracture occurred along the tortuous pyrite lamina embedded in calcite lamina, showing various tensile strengths due to the variation of pyrite content and the occurrence of the interface.
Our experimental results are in good agreement with that the fracture patterns impacted by bedding planes and calcite veins under the Brazilian tensile test [17,33,34]. Generally, for the shale loaded parallel to the bedding plane, the diametrical tensile fracture initiates in the middle region but does not follow a single layer or bedding plane and is also not perfectly straight but slightly curved [33]. When the contained calcite-filled discontinuities are along the applied force or crossed with small angle, the discontinuities are easy to activate and tensile failure mainly occurs along them [17]. However, different from the profoundly reduced tensile strength of shale with sealed fractures or calcite veins, nearly half of the sample without the fractures or veins [34,35], the tensile strength of shale containing calcite laminae in our tests are measured to be much higher than the shale solely composed with organic-rich laminae, as we have demonstrated in Table 1 and Figure 5. The probable reason for it is that the thickness of mineral laminae in our tests, approximately 1.5-4.0 mm, is much larger than the veins or sealed fractures, about 0.17 mm by Mighani et al. [34]. Consequently, the thicker calcite laminae have coarser crystal grains and the interface between calcite laminae and organic-rich matrix is a strong bond rather than a plane of weakness indicated by Padin et al. [36]. Overall, the fracture modes and tensile strength differences show that the calcite laminae have a great influence on the tensile fracture of shale. The geometry and the mechanical properties of calcite laminae should be taken into consideration during the design of shale gas exploration. The mechanical behavior of shale laminae in this study is mainly based on laboratory uniaxial tests on natural shale rocks. Integrating by experimental studies and numerical simulations, the future investigations on this topic mainly focus on the mechanical properties of shale laminae under various fracturing conditions, the process detection of crack initiation, propagation, and eventual failure during mechanical tests of shale laminae, and the in-situ mechanical response of shale laminae to fracturing in field-scale. These further research achievements would provide more optimization data of fracturing design for shale gas exploitation.

Conclusions
This paper presents a series of Brazilian tensile tests and angle-varied plate shear tests on the Longmaxi Shale to investigate the mechanical properties of shale laminae. With the mechanical strength and corresponding fracture surface morphology characterization, the main understandings of the tension and shear behavior of the Longmaxi Shale laminae are summarized as follows: