Study on the mechanism of anisotropic wellbore instability in continental shale in Songliao Basin

Wellbore instability is frequently encountered in the process of shale oil drilling. Due to the anisotropic strength characteristics of shale, most conventional models are not suitable for wellbore stability analysis in layered formation. In this paper, taking the continental shale in the north of Songliao Basin as the research object, the anisotropic elasticity and strength parameters of shale are measured experimentally. Based on pore elastic mechanics, an anisotropic wellbore stability model of layered shale is established, and the variation of collapse pressure under different formation conditions is analyzed. The results show that considering the elastic and strength anisotropy, the collapse pressure is the largest, and the influence of strength anisotropy is significantly greater than that of elastic anisotropy. The stability of horizontal wells drilled along the horizontal maximum principal stress is the optimal, followed by vertical wells and small angle directional wells. The actual drilling conditions of G101 well are basically consistent with the prediction results, which confirms the accuracy of the model. According to the traditional wellbore stability model, the equivalent density of collapse pressure in horizontal well section is 1.58 ~ 1.76 g/cm3. It is much lower than the collapse density of 1.86 g/cm3 calculated by the anisotropic model in this study. The findings of this study can help for better understanding of the mechanism of wellbore instability in horizontal wells in shale formations. The model can be used to guide the drilling engineering design of shale oil horizontal wells and reduce the losses caused by the instability of the wellbore.


Introduction
In recent years, with the development of conventional oil and gas resources has gradually failed to meet the demand of current industrial development, the development of unconventional oil and gas resources such as shale oil has gradually become a focus of attention. However, the wellbore instability of shale formation seriously restricts the efficient development of shale oil. After drilling the formation to form a wellbore, the drilling fluid column pressure replaces the support provided by the formation, breaks the original stress balance of the formation, and causes the stress redistribution of the rocks around the wellbore (Asaka and Holt 2020;Ma and Chen 2015;Ma et al. 2016). If the redistributed stress exceeds the maximum load that the shale can support, it will lead to wellbore instability. In addition, drilling fluid intrudes into the formation, causing the increase of formation pore pressure and the decrease of shale strength, which will further aggravate the instability of wellbore (Ding et al. 2020;Ma et al.2020a, b, c, d;Ma et al. 2020a, b, c, d). In the development of shale oil, the unstable layered shale is usually the target layer to drill long horizontal well, and there is a particularity of itself, and the secondary stress distribution in horizontal wells is quite different from that in straight wells. In addition, due to the complex in situ stress, bedding and fracture development of shale reservoir, the risk of wellbore collapse is aggravated (Liu et al. 2016;Zeynali 2012;Ibrahim 2021;Zarei and Nasiri 2021).
Aiming at the wellbore instability of layered shale, researchers have carried out a lot of research and achieved rich results. Chenevert and Gatlin (1965), Aadnoy (1988); Dusseault and Gray (1992), Lee (2004), Shuai et al. (2015) analyzed the effect of bedding on wellbore stability by adopting single weak plane strength theory and anisotropic strength theory. Liang et al. (2014), Kang et al. (2011), Zhi et al. (2016), Nguyen et al. (2020) studied the relationship between wellbore stability and bedding shale based on single weak plane strength theory and linear elastic wellbore stability model. Chen et al. (2015a, b), Ma et al. (2019), Ma et al. (2020a, b, c, d) studied the effect of bedding plane occurrence and bedding plane strength reduction on the stability of vertical wells.
However, at present, for the instability of shale wellbore, the comprehensive consideration of elasticity and strength anisotropy is not insufficient. Hard and brittle shale strata have weak planes, and the research on formation collapse pressure mostly focuses on the influence of weak planes or formation mechanics anisotropy alone, and does not comprehensively study the dual influence of formation mechanics and strength anisotropy. Moreover, the current research on shale wellbore stability in China mainly focuses on the marine shale of Longmaxi formation in Sichuan, and the wellbore instability of continental shale in the north of Songliao basin has not been studied. Therefore, in this paper, the anisotropic wellbore collapse pressure model is established for the study of layered shale oil horizontal wells. The variation of wellbore collapse pressure of continental shale is analyzed, and an example is analyzed to verify the correctness of this method. The regularity of wellbore collapse pressure of continental shale is analyzed, and the correctness of the method is verified by the analysis of a horizontal well G101 of shale oil in the north of Songliao Basin. The research results of this paper reveal the mechanical mechanism of wellbore instability in anisotropic shale, which can be used to guide the drilling design of shale oil horizontal wells and the prevention of wellbore instability.
The research route of this paper is shown in Fig. 1. First, rock mechanics experiments are carried out for continental shale in Songliao basin. Then, the elastic and strength anisotropy parameters of shale in this block are determined, respectively. Next, the anisotropic wellbore stability model of shale is established. Afterwards, the difference of shale wellbore collapse pressure under three anisotropic conditions is analyzed and discussed. Finally, a well in the oil field is used to verify the accuracy of the model.

Experimental core preparation
Through the observation of full-scale cores, it can be found that the shale bedding and fractures in this block are obviously developed, and the weak plane is clearly visible. Since the shale is prone to hydration expansion when encountering water, it cannot be cut and cored used in conventional processing methods, so the full-size core is processed by using the NC wire cutting method. According to the coring method in Fig. 2, the 0, 30°, 45°, 60°, 90° direction cores are processed, respectively, and both ends of the sample are cut flat and polished to make the length diameter ratio of the shale sample about 2.

X-RD diffraction experiment
The mineral composition of continental shale in Songliao basin is analyzed and tested by X-ray diffraction analysis technology. The main minerals of shale in this block are quartz and clay, the content of quartz is 23.4 ~ 29.4%, and the content of clay minerals is 33.3 ~ 45.6%. There is also a part of dolomite, calcite and plagioclase, which belongs to a typical hard brittle mineral. Quartz and feldspar are the main brittle minerals in shale. The brittle minerals of shale in this block generally exceed 50%, which is an important factor affecting the development degree of pores and microfractures in shale matrix. The higher the content of quartz and feldspar, the greater the brittleness of shale and the higher the possibility of fracture development (Liu et al. 2021). The clay minerals in the shale are mainly illite and chlorite, containing a small amount of illite-smectite mixed layer, without pure smectite minerals, and the content of smectite in illite-smectite mixed layer is not high, indicating that it has certain hydration characteristics, but its hydration expansion is not obvious. Illite content is 53 ~ 68.5%, chlorite content is 6.4 ~ 23.5%.
In the actual drilling operation, the filtrate of drilling fluid enters the formation, and the clay minerals are hydrated and dispersed, and the hydration expansion capacity is weak. However, different clay minerals have different hydration expansion characteristics along different directions. Different clay particles have different water absorption and expansion rates, resulting in different expansion pressures, resulting in unbalanced stress in the formation and reduced formation strength, which is easy to cause spalling and collapse of the formation along the bedding and fracture section. The polarity of water in the drilling filtrate is extremely strong, the hydroxyl group of water and the surface of the clay particles and the surface of the silt particles, replacing the original hydrogen bonds between clay particles, silt particles and clay particles and silt particles. It will reduce the binding force between rock particles, disperse rock particles, reduce the friction strength of bedding and fracture surface. As a result, the shale wellbore is unstable, resulting in serious borehole collapse.

Uniaxial and triaxial compression strength experiment
In order to clarify the mechanical anisotropy characteristics of shale in block G of Songliao Basin, the cores of the three wells are collected, and the uniaxial and triaxial experiments are carried out by TAW2000 microcomputer controlled triaxial stress testing machine. The core crushing morphology after loading in parallel bedding direction and vertical bedding direction is shown in Figs. 3 and 4.
The fracture morphology of shale is different after loading in parallel bedding and vertical bedding directions. When it is loaded in parallel to the bedding direction, it will generate multiple openings the bedding seam to form multiple groups of fragments along the bedding plane. However, the vertical direction of bedding loading will produce the shear crack first, and the shear crack and the micro-crack along the bedding are connected to each other to form a large volume fracture. It can be concluded that the rock loaded parallel to the bedding direction produces tensile failure along the structural plane, and the vertical bedding direction needs to overcome the shear strength of the shale body.
As can be seen from Table 1, since the shale in Block G is a bedding developmental formation, its uniaxial compressive strength is related to the coring angle of the core and is a function of the loading direction and the angle between the bedding planes. The strength of shale samples with sampling angle of 60° is the lowest, the strength of shale samples with sampling angle of 45° is slightly higher, while the strength of shale samples with sampling angle of 90° (perpendicular to bedding plane) is the highest, and the strength of shale samples with sampling angle of 0 (parallel to bedding plane) is lower than that of shale samples with sampling angle of 90°. With the increase of sampling angle, the strength is the highest when the loading direction is perpendicular to the shale bedding plane. According to Mohr Coulomb criterion, when the angle is 45° + φ/2 (φ is the internal friction angle), and then the strength increases gradually. These characteristics fully show that the shale in this block shows significant anisotropic mechanical characteristics.
As shown in Fig. 5, as the increase of confining pressure, the compressive strength of shale is enhanced. When the bedding direction and loading direction are at different angles, the failure form of shale also changes:   (1) When the angle between the loading direction and the bedding is 0° ~ 30°, the splitting failure mainly occurs along the bedding plane.
(2) When the angle between the loading direction and bedding is 30° ~ 75°, shear slip mainly occurs along the bedding plane. (3) When the loading direction and bedding are 75° ~ 90°, the shear failure of shale body mainly occurs, and the strength is the largest.

Determination of anisotropic elastic parameters of shale
Considering that the shale in block G of Songliao basin has significant hard brittle characteristics, the stress-strain curve shows significant elasticity and anisotropy. Therefore, shale can be assumed as transversely isotropic porous elastic medium, and the elastic constitutive equation of transversely isotropic porous elastic medium can be expressed as (Cao et al. 2017): where is the strain vector under stratum coordinate; C is the flexibility matrix under stratum coordinate; is the stress vector under stratum coordinate, MPa; E h , E v are the elastic modulus parallel to the bedding direction and perpendicular to the bedding direction, respectively, MPa;v h , v v are the Poisson's ratio to the bedding direction and perpendicular to the bedding direction, respectively; G h , G v are the shear modulus parallel to the bedding direction and perpendicular to the bedding direction, respectively, MPa. Therefore, according to Eqs. (1) -(4), combined with the uniaxial shale mechanics test results of continental shale in Songliao Basin shown in Table 1, the elastic parameters of the block shale are shown in Table 2:

Determination of anisotropic strength parameters of shale
According to the uniaxial and triaxial compression test results of shale, the continental shale in Songliao basin shows significant anisotropic strength characteristics. If the traditional isotropic Mohr Coulomb criterion is used to characterize the strength of this shale, it will not accurately reflect the strength anisotropy of shale. In this paper, the single weak plane strength theory is used to describe the strength anisotropy of shale. The strength failure criterion of rock can be divided into discontinuous model and continuous model. Because the bedding shale has a large number of discontinuous weak planes, the discontinuous model can effectively predict the anisotropic strength of shale. Chenevert et al. (1965) confirmed through experiments that when the angle between the normal of the bedding plane and the maximum principal stress is 20° ~ 30°, its strength is about 40% lower than that of coring along the direction perpendicular to the bedding plane. In this paper, Jaeger's single weak plane strength failure criterion is used to describe the strength anisotropy of layered shale. A group of weak planes AB are developed in shale. It is assumed that the angle between AB plane (referring to its normal direction) and the direction of maximum principal stress is . Jaeger (1960) studied layered shale and established a single weak plane criterion. There are two forms of shear failure of layered stratum: shear failure along bedding plane and shear failure of rock matrix, which is expressed as: Failure model of shale matrix: Failure model of weak plane: where 1 , 3 are the maximum principal stress and the minimum principal stress, respectively, MPa; c 0 , c w are the cohesion of shale matrix and weak plane, respectively, MPa; 0 , w are the internal friction angle of shale matrix and weak plane, respectively, • ; 0 , are the angle between the failure plane of shale body and the maximum principal stress, and between the normal of weak plane and the maximum principal stress, respectively, • ; 0 = 4 + 0 .
(5)  According to the uniaxial and triaxial compression test results of shale, the mechanical parameters of shale are obtained by grid search method (Ma et al. 2020a, b, c, d). The theoretical strength under different confining pressures is calculated and compared with the measured results. The results are shown in Fig. 7.
As a typical structure of shale, the weak surface will affect the overall rock strength. However, rock strength is not only related to structure, but also related to its stress state. Therefore, by setting different confining pressure conditions, the relation between the strength of shale and the angle of weak plane is obtained, as shown in Fig. 6. As the increase of confining pressure, the rock strength is increased. As the change of dip angle of weak plane, weak plane failure occurs and the strength decreases obviously at some dip angles. However, under high confining pressure, the angle range is reduced along the weak plane, indicating that the rock is not easily destroyed along the weak plane under high confining pressure.
In Fig. 7, the method for obtaining the anisotropic strength parameters of shale is as follows: the grid search calculation (Ma et al. 2020a, b, c, d) is carried out on the test results of shale samples with sampling angles of 0 and 90°, and the optimal cohesion of shale matrix and the optimal internal friction angle of matrix are 10.00 MPa and 33°, respectively. The optimal cohesion and internal friction angle of the weak plane are 3 MPa and 27°, respectively, according to the grid search calculation of the test results of shale samples with sampling angle of 30°, 45° and 60°.
Anisotropic wellbore collapse pressure model of shale

Wellbore stress distribution model
Considering that shale is a hard brittle rock, the stress-strain curve shows strong elasticity and significant anisotropy, so shale can be assumed to be a transversely isotropic medium. Once the borehole is drilled, the secondary distribution of stress around the well will occur. In order to accurately calculate the stress distribution around the well in transversely isotropic medium, necessary coordinate transformation is adopted. Four coordinate systems are involved here (Fig. 8): geodetic rectangular coordinate system, stratigraphic rectangular coordinate system, in situ stress rectangular coordinate system and borehole rectangular coordinate system.
The secondary distribution stress is closely related to the in situ stress and wellbore pressure, including three components: the in situ stress component before the wellbore is drilled, the stress component caused by wellbore formation, and the stress component caused by wellbore fluid pressure. Among them, the in situ stress component is solved by coordinate transformation, which is to convert the in situ stress to borehole rectangular coordinate system; the other two components are solved by the analytical formula derived by Aadnøy (1987). Finally, the three stress components can be linearly superimposed, and the borehole stress distribution model under the borehole rectangular coordinate system of the transversely isotropic medium can be obtained.
where xx , yy , zz , xy , yz , xz are wellbore circumferential stress in rectangular coordinate system, MPa; xx0 , yy0 , zz0 , xy0 , yz0 , xz0 are component of in situ stress in borehole rectangular coordinates, MPa; a ij is the flexibility matrix coefficient of formation in borehole local coordinates, i = 1, 2, 3, j = 1, 2, ⋅ ⋅ ⋅, 6; � 1 (z 1 ), � 2 z 2 , � 3 z 3 are the three  In order to facilitate the calculation and solution of collapse pressure, the stress component in borehole cylindrical coordinate should be calculated first. According to the coordinate relationship shown in Fig. 8, through the transformation of the rotation axis formula, the wellbore stress component in cylindrical coordinates can be obtained as follows: where r , , z , r , z , rz are wellbore stress component in cylindrical coordinates, MPa; is well circumference angle, • .

Calculation method of collapse pressure
In order to determine whether the wellbore stress exceeds the shale strength, the maximum and minimum principal stresses around the wellbore need to be obtained. According to the analysis of wellbore stress state, the maximum and minimum principal stresses can be obtained as follows: = xx sin 2 + yy cos 2 − xy sin 2 z = zz r = −0.5 xx sin 2 + 0.5 yy sin 2 + xy cos 2 z = yz cos − xz sin rz = yz sin + xz cos The angle between the maximum principal stress and the borehole axis is as follows: After the maximum and minimum principal stresses of the wellbore are obtained, the wellbore stability can be calculated in combination with the anisotropic strength criterion. Due to the existence of weak plane in shale, comprehensive calculation needs to be carried out according to borehole trajectory, bedding occurrence, angle between maximum principal stress and borehole axis, etc. It can be used to judge whether the weak plane failure or the main body failure occurs in the shale. The angle between the maximum principal stress and the normal direction of bedding plane is as follows: where n is the normal vector of bedding plane; N is the direction vector of maximum principal stress in wellbore.
If the angle calculated by Eq. (11) satisfies 1 ≤ ≤ 2 , substitute the maximum principal stress and minimum principal stress of the wellbore calculated by Eq. (9) into Eq. (6) to judge the wellbore stability. The critical collapse pressure along bedding failure can be obtained by solving the nonlinear equation with p w . If the angle calculated by Eq. (11) does not satisfies 1 ≤ ≤ 2 , substitute the maximum principal stress and minimum principal stress of the wellbore calculated by Eq. (9) into Eq. (5) to judge the wellbore stability. The critical collapse pressure along body failure of shale can be obtained by solving the nonlinear equation with p w .

Analysis of wellbore collapse pressure of continental shale in Songliao Basin
For continental shale in Songliao Basin, the collapse pressure is closely related to formation occurrence, borehole trajectory, shale rock mechanical parameters, in situ stress and pore pressure conditions, etc. In order to analyze the influence of anisotropic characteristics on the wellbore instability of continental shale, based on the measured data in the block, the equivalent density of collapse pressure is calculated under isotropic, elastic anisotropic, elastic and strength anisotropic conditions. The basic parameters are shown in Table 3.
Calculation results of isotropic model Figure 9 shows the hemispherical projection of the wellbore trajectory to the collapse pressure calculated by the isotropic model. The center of the contour line of equivalent density of collapse pressure is located in the projection center of the hemisphere, which is approximately symmetrically distributed along the direction of in situ stress, and the equivalent density of overall collapse pressure is relatively low. The collapse pressure is the lowest when the deviation angle is 0 and the azimuth is 90° or 270°. At the moment, the equivalent density of collapse pressure is 1.11 g/cm 3 ; As the increase of the deviation angle, the collapse pressure gradually increases. When the deviation angle is 60° and the azimuth is 0 or 180°, the collapse pressure reaches the maximum value of 1.58 g/cm 3 . When the deviation angle is small, the borehole orientation has little effect on the collapse pressure. However, with the increase of deviation angle, the influence of borehole orientation on collapse pressure becomes more and more significant. The equivalent density of collapse pressure first decreases and then increases in the process of wellbore orientation from the minimum horizontal principal stress to the maximum horizontal principal stress. Figure 10 shows the hemispherical projection of the wellbore trajectory to the collapse pressure calculated by the elastic anisotropy model. Under the condition of elastic anisotropy, the hemispherical projection of collapse pressure is  Azimuth/°E quivalent density of collapse pressure (g/cm 3 ) Fig. 9 Hemispherical projection of collapse pressure of isotropic model no longer symmetrically distributed along the direction of in situ stress. The collapse pressure under elastic anisotropy is the lowest when the deviation angle is 0 ~ 30° and the azimuth is near 0 or 180° and the equivalent density of collapse pressure is 1.13 g/cm 3 . As the increase of the deviation angle, the collapse pressure gradually increases. When the deviation angle is 40° ~ 65°, the equivalent density of collapse pressure at azimuth angles of 50°, 130°, 230° and 310° reaches the maximum value of 1.71 g/cm 3 . Under certain specific orientations, the influence of well deviation angle on collapse pressure is not obvious which can be inclined here. When the borehole orientation is parallel to the direction of minimum horizontal in situ stress, the collapse pressure is the largest. In the process of borehole orientation from the minimum horizontal principal stress to the maximum horizontal principal stress, the equivalent density of collapse pressure will also decrease first and then increase. Figure 11 is a hemispherical projection of borehole trajectory to collapse pressure calculated by elastic and strength anisotropy model. When considering the elasticity and strength anisotropy at the same time, the calculated results of collapse pressure are exceedingly different from those of isotropic model and elastic anisotropic model. Therefore, when analyzing the collapse and instability behavior of continental shale in Songliao Basin, the elasticity and strength anisotropy of shale cannot be ignored. Taking into account the characteristics of elasticity and strength anisotropy, the highest collapse pressure appears in a horizontal well, and its borehole orientation is parallel to the direction of the minimum horizontal principal stress.

Calculation results of elastic and strength anisotropy model
The equivalent density of the maximum collapse pressure is 1.86 g/cm 3 . When the borehole orientation is parallel to the direction of minimum horizontal in situ stress, the collapse pressure increases with the increase of well deviation angle. When the borehole orientation is parallel to the maximum horizontal in situ stress direction, the collapse pressure first decreases and then increases with the increase of well deviation angle, and the minimum collapse pressure equivalent density is 1.64 g/cm 3 . It can also be found from Fig. 10 that the equivalent density of collapse pressure of horizontal wells under different borehole orientations is quite different. Therefore, in the process of actual drilling, the borehole trajectory should be optimized and the variation range of borehole orientation should be reduced as far as possible to ensure that the drilling fluid density window in the whole horizontal section satisfies the requirements of wellbore stability.

Field application
Wellbore instability frequently occurs in a shale oil block in the north of Songliao Basin during drilling, resulting in serious downhole accidents such as wellbore collapse and sticking. The shale in the area is a typical hard brittle shale. The conventional isotropic model is used to analyze the wellbore collapse, and the results are quite different from the actual situation. Therefore, the wellbore collapse of a horizontal well G101 in the block is studied. The well depth is 4082 m, the horizontal section length is 1800 m, the actual drilling depth is 3087 m, and the final horizontal section length is 1300 m. In the process of well design, the conventional isotropic wellbore stability analysis model is used for the drilling fluid density. Therefore, due to the wellbore collapse and instability in the horizontal section, it takes   Fig. 11 Hemispherical projection of collapse pressure of elastic and strength anisotropy model 43.46 days to deal with downhole complex accidents, and the time to deal with downhole complex accidents reaches 58.34% of the whole well. The wellbore stability analysis model considering the influence of shale anisotropy and the conventional wellbore stability analysis model are used for comparative analysis. It can be seen that the bedding plane has no obvious effect on the stability of vertical wellbore, and the formation collapse density is 1.02 g/cm 3 .The impact on horizontal wells is more serious. If the influence of bedding plane is not considered, the collapse density to maintain the stability of wellbore is 1.57 g/cm 3 . Considering the influence of bedding plane, the collapse density to maintain the stability of wellbore is 1.86 g/cm 3 . It is recommended to optimize the drilling orientation when designing horizontal wells in the block in the later stage. The influence of shale bedding on horizontal wells and deviated wells is very significant, especially when the well deviation angle is greater than 30° and the borehole orientation is within 90° of the minimum horizontal in situ stress direction. In the drilling of shale formation in the vertical section of well G101, the density of 1.12 ~ 1.30 g/cm 3 is adopted according to the design. Wellbore collapse was not found during drilling, indicating that the analysis results of shale wellbore collapse pressure in vertical section are in good agreement with the actual situation. In the drilling of shale formation in horizontal well section, the density is 1.58 ~ 1.76 g/cm 3 according to the design. The drilling fluid is much lower than the collapse density of 1.86 g/cm 3 calculated by the anisotropic model in this paper. During this period, serious wellbore instability occurred, and the problem of wellbore collapse was not alleviated until the drilling fluid density was increased from 1.68 to 1.94 g/cm 3 after 2876 m. The results of wellbore stability risk analysis are basically consistent with the actual situation.

Summary and Conclusions
(1) The mineral composition test, uniaxial and triaxial compressive strength tests of shale in the northern block of Songliao Basin are carried out. The hard and brittle properties of shale in this block are defined, and the anisotropic elasticity and strength parameters of shale oil block in the north of Songliao basin are determined. The experimental results show that there are significant differences in shale strength at different sampling angles, and the shale strength is the lowest when the sampling angle is 60°. The failure modes of rock samples under triaxial compression include matrix shear failure and weak plane shear failure. The anisotropic elastic parameters of shale are as follows: the elastic modulus of parallel bedding plane is 11.258 GPa, the Poisson's ratio of parallel bedding plane is 0.1875, the elastic modulus of vertical bedding plane is 3.373 GPa, and the Poisson's ratio of vertical bedding plane is 0.2445. The anisotropic strength parameters of shale are as follows: the cohesion of shale matrix is 10.0 MPa, the internal friction angle of shale matrix is 33°, the cohesion of weak plane is 3 MPa, and the internal friction angle of weak plane is 27°.
(2) According to the anisotropic characteristics of shale in the north of Songliao Basin, a calculation model of wellbore collapse pressure considering the anisotropy of shale elasticity and strength is established. Based on the measured rock mechanics and geomechanics parameters, the wellbore collapse pressures under isotropy, elastic anisotropy, elasticity and strength anisotropy are analyzed. The results show that the collapse pressure is the highest when elastic and strength anisotropy are taken into account, and the effect of strength anisotropy is obviously greater than that of elastic anisotropy. While the deviation angle is small, the borehole orientation has little effect on the collapse pressure. With the increase of well deviation angle, the influence of borehole orientation on collapse pressure becomes more and more significant. The stability of drilling horizontal wells along the direction of horizontal maximum principal stress is the best, followed by vertical wells and small angle directional wells, and the stability of drilling horizontal wells along the direction of horizontal maximum principal stress is the worst. (3) Taking a horizontal well G101 in the north of Songliao Basin as an example, the analysis shows that the analysis results are basically consistent with the actual situation by using the wellbore collapse instability analysis method considering the effects of shale elasticity and strength anisotropy. Therefore, it is recommended to design the drilling fluid density in the shale formation of vertical well section according to the conventional wellbore stability analysis method. In the shale formation of inclined well section and horizontal section, the method described in this paper is recommended, that is, the drilling fluid density is designed for the wellbore collapse instability analysis method considering the influence of shale elasticity and strength anisotropy. So as to reduce the risk of wellbore collapse and instability and ensure safe and rapid drilling. (4) The advantage of this paper is that the wellbore stability model of layered shale anisotropy is established, and the dual effects of shale elastic anisotropy and strength anisotropy are considered. The disadvantage is that this paper mainly analyzes from the perspective of mechanics, without considering chemistry and thermology. The angle between the failure plane of shale body and the maximum principal stress, • The angle between the normal of weak plane and the maximum principal stress, •