Experimental Study on the Difference of Shale Mechanical Properties

)is paper studies the anisotropic characteristics of shale and the difference in mechanical performance between deep shale and outcrop shale. )e outcrop shale was collected from the Shuanghe section in Changning County, southern Sichuan, and the deep shale was collected from the Wells Yi201 and Lu202. Study their basic mechanical parameters, failure modes, and wave velocity responses through laboratory tests. Research shows that with the increase of bedding angle, the deformation mode has the trend from elastic deformation to plastic deformation in high-stress state.When the bedding angles are 0°, 30°, and 45°, the weak bedding surface plays a leading role in the formation of the failure surface trend. As the bedding angle increases to 60° and 90°, its influence is weakened.)e tensile strength, elastic modulus, and wave velocity decrease with the increase of bedding angle.)e compressive strength and Poisson’s ratio have the law of U-type change, there are higher values at 0° and 90°, and the lowest values are at 30°. )e brittleness index first increases and then decreases with the increase of the bedding angle. )e tensile strength and Poisson’s ratio of outcrop shale and deep shale are close, but the compressive strength of deep shale is only 1/3 of outcrop shale, the elastic modulus is only 3/4 of outcrop shale, and the failure of deep shale is accompanied by instability failure.


Introduction
Shale gas will play an important role in the global energy structure in the future [1,2]. Shale is not only the source rock but also the reservoir of shale gas. Its mechanical characteristics in different directions (that is anisotropy) are one of the important factors to be considered in shale gas drilling stability, hydraulic fracturing design, and field construction [3][4][5][6].
At present, many scholars have made great progress in the study of shale anisotropy. Jaeger et al. [7] has carried on the pioneering research to the anisotropic destruction criterion of rock materials. Lekhnitskii et al. [8] general equation of anisotropic elastic theory is derived, which provides a theoretical basis for the study of anisotropic problems. On these bases, Nova [9] and Cazacu et al. [10] and Lee and Pietruszczak [11] put forward different anisotropic strength criteria. Vernik and Nur [12] studied the wave velocity anisotropy of black organic shale and concluded that the anisotropy of black shale is mainly caused by its own microstructure. Niandou et al. [13] have carried out triaxial tests and addition and unloading tests on shale. e test results show that shale has obvious anisotropic characteristics and has two failure modes: shear failure and tension failure. Cho et al. [14] studied the anisotropy of elastic parameters and strength of gneiss, shale, and schist at different angles through uniaxial compression and Brazilian splitting tests. Sone and Mark [15], through the comparison of shale anisotropy in different blocks, considered that shale anisotropy increases with the increase of clay and organic matter content. Chalmers et al. [16] investigated the nanopore system of gas shale by using scanning electron microscope and got the rule that the pores parallel to the bedding direction are larger. Li et al. [17] defined the stress threshold to study the crack evolution process and gave the conclusion that crack damage stress, crack initiation strain, and crack damage strain are strongly dependent on the bedding layer inclined angle. Wang et al. [18] concluded that the maximum shear strength of shale is 60°based on the experimental data. rough the Brazilian split test and threepoint bending test, Jin et al. [19] discussed the anisotropy of shale fracture behavior. Tang and Wu [20] gave the anisotropic crack evolution characteristics of shale and mudstone in laboratory experiments. Ma et al. [21] proposed a novel theoretical method to underpin the failure mechanisms of anisotropic rocks under Brazilian test conditions. Gui et al. [22,23] reported on anisotropic damage of shale stress coupled with hydration and proposed the conclusion that damage parameters parallel to bedding are greater than those perpendicular to bedding. Wan et al. [24] analyzed the anisotropy of the diffusivity of gas shale, and the diffusivity is much higher at the parallel to the bedding direction than that at perpendicular direction. e shale has obvious bedding characteristics under the influence of deposition mode. e existing studies show that the shale mechanical parameters of different blocks are diversified, coupled with the heterogeneity of shale and the difficulty of sample preparation far, and there is no conclusion on the complete system of shale anisotropy. Based on the previous research experience, this paper makes an indepth study of outcrop shale of Longmaxi Formation in Sichuan Basin on the difference of multidirectional mechanical properties of shale. e stress-strain curve, compressive strength, tensile strength, elastic modulus, Poisson's ratio, brittleness index, and failure mode of shale during shale failure are analyzed. e anisotropy of mechanical parameters is compared, and the elastic constants of Longmaxi shale based on transverse isotropic model are obtained.
In addition, the basic research of shale mechanical parameters is used for underground shale gas exploration. Most of the existing research data are the results of outcrop shale measurements. e mechanical properties of outcrop/ deep shale are not completely the same due to the influence of stratum structure, light, and temperature and cannot be applied arbitrarily. erefore, this paper selects the deep shale of Longmaxi Formation in Sichuan Basin, China, to carry out general mechanical experiments, and makes a comparative analysis and summary of the similarities and differences between outcrop and deep shale mechanical parameters, in order to better match and apply the existing research results of outcrop shale anisotropy to deep shale gas reservoir exploitation in the future and provide basic data support for shale gas exploration and development.

Experimental Sample.
e southern Sichuan Basin is the largest shale gas production base in China. In recent years, PetroChina has vigorously promoted the exploration and development of shale gas resources in southern Sichuan, in which the Longmaxi Formation shale is the main exploration block [25]. erefore, this paper selects the Longmaxi Formation shale as the research object. e outcrop shale sample is located in Shuanghe section of Changning County, Yibin City, Sichuan Province.
is section is exposed to relatively complete and easy to sample around the quarry. e deep shale samples were collected from Well Yi201 in Changning shale gas exploitation block and from Well Lu202 in Luzhou block, respectively, two groups of samples are located in Long 1 1 1 strata of Longmaxi Formation (see document [26] for stratigraphic division), and the specific location is shown in Figure 1. Large shale samples without weathering and with obvious bedding characteristics were selected from the surface. Establish a Cartesian coordinate system based on bedding plane, as shown in Figure 2(a), and prepare standard shale sample with diameter of 25 mm and height diameter ratio of 2 : 1 according to GB/T 50266-2013 standard practice [27]. According to different azimuth angle α and bedding angle β (the angle between axial loading direction and sedimentary bedding plane of sample is shown in Figures 2(c) and 2(d)), shale samples in 12 directions were prepared as shown in Figure 2(b). To ensure that the samples are located in the same coordinate system, all standard shale samples are made on the same large shale sample. A total of 72 outcrop shale samples and 6 deep shale samples are obtained in this test, as shown in Figure 3.

Experimental Equipment and Scheme.
Shale uniaxial compression test and Brazilian splitting test instruments use SANS CMT-5504 electronic universal testing machine with maximum axial load up to 1000 kN and strain data recorded by the DH5923N dynamic signal test and analysis system with accuracy less than 0.5%. Displacement loading control is used in loading mode, and the rate is set 0.02 mm/s. e instrument for measuring wave velocity is based on DS5-8B series full information acoustic emission signal analyzer, whose probe frequency is 50∼400 KHz and the center frequency is 150 KHz. e specific research programme is shown in Table 1.

Analysis of the Variation of Stress-Strain Curve with
Bedding Angle. e stress-strain curves of outcrop shale with square angle α � 0°, 30°, 45°, 60°, and 90°(coring directions 1, 9, 10, 11, and 12) are shown in Figure 4. e stress-strain curves of outcrop shale of Longmaxi Formation show no obvious distinction in each stage, which indicates that shale is typical brittle rock. e initial positions of stressstrain curves of five different bedding angles are all similar to straight lines, and the difference between the initial compaction stage and the linear elastic stage is not obvious, which indicates that there is no obvious compaction stage; the elastic stage of the curve is long, and the distinction between the elastic limit and the plastic limit is not obvious, especially the curves of 0°, 30°, and 45°; after the curve reaches its peak value, the curve drops to its lowest value immediately, and it has typical brittle failure characteristics. e stress-strain curves of 5 different bedding angles coincide at the initial compaction stage, and the strain increases gradually with the increase of axial stress. e turning point of 5 curve strike occurs near the axial strain ε � 4.22 × 10 −3 . e five curves show different trends, from top to bottom, which are 0°, 30°, 45°, 60°, and 90°bedding angle. e larger bedding angle, the lower slope, indicating that the elastic modulus decreases gradually with the increase of bedding angle. 2 Advances in Civil Engineering  e change of bedding angle will affect the deformation mode before the sample is destroyed. e elastic stage of the stress-strain curve is longer, especially the curve of 0°, 30°, and 45°approaches a straight line, indicating that most of the deformation is elastic after compaction and before the failure of the specimens. When the bedding angle increases to 60°a nd 90°, the slope of the curve decreases in the high-stress section before the sample is destroyed, which indicates that the curve has plastic deformation after reaching the elastic limit. With the increase of bedding angle, the loading deformation of shale under high stress tends to change from elastic deformation to plastic deformation. e larger the bedding angle, the longer the plastic deformation stage.

Analysis of the Variation of Shale Failure Type with
Bedding Angle. Because of the tight nature of shale, bedding as a weak surface is the main factor affecting the type of shale failure. e experimental results show that there are two main failure modes of shale under uniaxial load: tensile splitting failure and shear failure. e failure modes with different bedding angles are shown in Figure 5.
(1) When β � 0°, the failure type is vertical splitting tensile failure along the bedding surface. After failure, multiple tensile failure planes parallel to the bedding plane are formed through the two ends of the specimen. e specimens are separated into a plurality of vertical long strip flake rock blocks. (2) When β � 30°, the failure type is shear failure along the bedding plane. After failure, a single shear failure surface along the weak plane of bedding is formed, accompanied with vertical microcracks. ere is obvious shear slip along the bedding plane direction.
(3) When β � 45°, the failure type is a mixture of shear failure along bedding plane and Y-type shear failure. After failure, a shear failure surface along the weak plane of bedding is formed. Under the influence of the end effect, a Y-type shear failure surface is formed at the end of the specimen. (4) When β � 60°, the failure type is a mixture of vertical splitting failure and Y-type shear failure. After fracture, inclined shear failure surface and vertical splitting failure surface are formed passing through many bedding planes. e failure is accompanied by loud noise, and the sample is easy to burst into small rock blocks. (5) When β � 90°, the failure type is a mixture of splitting tension failure of the cross-cut bedding plane and shear failure along the direction of the bedding. After failure, the vertical splitting failure surface and the transverse failure surface along the bedding plane are formed. e two failure surfaces are intersecting.
e experimental results show that with the increase of bedding angle, the failure mode of shale firstly changes from vertical splitting tensile failure to shear slip failure along the weak bedding plane and finally evolves into splitting tension failure of the cross-cut bedding plane. When the bedding angles are 0°, 30°, and 45°, the weak bedding surface plays a leading role in the formation of the failure surface trend. As the bedding angle increases to 60°and 90°, its influence is weakened.

Analysis of Variation of Shale Mechanical Parameters with
Bedding Angle. Figure 6 shows the compressive strength of shale in different directions. Dotted lines represent the sampling angle, and underlined numbers represent the average compressive strength of the group. It can be seen from the figure that the uniaxial compressive strength of outcrop shale of Longmaxi Formation is highly discrete. e uniaxial compressive strength ranges from 126 to 189 MPa, which belongs to high-strength rock.  Advances in Civil Engineering Figure 7 shows the variation of compressive strength with bedding. e curves are plotted with the compressive strength data of bedding angle β � 0°, 30°, 45°, 60°, and 90°at azimuth angle α � 0°. From the curve, it can be seen that the compressive strength of shale has a U-type change law with the increase of bedding angle. When the bedding angle is 0°a nd 90°, the compressive strength is higher, and the maximum value of 90°is 189 MPa. When the bedding angle is 30°, the minimum compressive strength is 126 MPa. e difference between the maximum and minimum compressive strength is 63 MPa. When the bedding angle is 45°∼60°, the slope of the curve decreases, which indicates that in this interval the effect of bedding angle on compressive strength is weakened. Figure 8 shows the variation of elastic modulus and Poisson's ratio with bedding. e curves are plotted with the elastic modulus and Poisson's ratio data of bedding angle β � 0°, 30°, 45°, 60°, and 90°at azimuth angle α � 0°. e elastic modulus decreases with the increase of bedding angle, with a maximum value of 7.51 GPa at 0°and a minimum value of 5.22 GPa at 90°. e changing trend of elastic modulus is slow at both ends and fast in the middle. e change trend of Poisson's ratio is similar to that of compressive strength, showing a U-type change rule of high, middle, and low between two ends. e maximum Poisson's ratio is 0.4 at 0°, and the minimum value is 0.19 at 30°. Table 2 records the tensile strength of outcrop shale with different bedding angles. e curve of tensile strength changing with bedding angle is drawn as shown in Figure 9.
e tensile strength has strong discreteness, and the value of each group of three samples has a difference of about 1 MPa even under the same angle. In particular, the values of tensile  Advances in Civil Engineering 5 strength are more discrete at 0°and 15°. e curve shows that the tensile strength is significantly affected by bedding angle. When the bedding angle is 0°∼15°, the dispersion of tensile strength value is higher but has no obvious change rule with angle. When the bedding angle increases to 30°, the tensile strength increases with the increase of bedding angle. e tensile strength value has a minimum average value of 3.82 MPa at 0°∼15°and has a maximum average value of 7.59 MPa at 90°.

Analysis of Variation of Wave Velocity with Bedding
Angle. Due to the energy attenuation of ultrasonic wave propagating through the weak bedding plane in shale, the velocity of ultrasonic wave propagating in shale shows anisotropic characteristics. e wave velocity of ultrasonic wave passing through bedding in different directions is measured as shown in Figure 10. e test results and the variation of wave velocity with bedding angle are shown in Table 3 and Figure 11. Figure 11 shows the variation of ultrasonic wave velocity with bedding angle. e wave velocity is obviously affected by bedding angle and decreases with the increase of bedding angle. e transverse wave velocity has a maximum value of 2586 m/s at 0°and a minimum value of 2280 m/s at 90°. ere is a difference of 306 m/s between the maximum and minimum. e vertical wave velocity has a maximum value of 4820 m/s at 0°and a minimum value of 4065 m/s at 90°. ere is a difference of 755 m/s between the maximum and minimum.

Difference Analysis of Mechanical Performance between
Outcrop Shale and Deep Shale. In order to directly compare the differences between deep shale and outcrop shale, this paper measures the mineral composition and microstructure of deep shale and outcrop shale. Figure 12 shows the mineral content analysis results of deep shale and outcrop shale. e deep shale is rich in brittle minerals, including quartz and feldspar, accounting for           Figure 13 shows the microstructure of deep shale and outcrop shale. ere is no significant difference in structure between the two at 800 magnification. No large pores or natural cracks are seen in the picture, suggesting that both shales have dense properties. At 3000 magnification, quartz and pyrite grains can be clearly observed embedded in the shale matrix. Deep shale has more lamellar foil structure than outcrop shale, and the area of asphaltene is larger than that of outcrop shale. e mineral components of outcrop shale are arranged more closely.
In this paper, the basic mechanical parameters of deep shale are measured. Since the preparation of deep samples is difficult, only samples with bedding angle 90°are obtained in this experiment, and six samples as shown in Figure 3 e results of uniaxial compression and Brazilian splitting tests on deep shale samples are shown in Table 4. e data in Table 4 show that there are differences between the mechanical characteristics of deep shale and outcrop shale. Specifically, the compressive strength of deep shale is 57 MPa, which is about 1/3 of outcrop shale, and the elastic modulus is 3.56 GPa∼ 4.80 GP, which is about 3/4 of outcrop shale. In addition, the tensile strength and Poisson's's ratio of deep shale and outcrop shale are close. e tensile strength ranges from 4.24 to 7.88 MPa. Poisson's ratio ranges from 0.28 to 0.36. Figure 14 shows the stress-strain curve and postfailure morphology of deep shale in Well Lu202. When strain is less than 0.56 × 10 −2 , the curvature of the curve increases gradually, which indicates that there is an obvious compaction stage in the process. When the strain is greater than 1.64 × 10 −2 , the curve enters the plastic failure zone, and the curve fluctuates violently and forms two stress drops. e second drop value is relatively large, and it decreases from the peak value of 56.6 MPa to 53 MPa. At this time, the specimen does not lose the supporting force. With the increase of vertical load, when the stress gradually increases to 56 MPa, the specimen loses its supporting force instantly after failure. e whole process is characterized by plastic failure. Figure 14(b) shows the specimen morphology after failure. e specimen has splitting failure along the direction of vertical bedding plane. e long strip broken rock block is scattered on the right side of the shale, and the left side of the failure surface is stripped of the main body to form a transverse instability failure surface in the middle part. It is considered that these two failure modes are the main reasons for the sharp fluctuation of stress-strain curve in the plastic failure zone. With the axial load gradually increases to the ultimate load-bearing capacity, the vertical splitting failure occurs in the local area around the cylindrical specimen. At this time, it is accompanied by a stress drop, but the whole specimen is not completely destroyed and still has a supporting force. With the increase of axial stress, the specimen reaches the ultimate bearing capacity and loses its supporting force. is process consists of two parts, namely, the splitting failure of the undamaged part and the instability failure of the slender shale formed by the splitting failure in the previous process.
In conclusion, there are similarities and differences in mechanical performance between deep shale and outcrop shale. e compression stage of deep shale is obvious while that of outcrop shale is not obvious. However, both of them have a long elastic stage and tend to change from elastic deformation to plastic deformation in high-stress state. Compared with outcrop shale, deep shale is accompanied by instability failure when the stress reaches the maximum bearing capacity. However, the angle between loading direction and bedding is the main factor affecting deep/outcrop shale failure mode. In terms of mechanical parameters, the compressive strength and elastic modulus of deep shale are smaller than that of outcrop shale. However, the tensile strength and Poisson's ratio in deep/outcrop shale are close to each other.
Because they have similar mineral composition, it is thought that the difference in microstructure may lead to the different mechanical performance of deep/outcrop shale. e experimental results indicate that it is feasible to study deep shale with outcrop shale as reference, but the differences of mechanical performance between outcrop shale and deep shale should be considered when applying outcrop shale research results to deep shale.

Anisotropic Analysis of Mechanical
Parameters of Shale e property that the mechanical parameters of shale change in one direction with the change in this direction is called anisotropy.
e third section shows that the mechanical parameters of shale change with the bedding angle. e following is an analysis of the anisotropy of the mechanical parameters including compressive strength, tensile strength, elastic modulus, Poisson's ratio, brittleness, and velocity of ultrasonic wave propagation with the change of bedding angle.

Anisotropic Analysis of Compressive/Tensile Strength.
Referring to the degree of strength anisotropy defined by the Saroglou and Tsiambaos [28], this paper introduces a new degree of anisotropy "A" which can measure the degree of variation of mechanical parameters with bedding angle. We define the degree of anisotropy of compressive strength as where A c is the degree of anisotropy of compressive strength; σ cmax is the maximum value of compressive strength; and σ cmin is the minimum value of compressive strength.

Advances in Civil Engineering
According to formula (1), the degree of anisotropy of compressive strength A C � 1.5. rough analysis, it is considered that the angle between the bedding plane and the principal stress leads to the anisotropic characteristics of compressive strength. When the bedding angle is 30°, the shear slip failure occurs along the direction of the weak bedding plane. e compressive strength largely depends on the shear capacity of the weak bedding plane. e low shear strength of the bedding plane directly leads to the reduction of the compressive strength. When the bedding angle is 0°and 90°, the compressive strength mainly depends on the shale matrix itself, so the compressive strength is higher than that at the bedding angle of 30°.
Define the degree of anisotropy of tensile strength as where A T is the degree of anisotropy of tensile strength; T max is the maximum value of tensile strength; and T min is the minimum value of tensile strength. According to formula (2), the degree of anisotropy of shale tensile strength A T � 2.0. It is considered that the   difference of failure modes under different bedding angles leads to the anisotropy of tensile strength. When the bedding angle is small, the failure mode is splitting failure along the direction of bedding plane. In this case, the tensile strength is mainly dependent on the tensile strength of weak bedding plane, so the tensile strength is lower. When the bedding angle is 90°, the failure mode of shale is splitting failure through the bedding plane. Its tensile strength is affected by the shale matrix; thus, the tensile strength increases by two times than 0°.

Analysis of Anisotropy of Elastic Modulus, Poisson's Ratio, and Brittleness Index.
Similarly, the anisotropy of elastic modulus and Poisson's ratio are defined as where A E is the degree of anisotropy of elastic modulus; E max and E min are the maximum and minimum values of elastic modulus; A υ is the degree of anisotropy of Poisson's ratio; and υ max and υ min are the maximum and minimum values of Poisson's ratio. According to formula (3), the degree of anisotropy of elastic modulus A E � 1.4 and Poisson's ratio A υ � 2.1. e elastic modulus and Poisson's ratio reflect the brittleness of the specimen. Using the normalized average method of elastic modulus and Poisson's ratio proposed by Rickman and Mullen [29], the brittleness index of each bedding angle direction is obtained: where BI is the brittleness coefficient; E st is the static elastic modulus in 10 GPa units; υ st is the static Poisson's ratio; E no is the normalized elastic modulus; and υ no is the normalized Poisson's ratio. Figure 15 shows the variation of brittleness index with bedding angle. e curves are plotted with the brittleness index data of bedding angle β � 0°, 30°, 45°, 60°, and 90°at azimuth angle α � 0°.With the change of bedding angle, the brittleness index increases first and then decreases. ere is a maximum brittleness index of 0.33 in the direction of 30°, and a small brittleness index of 0.05 and 0.03 in the direction of 0°and 90°.
Define the degree of anisotropy of brittleness index as where A BI is the degree of anisotropy of brittleness index and BI max and BI min are the maximum and minimum values of brittleness index. According to formula (5), the anisotropy of brittleness index A BI � 11, which indicates that bedding angle has great influence on shale brittleness index.

Analysis of Anisotropy of Ultrasonic Wave Velocity.
Define the degree of anisotropy of wave velocity as where A v is the degree of anisotropy of wave velocity; V 0°i s the velocity of ultrasonic wave in the direction of 0°; V 90°i s the velocity of ultrasonic wave in the direction of 90°. According to formula (6), the anisotropy of transverse wave velocity A Vb � 1.13, and the anisotropy of vertical wave velocity A Vp � 1.19. After analysis, it is considered that the anisotropy of wave velocity is caused by the blocking effect of weak bedding plane. When the bedding angle is 0°, the direction of ultrasonic propagation is consistent with the direction of bedding and does not need to cross the bedding plane. When the bedding angle is 90°, the ultrasonic propagation direction is vertical to the bedding plane. In this case, the mineral components, clay, and complex spatial morphology of the bedding plane impede the propagation of waves, so the wave velocity is slower.
When the bedding angle increases from 0°to 90°the propagation direction of ultrasonic wave gradually changes from parallel bedding direction to vertical bedding direction. e number of layers of bedding plane needed to be penetrated increases, and the energy attenuation in the propagation process increases. So, the propagation velocity of wave tends to slow down gradually.
In summary, due to the differences in composition and spatial structure between the bedding plane and shale matrix, shale has different mechanical characterization along different bedding directions. In this paper, the degree of

Shale Constitutive Model and Parameter Calculation
Taking the bedding plane as the basic plane, the spatial coordinate system is established, as shown in Figure 2. In order to obtain the influence of the azimuth on the behavior of shale mechanics, this paper conducted two groups of experiments for comparison. e first group of experiments measured the mechanical parameters in different azimuth directions of α � 0°, 30°, 45°, 60°, and 90°in the bedding plane (directions 1, 2, 3, 4, and 5). e second group of experiments measured the direction of the same bedding angle in two vertical planes (direction 1 and 5, 12 and 6, 11 and 7, and 10 and 8). Figure 16 shows the variation of compressive strength with azimuth angle. e compressive strength in the direction of bedding plane is discrete, even if the three compressive strengths in the same direction are not exactly the same. Removing the deviation of compressive strength values of individual samples caused by shale discrete characteristics, the values of compressive strength fluctuate slightly with azimuth change from 163 to 183 MPa. e fluctuation range is within 12.27%, and the change trend along the azimuth direction is not obvious. Figure 17 shows the variation of compressive strength with bedding in two vertical planes. In the direction of two identical bedding angles perpendicular to each other, the compressive strength of the shale is always similar. e compressive strength in two vertical planes has the same U-type change rule with bedding angle. e higher coincidence rate of the two U-type change rule curves indicates that the change of azimuth angle does not affect the anisotropy.
Like shale, this material has the same elastic parameters in all directions (called transverse) of the parallel bedding plane, unlike the elastic parameters in the vertical direction (called longitudinal), which is called transverse isotropic in engineering. e number of independent elastic constants for transverse isotropies decreases from 21 to 5 compared to extreme anisotropic elastomers [30].
As shown in Figure 18, a transverse isotropic coordinate system consists of a global coordinate system (x, y, z) and a local coordinate system (x', y', z'). In the local coordinate system (x', y', z'), the z' axis and the normal direction of the bedding plane are the same, while the x' axis and y' axis are located in the bedding plane, and the stressstrain relationship can be expressed by the following formula [31]: where ε ′ is the strain tensor; D ′ − 1 is the flexibility matrix; and σ ′ is the stress tensor: where E 1 is the elastic modulus in the bedding plane; E 2 is the elastic modulus perpendicular to bedding direction; υ 1 is Poisson's ratio in the bedding plane; υ 2 is Poisson's ratio perpendicular to bedding direction; G 1 is the shear modulus in the bedding plane; and G 2 is the shear modulus perpendicular to bedding direction. Because G 1 � E 1 /2(1+υ 1 ), G 1 is not an independent elastic constant. To describe the transversely isotropic elastic body of shale, there are only 5 independent elastic constants E 1 , E 2 , υ 1 , υ 2 , and G 2 .
e first four parameters can be directly measured by the bedding angle β � 0°and 90°samples. e shear modulus perpendicular to the bedding direction G 2 can be obtained by fitting the elastic modulus in the direction of multiple bedding angles.
By means of coordinate system transformation, Niandou et al. [13] obtained the elastic modulus E β of transversely isotropic body under any bedding angle: e following formula can be obtained from formula (9): (10) Formula (10), it can be rewritten into the following form: As shown in Figure 19, the reciprocal elastic modulus of Longmaxi Formation outcrop shale under different bedding angles is fitted by the least square method. e following where E β is the elastic modulus in the direction of β angle with bedding plane. By comparing the fitting curve formula (12) with formula (11), the elastic constant G 2 can be retrieved. 5 independent elastic constants of outcrop shale of Longmaxi Formation can be obtained as shown in Table 5.

Conclusion
(1) e stress-strain curve of outcrop shale has no obvious compaction stage. e elastic stage is longer, and the distinction between elastic stage and plastic stage is not obvious. With the increase of bedding angle, the deformation mode trends to from elastic deformation to plastic deformation in high-stress state. e larger the bedding angle, the longer the plastic deformation stage. (2) When the bedding angle is 0°, 30°, and 45°, the direction of failure plane is basically the same as that of bedding plane. Tensile failure and shear failure occur along the bedding direction, and the bedding plane plays a leading role in the formation of failure surface. When the bedding angle is 60°and 90°, the influence of the bedding plane on the trend of failure surface is weakened. Under the action of the maximum principal stress, the inclined shear failure surface along the bedding plane and the vertical tensile splitting failure surface along the direction of maximum principal stress are produced.
(3) e anisotropy of mechanical parameters of shale along different directions of bedding angle is obvious. Tensile strength, elastic modulus, and ultrasonic wave velocity decrease with the increase of bedding angle. Compressive strength and Poisson's ratio show a U-type change pattern, which are higher at 0°/ 90°and lower at 30°. e angle of anisotropy of mechanical parameters affected by bedding from high to low is brittleness index, Poisson's ratio, tensile strength, compressive strength, elastic modulus, vertical wave velocity, and transverse wave velocity. (4)

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.