Effects of Hydration during Drilling on Fracability of Shale Oil Formations: A Case Study of Da’anzhai Section Reservoir in Sichuan Basin, China

: Accurate evaluation of shale oil reservoir fracability helps avoid blind fracturing and ensures efﬁcient fracturing. However, the current evaluation of the fracability index rarely considers the impact of hydration caused by drilling ﬂuid invasion during drilling. The results of a rock triaxial mechanical test conducted to evaluate the mechanical properties of shale oil reservoirs are reported in this paper. Based on the results, we developed a comprehensive evaluation method of shale oil reservoir fracability that considers hydration; the effects of the brittleness index, horizontal difference stress, and fracture toughness; and the law of water phase intrusion into shale oil reservoirs. The research results show that the average compressive strength decreased by 37.99%, the average elastic modulus decreased by 53.36%, and the average Poisson’s ratio increased by 68.75% after being soaked for 48.00 h at 80 ◦ C and 30.00 MPa. The water saturation rate at the borehole wall was the highest; with the extension to the stratum, it gradually decreases to the original water saturation rate of the formation, while the ﬂuctuation radius gradually increases with time. The Young’s modulus and fracture toughness decrease, the Poisson’s ratio increases, and the fracability index reaches a maximum value at the wellbore (i.e., the highest water saturation rate), indicating that the strength of the hydrated rock decreases and it can


•
The compressive strength of dry rock samples and soaked rock samples were tested by triaxial mechanical experiments. • According to the stress-strain curves, the changes of compressive strength, elastic modulus and Poisson's ratio before and after hydration were analyzed.
2. How to establish a mathematical model of shale fracability that comprehensively considers brittleness index, horizontal differential stress, and fracture toughness is proposed?
• Fracability index is proportional to brittleness index and horizontal differential stress, and inversely proportional to fracture toughness. • Young 's modulus, Poisson's ratio and tensile strength in fracability index model are related to rock water saturation, so the fracturing index is also related to rock hydration. 3. What is the effect of shale hydration on fracability and fracture pressure?
• As the water saturation rate of the rock gradually decreases, the FI of the formation rock de-creases • Combined with the variation law of water saturation rate, it shows that after the shale around the well is hydrated, the fracture pressure of the formation decreases Abstract: Accurate evaluation of shale oil reservoir fracability helps avoid blind fracturing and ensures efficient fracturing. However, the current evaluation of the fracability index rarely considers the impact of hydration caused by drilling fluid invasion during drilling. The results of a rock triaxial mechanical test conducted to evaluate the mechanical properties of shale oil reservoirs are reported in this paper. Based on the results, we developed a comprehensive evaluation method of shale oil reservoir fracability that considers hydration; the effects of the brittleness index, horizontal difference stress, and fracture toughness; and the law of water phase intrusion into shale oil reservoirs. The research results show that the average compressive strength decreased by 37.99%, the average elastic modulus decreased by 53.36%, and the average Poisson's ratio increased by 68.75% after being soaked for 48.00 h at 80 • C and 30.00 MPa. The water saturation rate at the borehole wall was the highest; with the extension to the stratum, it gradually decreases to the original water saturation rate of the formation, while the fluctuation radius gradually increases with time. The Young's modulus and fracture toughness decrease, the Poisson's ratio increases, and the fracability index reaches a maximum value at the wellbore (i.e., the highest water saturation rate), indicating that the strength of the hydrated rock decreases and it can be easily fractured. The case analysis shows that the optimal fracturing position of the Da'anzhai Section of Well NC2H is around 2600 m deep. After the

Introduction
In 2018, Sinopec estimated that the technically recoverable resources of shale oil in China were (74-372)10 8 t. The continental shale oil and gas resources in the Sichuan Basin accounted for 30% of the national shale oil and gas geological reserves and 40% of the technically recoverable reserves. Shale oil reservoirs generally do not have natural productivity, while their fracability depends on the degree of natural fracture development, the mineral composition of the rock [1] (Fathy 2021), and the mechanical properties of the rock [2,3]  Studies have shown that the reservoir space type of the shale reservoirs in the Da'anzhai Section of the Ziliujing Formation in the Yuanba and Fuling areas is primarily inorganic pores, and the shale porosity ranges from 1.5% to 6.7%, with an average of 3.9% [4] (Zhou 2020). Affected by the drilling pressure difference and chemical potential, the drilling fluid intrudes into the rock along the pores and bedding fractures, reducing the strength of the rock body and weak plane. The concentration of clay minerals in the shale is relatively high; illite, montmorillonite, and montmorillonite-mixed layers have strong hydration ability. Under the influence of hydration stress, ion transport, film effect, and solute diffusion, clay minerals are easily hydrated, thus reducing the rock strength [5][6][7] (Muniz 2005; van Oort E 1996; Wang 2006). Many studies have attempted to describe this phenomenon: Yew and Chenevert [8] (1990) proposed a model for calculating the amount of adsorbed water after quantitatively analyzing the mechanochemical coupling. Hale and Mody [9] (1993) pointed out that the chemical kinetic factors that cause rock deformation and failure mainly include [10] (2011) van der Waals force, electrostatic attraction, plate potential energy, capillary tension, matrix suction, osmotic pressure, and concentration polarization boundary resistance. He [11] (2020) explored the deformation and failure law of rational shale under the action of mechanical-chemical coupling. Therefore, the shale mechanical-chemical coupling methods can be divided into the following: water transport mechanism in shale, equivalent pore pressure method, total adsorbed water correlation method, and incremental elasticity theory of total water potential [12] (Deng 2003).
To study shale fracability, scholars have established different evaluation methods based on different physical parameters and mathematical theories. Many scholars [13][14][15][16][17][18] [19,20] (2013; 2019) believed that the brittleness index and the plasticity coefficient should both be used to evaluate the shale fracability. To this end, they carried out numerous experiments and proposed a rock-brittleness evaluation method based on triaxial compression energy dissipation. Mike J. Mullen [21] (2012) deemed that the fracability should consider reservoir sedimentary characteristics and mineral content composition. Chen [22] (2001) evaluated the fracability of shale reservoirs in terms of rock mechanics. Yuan [17] (2013) established a fracability evaluation model using three parameters: the elastic modulus, Poisson's ratio, and uniaxial tensile strength. Zhao [23] (2015) used shale brittleness and fracture toughness to evaluate shale fracability. Sui [24] (2016) discussed how fracability is affected by the clay mineral content, cohesion, and internal friction angle. Yue [25] (2019) established a fuzzy comprehensive evaluation model of fracability based on the fuzzy mathematics theory and the analytic-hierarchy-process-entropy method. Ma [26] (2019) focused on the lacustrine shale reservoirs located in the Dongying Sag and the Subei Basin. The results showed that the macrostructure, mineral content, lithology, and the presence of microfractures in the shale are the main factors affecting its brittleness. Considering the fracture toughness of carbonate minerals and the plastic effect of organic matter, a formula for calculating the brittleness index based on the mineral composition method was proposed. Jia [27] (2020) explored the mineralogical composition, microstructure formation, static elastic properties, fracture toughness, and mechanical anisotropy of the Longmaxi shale through a series of laboratory geomechanical and petrophysical experiments. Four factors, namely brittleness, quartz content, diagenesis, and natural fracture, were analyzed to determine the fracability of the shale in the Lower Cambrian Niutitang Formation, Northwestern Hunan [28] (Wu 2018). Brittle rocks can have higher fracture toughness, making their fracture more difficult. Iyare [29] (2022) used an existing fracability evaluation model describing the brittleness, fracture toughness, and minimum horizontal in situ stress to evaluate the fracability of the four lithofacies within the Naparima Hill Formation. Brittleness indices based on the rock elastic properties, rock strength, and mineral concentrations were used to evaluate the rock brittleness, while the fracture toughness was determined from the tensile strength.
However, the previous fracability evaluation methods considered a single factor or a simple superposition of multiple factors, failing to capture the physical and chemical reactions between the working fluid and shale. Studies have shown that shale contains a larger number of clay minerals and has a strong hydration ability. After hydration, the rock strength of shale will be greatly reduced, and creep will occur simultaneously, changing the crack structure. Therefore, the chemical interaction between the working fluid and the shale cannot be ignored. By considering the influence of hydration, the fracturability of shale oil reservoirs can be characterized more comprehensively.
The research in this paper mainly focuses on the following three parts: (1) Through experiments, the changes in the rock mechanical properties of the shale oil reservoirs before and after soaking are analyzed, and the changing laws of rock mechanical parameters are recorded. (2) We analyze the law of water-based drilling fluid invading shale oil reservoirs under the action of chemical potential. (3) Considering the rock-brittleness index, horizontal difference stress, and fracture toughness, a set of fracability evaluation methods suitable for shale oil reservoirs is formed. (4) Taking Well NC2H as an example, we predicted the fracability index of the Da'anzhai Section at different depths, selected the most suitable sweet spot position for fracturing, and calculated the corresponding fracturing pressure.

Shale Mechanics Experiment
Young's modulus and Poisson's ratio are the principal rock mechanical parameters that are used to characterize shale brittleness [13] (RICKMAN 2008). The Young's modulus reflects the ability of shale to maintain fractures after being fractured, while the Poisson's ratio reflects the ability of shale to fracture under pressure. The higher the Young's modulus of shale and the lower the Poisson's ratio, the stronger the brittleness is. By comparing the changes in the mechanical parameters of the shale before and after soaking, we shed further insight into the effect of hydration on rock fracability.

Sample Preparation
In this experiment, eight standard rock samples were prepared and used to test the mechanical properties of the shale before and after soaking, labeled A, B, C, D, E, F, G, and H, and the size of the rock samples was Φ 25 mm × 50 mm, where A, B, C, and D are dry rock samples, and four rock samples, E, F, G, and H, are soaked. The soaking conditions are as follows: pressure of 30 MPa, time of 48 h, and temperature of 80 • C.

Methodology
To study the mechanical parameters of the Da'anzhai shale formation in Central Sichuan in the downhole stress environment, the triaxial mechanical test was carried out. The GCTS triaxial mechanical testing machine was used to test and analyze the triaxial compressive strength, elastic modulus, and Poisson's ratio of the Da'anzhai formation shale in the downhole stress environment. Experimental scheme: First, a triaxial mechanical test was performed on four groups of dry rock samples (A, B, C, and D); the experimental confining pressure was set to 30 MPa, and the triaxial compressive strength, elastic modulus, and Poisson's ratio of the dry rock samples were calculated according to the stress-strain curve. Similarly, triaxial mechanical tests were performed on four groups of soaked rock samples (E, F, G, and H), and the experimental confining pressure was set to 30 MPa. According to the stress-strain curve, the triaxial compressive strength, elastic modulus, and Poisson's ratio of the soaked shale were calculated. Finally, the experimental data were counted, and the test results were compared and analyzed.

Test Results and Analysis
The results of the triaxial mechanical experiments are shown in Figure 1. stress-strain curve. Similarly, triaxial mechanical tests were performed on four groups of soaked rock samples (E, F, G, and H), and the experimental confining pressure was set to 30 MPa. According to the stress-strain curve, the triaxial compressive strength, elastic modulus, and Poisson's ratio of the soaked shale were calculated. Finally, the experimental data were counted, and the test results were compared and analyzed.

Test Results and Analysis
The results of the triaxial mechanical experiments are shown in Figure 1.  When comparing the dry rock sample and the soaked rock sample, we found that, under the same confining pressure, the average compressive strength of the dry rock sample was 110.50 MPa, the average elastic modulus was 15.719 GPa, and the average Poisson's ratio was 0.16. The average compressive strength of the soaked rock samples was 68.52 MPa, a decrease of 37.99%; the average elastic modulus was 7.332 GPa, a decrease of 53.36%; and the average Poisson's ratio was 0.27, an increase of 68.75%. After soaking, the compressive strength and elastic modulus of the shale decreased sharply, while the Poisson's ratio increased, indicating that the mechanical strength and elastic parameters of the shale rock in the Da'anzhai Section were obviously affected by the drilling fluid soaking. This also showed that the water saturation rate of the shale also changed, causing changes in formation fracability.
The triaxial test curves of rock samples A, B, G, and H are shown in Figure 2. Ev is the vertical strain, Ea is the volumetric strain, and Er is the radial strain. The figure shows that the compressive strength of the shale after soaking is significantly reduced, the rock brittleness is reduced, the plasticity is enhanced, and the residual strength can still be maintained after failure. When comparing the dry rock sample and the soaked rock sample, we found that, under the same confining pressure, the average compressive strength of the dry rock sample was 110.50 MPa, the average elastic modulus was 15.719 GPa, and the average Poisson's ratio was 0.16. The average compressive strength of the soaked rock samples was 68.52 MPa, a decrease of 37.99%; the average elastic modulus was 7.332 GPa, a decrease of 53.36%; and the average Poisson's ratio was 0.27, an increase of 68.75%. After soaking, the compressive strength and elastic modulus of the shale decreased sharply, while the Poisson's ratio increased, indicating that the mechanical strength and elastic parameters of the shale rock in the Da'anzhai Section were obviously affected by the drilling fluid soaking. This also showed that the water saturation rate of the shale also changed, causing changes in formation fracability.
The triaxial test curves of rock samples A, B, G, and H are shown in Figure 2. E v is the vertical strain, E a is the volumetric strain, and E r is the radial strain. The figure shows that the compressive strength of the shale after soaking is significantly reduced, the rock brittleness is reduced, the plasticity is enhanced, and the residual strength can still be maintained after failure. soaking. This also showed that the water saturation rate of the shale also changed, causing changes in formation fracability.
The triaxial test curves of rock samples A, B, G, and H are shown in Figure 2. Ev is the vertical strain, Ea is the volumetric strain, and Er is the radial strain. The figure shows that the compressive strength of the shale after soaking is significantly reduced, the rock brittleness is reduced, the plasticity is enhanced, and the residual strength can still be maintained after failure. Stress/MPa

Fracability Evaluation Method Considering Hydration
First, the three indexes, namely brittleness index, in situ stress, and fracture toughness, were selected to establish a fracability evaluation model. Then the mathematical model of water phase intrusion into shale oil reservoir was deduced, and both the law of influence of hydration on the brittleness index (Young's modulus and Poisson's ratio) and fracture toughness were obtained. Finally, the effect of hydration on the fracability index was calculated.

Fracability Index Model
Fracability reflects the shale geology and reservoir characteristics, while it is affected by many factors. According to the research experience of scholars on the fracability of shale gas reservoirs, the brittleness index represents the difficulty of shale fracturing. At present, the commonly used calculation methods of the brittleness index include the Young's modulus and Poisson's ratio [13] (RICKMAN 2008). This method makes it easy to obtain data, and it is more closely integrated with the field and more widely used. In this paper, the elastic modulus of rock refers to the Young's modulus for calculation: where Brit is the brittleness index, 0-1; YM is the static Young's modulus; YMcmax and YMcmin are the maximum and minimum static Young's modulus, respectively, GPa; YMBRIT is the normalized Young's modulus, 0-1; PR is the static Poisson's ratio; PRcmax and PRcmin are, respectively, the largest and smallest static Poisson's ratio in the region, dimensionless; and PRBRIT is the normalized Poisson's ratio, 0-1.

Fracability Evaluation Method Considering Hydration
First, the three indexes, namely brittleness index, in situ stress, and fracture toughness, were selected to establish a fracability evaluation model. Then the mathematical model of water phase intrusion into shale oil reservoir was deduced, and both the law of influence of hydration on the brittleness index (Young's modulus and Poisson's ratio) and fracture toughness were obtained. Finally, the effect of hydration on the fracability index was calculated.

Fracability Index Model
Fracability reflects the shale geology and reservoir characteristics, while it is affected by many factors. According to the research experience of scholars on the fracability of shale gas reservoirs, the brittleness index represents the difficulty of shale fracturing. At present, the commonly used calculation methods of the brittleness index include the Young's modulus and Poisson's ratio [13] (RICKMAN 2008). This method makes it easy to obtain data, and it is more closely integrated with the field and more widely used. In this paper, the elastic modulus of rock refers to the Young's modulus for calculation: where B rit is the brittleness index, 0-1; YM is the static Young's modulus; YM cmax and YM cmin are the maximum and minimum static Young's modulus, respectively, GPa; YM BRIT is the normalized Young's modulus, 0-1; PR is the static Poisson's ratio; PRcmax and PRcmin are, respectively, the largest and smallest static Poisson's ratio in the region, dimensionless; and PR BRIT is the normalized Poisson's ratio, 0-1.
The research results show that, in addition to overcoming the tensile strength of the rock body, the net pressure in the fracture also must overcome the horizontal stress difference. This should enable the fracture to open and extend, forming a complex fracture network [30,31] (Chen 2020; Lian 2022). The horizontal difference stress coefficient, K h , is expressed as follows: where K h is the horizontal difference stress coefficient, dimensionless; σ H is the maximum horizontal in situ stress, MPa; and σ h is the horizontal minimum ground stress, MPa. Fully determining the fracability of the reservoir does not rely only on the elastic modulus and Poisson's ratio. For example, the elastic modulus and Poisson's ratio of some rocks are similar, but the brittleness greatly differs. The important difference is the fracture toughness. Fracture toughness is crucial for characterizing the difficulty of reservoir fracturing, which reflects the ability to maintain the forward extension of the fracture after its formation.
The type I fracture toughness (K Ic ) and type II fracture toughness (K IIC ) of shale can be calculated from the confining pressure and uniaxial tensile strength (S t ) on the fracture surface [22,32,33] (Chen 1997;Jin 2001Jin ,2011. The calculation methods are as shown in Formulas (3) and (4): where σ n is the confining pressure on the crack surface, MPa; S t is the uniaxial tensile strength, MPa; and K Ic is the fracture toughness, MPa·m 0.5 .
The rock mechanical parameters involved in the calculation method of types I and II fracture toughness, such as uniaxial tensile strength and confining pressure, can be obtained from logging data. The calculation method is as follows: where V cl is the clay volume, dimensionless; E d is the dynamic elastic modulus, GPa; and K is a constant, K = 12.26: To determine the confining pressure, the confining pressure acting on the fracture surface of the n-th branch fracture of shale gas fracturing can be expressed as [34] follows: where n = [cos βsin θ, sin βcos θ, sin β]; β is the fracture dip angle; θ is the angle between the fracture strike and the well axis; and σ v is the overburden pressure, MPa; According to the above studies, the fracability of shale gas reservoirs is obviously positively correlated with the brittleness index and negatively correlated with the fracture toughness. Considering the brittleness index, in situ stress, and rock fracture toughness, a new model for fracability evaluation is established in this paper, as shown in Equation (8). where FI is the fracability index; a and b are constant coefficients, dimensionless. According to the classical rock mechanics theory, the fracture of the rock during hydraulic fracturing is caused by the excessive density of the fracturing fluid in the well, which causes the stress on the rock to exceed the tensile strength. The fracture initiation pressure when the rock produces tensile failure is as follows [35]  ): When the formation permeability is low, δ equals 0, as shown with the following equation:

Analysis of the Law of Water Phase Intrusion into Shale Oil Reservoirs
During the drilling of shale oil reservoirs, the pressure difference and capillary force will jointly affect the entry of the water-based drilling fluid into the oil layer and change the water saturation rate of the formation rocks around the well. Therefore, for different shale oil reservoirs, different curve characteristics and different J 0 functions will be obtained according to the experiment [36] (Buckley 1942). The expression of the J 0 function is as follows: where σ ow is the surface tension of the oil-water interface, N m; θ is the contact angle of the oil-water interface, • ; p c (S) is the function of the capillary force with saturation; and ϕ is the shale oil reservoir porosity, %. According to the J 0 function, the expression of the capillary force function can be derived as follows [37] (Jiang 2014): According to Formula (12), the capillary force is related to formation wettability. When the size of the wetting contact angle changes, the positive and negative capillary force can be directly changed; that is, the wetting reversal occurs. When a water-based fluid is used to drill a shale oil reservoir, the fluid flow near the wellbore is an oil-water two-phase flow process. Therefore, the model assumes the following conditions:
During drilling fluid intrusion into shale oil layer, the fluid viscosity and volume coefficient remain unchanged; 3.
During the whole two-phase flow process, the gravity effect is not considered, and the oil and water are immiscible when in contact with each other.
Based on the above model assumptions, according to Darcy's law of oil-water twophase flows, we obtain the following: (13) where K rw is the relative permeability of the water phase of shale oil reservoir; v w is the water phase velocity, m/s; µ w is the water phase viscosity, mPa·s; p w is the water phase pressure, Pa; K ro is Oil phase relative permeability, dimensionless; v o is the oil phase velocity, m/s; µ o is the oil phase viscosity, mPa·s; and p o is oil phase pressure, Pa. In the process of drilling fluid intrusion into the reservoir, the oil-water two-phase continuity equation is as follows: where is the Laplace operator; K w describes the water permeability of the shale oil reservoir, mD; S w is the water saturation rate, %; K o describes the oil permeability of the shale oil reservoir, mD; S o is the oil saturation, %; and t is the oil-water two-phase seepage time, s.
If the oil phase and water phase volume fracability factors are considered, the drilling fluid invasion model of the oil layer can be obtained: In the above model, to combine water phase pressure and oil phase pressure, the water saturation rate and oil saturation, and the absolute permeability and relative permeability of the two phases, the following three auxiliary equations are required [38] (Naseri 2012): where A, B, and C are all constants which can be obtained from experiments according to the reservoir characteristics. In this paper, the empirical values of a certain oilfield are taken: A = 0.9, B = 0.8, and C = 4; S or and S wi are the residual oil saturation and original water saturation rate of the reservoir, respectively, dimensionless. According to Equations (16)- (18), there are a total of six dependent variables, namely S w , S o , p w , p o , v w , and v o . The capillary force p c (S w ) can be obtained according to the J 0 function. This set of equations can be solved if the initial and boundary conditions are given.
Initial time: t = 0, r w < r < r e , S w = S wi ; Inner boundary conditions of shale oil reservoirs: External boundary conditions of shale oil reservoirs: where ρ b is drilling fluid pressure in wellbore, MPa; ρ e is original oil reservoir pressure, MPa. The above equations are strongly nonlinear and cannot be analyzed by using conventional analytical methods. Therefore, in this paper, the staggered grid is used to conduct double implicit difference on the above equations. The following differential equation of water-based drilling fluid intrusion into the shale oil reservoir model is obtained: Processes 2022, 10, 2313 9 of 16 According to the mass conservation theory, the formula of water content in the shale oil layer is theoretically derived: where ρ wf and ρ s are, respectively, the density of the drilling fluid and shale framework, g/cm 3 ; V W and V S are, respectively, the volume of water in the rock and the volume of the shale framework, cm 3 . Scholars have conducted much experimental research on water absorption [39] (Huang 1995) and concluded that the elastic modulus and Poisson's ratio change regularly with the water saturation rate of shale: where E is the elastic modulus of rock after water absorption, MPa; E 1 and E 2 are coefficients, which can be obtained from experiments; the empirical values are 4 × 10 4 and −11; W and W a are the mass percentage of adsorbed water and the mass percentage of original water content, %; and ν is the Poisson's ratio of rock formation after water absorption. According to Equations (20) and (21) and data listed in Table 1, from the initial water saturation rate, S wi = 0.25; the residual oil saturation, S or = 0.2; the porosity, ϕ = 0.25; and the oil layer permeability, K = 1.0 mD, obtained Figure 3, which shows the variation law of rock water saturation rate, S w (dimensionless), with time and dimensionless distance from the borehole axis (r/r w ). Note that, for the convenience of calculations and expressions, the dimensionless distance is used to represent the distance between the calculated formation rock position and the borehole axis, which is represented by r/r w . According to the mass conservation theory, the formula of water content in the shale oil layer is theoretically derived: where wf and are, respectively, the density of the drilling fluid and shale framework, g/cm 3 ; VW and VS are, respectively, the volume of water in the rock and the volume of the shale framework, cm 3 .
Scholars have conducted much experimental research on water absorption [39] (Huang 1995) and concluded that the elastic modulus and Poisson's ratio change regularly with the water saturation rate of shale: where E is the elastic modulus of rock after water absorption, MPa; E1 and E2 are coefficients, which can be obtained from experiments; the empirical values are 4 × 10 4 and −11; W and Wa are the mass percentage of adsorbed water and the mass percentage of original water content, %; and is the Poisson's ratio of rock formation after water absorption. According to Equations (20) and (21) and data listed in Table 1, from the initial water saturation rate, Swi = 0.25; the residual oil saturation, Sor = 0.2; the porosity, φ = 0.25; and the oil layer permeability, K = 1.0 mD, obtained Figure 3, which shows the variation law of rock water saturation rate, Sw (dimensionless), with time and dimensionless distance from the borehole axis (r/rw). Note that, for the convenience of calculations and expressions, the dimensionless distance is used to represent the distance between the calculated formation rock position and the borehole axis, which is represented by r/rw.  In Figure 3, the water saturation rate at the wellbore is the highest, and with the extension to the formation, the water saturation rate gradually decreases to the original water saturation rate of the formation. As time passes, the fluctuation radius (dimensionless In Figure 3, the water saturation rate at the wellbore is the highest, and with the extension to the formation, the water saturation rate gradually decreases to the original water saturation rate of the formation. As time passes, the fluctuation radius (dimensionless distance corresponding to the change in the formation water saturation rate) gradually increases. After 350 h, the fluctuation radius increases to r/r w = 38. As the water saturation rate changes, and according to Equations (21) and (22), the Young's modulus and Poisson's ratio of the rock will change, and this, in turn, affects the rock's tensile strength (Equation (5)), fracture toughness (Equations (3) and (4)), and fracture pressure (Equations (9) and (10)).

Variation of Rock Parameters and Fracability after Hydration
According to the calculation results of the water saturation rate, combined with the fracability model, the variation law of each parameter in the FI was calculated. First, the variation laws of YM and µ with dimensionless radial distance at different times were calculated. The results are shown in Figures 4 and 5, respectively.
Processes 2022, 10, x FOR PEER REVIEW 10 of 16 distance corresponding to the change in the formation water saturation rate) gradually increases. After 350 h, the fluctuation radius increases to r/rw = 38. As the water saturation rate changes, and according to Equations (21) and (22), the Young's modulus and Poisson's ratio of the rock will change, and this, in turn, affects the rock's tensile strength (Equation (5)), fracture toughness (Equations (3) and (4)), and fracture pressure (Equations (9) and (10)).

Variation of Rock Parameters and Fracability after Hydration
According to the calculation results of the water saturation rate, combined with the fracability model, the variation law of each parameter in the FI was calculated. First, the variation laws of YM and μ with dimensionless radial distance at different times were calculated. The results are shown in Figures 4 and 5, respectively.  According to the data in Figure 4, the Young's modulus also changed with the change of water saturation rate. Near the wellbore with the highest water saturation rate, the Young's modulus demonstrated the lowest value. As the radial distance increases, the Young's modulus gradually increases and approaches the original value at the edge of the fluctuation radius. By comparing the maximum and the minimum values of the Young's modulus, we found that the change in water saturation rate greatly affects the Young's modulus of rock, from 15.83 GPa to 5.68 GPa, with a decrease of 64.6%.
As shown in Figure 5, with the increase in the water saturation rate, the Poisson's ratio of the rock gradually increases. The original Poisson's ratio of the rock is 0.226, and after hydration, it demonstrates a maximum value of 0.347 near the wellbore, an increase of 53.5%. distance corresponding to the change in the formation water saturation rate) gradually increases. After 350 h, the fluctuation radius increases to r/rw = 38. As the water saturation rate changes, and according to Equations (21) and (22), the Young's modulus and Poisson's ratio of the rock will change, and this, in turn, affects the rock's tensile strength (Equation (5)), fracture toughness (Equations (3) and (4)), and fracture pressure (Equations (9) and (10)).

Variation of Rock Parameters and Fracability after Hydration
According to the calculation results of the water saturation rate, combined with the fracability model, the variation law of each parameter in the FI was calculated. First, the variation laws of YM and μ with dimensionless radial distance at different times were calculated. The results are shown in Figures 4 and 5, respectively.  According to the data in Figure 4, the Young's modulus also changed with the change of water saturation rate. Near the wellbore with the highest water saturation rate, the Young's modulus demonstrated the lowest value. As the radial distance increases, the Young's modulus gradually increases and approaches the original value at the edge of the fluctuation radius. By comparing the maximum and the minimum values of the Young's modulus, we found that the change in water saturation rate greatly affects the Young's modulus of rock, from 15.83 GPa to 5.68 GPa, with a decrease of 64.6%.
As shown in Figure 5, with the increase in the water saturation rate, the Poisson's ratio of the rock gradually increases. The original Poisson's ratio of the rock is 0.226, and after hydration, it demonstrates a maximum value of 0.347 near the wellbore, an increase of 53.5%. According to the data in Figure 4, the Young's modulus also changed with the change of water saturation rate. Near the wellbore with the highest water saturation rate, the Young's modulus demonstrated the lowest value. As the radial distance increases, the Young's modulus gradually increases and approaches the original value at the edge of the fluctuation radius. By comparing the maximum and the minimum values of the Young's modulus, we found that the change in water saturation rate greatly affects the Young's modulus of rock, from 15.83 GPa to 5.68 GPa, with a decrease of 64.6%.
As shown in Figure 5, with the increase in the water saturation rate, the Poisson's ratio of the rock gradually increases. The original Poisson's ratio of the rock is 0.226, and after hydration, it demonstrates a maximum value of 0.347 near the wellbore, an increase of 53.5%.
According to the calculation model of fracture toughness (Formulas (3) and (4)), the calculation of fracture toughness is closely related to the tensile strength and elastic modulus of the rock. In this paper, the change law of rock fracture toughness is calculated and analyzed. The results are shown in Figure 6.
Processes 2022, 10, x FOR PEER REVIEW 11 of 16 According to the calculation model of fracture toughness (Formulas (3) and (4)), the calculation of fracture toughness is closely related to the tensile strength and elastic modulus of the rock. In this paper, the change law of rock fracture toughness is calculated and analyzed. The results are shown in Figure 6. Owing to the change in water saturation rate, fracture toughness also exhibits a similar variation law to water saturation rate and Young's modulus. In the near wellbore zone, the rock has a high water saturation rate and low fracture toughness. With the extension from the periphery of the well to the formation, the water saturation rate gradually decreased to the original water saturation rate of the formation, while the fracture toughness also gradually returned to the unhydrated state of the rock.
Based on the above parameters and according to the fracability model, the calculation and analysis of the FI were carried out. The results are shown in Figure 7. r/r w =100 r/r w =200 r/r w =300 r/r w =400 Figure 6. Variation of fracture toughness with dimensionless distance and time ((left) type I fracture toughness and (right) type II fracture toughness).
Owing to the change in water saturation rate, fracture toughness also exhibits a similar variation law to water saturation rate and Young's modulus. In the near wellbore zone, the rock has a high water saturation rate and low fracture toughness. With the extension from the periphery of the well to the formation, the water saturation rate gradually decreased to the original water saturation rate of the formation, while the fracture toughness also gradually returned to the unhydrated state of the rock.
Based on the above parameters and according to the fracability model, the calculation and analysis of the FI were carried out. The results are shown in Figure 7.
Processes 2022, 10, x FOR PEER REVIEW 11 of 16 According to the calculation model of fracture toughness (Formulas (3) and (4)), the calculation of fracture toughness is closely related to the tensile strength and elastic modulus of the rock. In this paper, the change law of rock fracture toughness is calculated and analyzed. The results are shown in Figure 6. Owing to the change in water saturation rate, fracture toughness also exhibits a similar variation law to water saturation rate and Young's modulus. In the near wellbore zone, the rock has a high water saturation rate and low fracture toughness. With the extension from the periphery of the well to the formation, the water saturation rate gradually decreased to the original water saturation rate of the formation, while the fracture toughness also gradually returned to the unhydrated state of the rock.
Based on the above parameters and according to the fracability model, the calculation and analysis of the FI were carried out. The results are shown in Figure 7. r/r w =100 r/r w =200 r/r w =300 r/r w =400 Figure 7. Variation of fracability index with dimensionless distance and time. Figure 7 shows that the change in water saturation rate changes rock fracability. As the drilling fluid invades the rock, the hydration becomes gradually obvious. The initial FI of the rock at r/r w = 10 is 0.39. After hydration occurs, the FI rises rapidly to 0.71. This shows that hydration has a significant effect on the FI. The rock strength after hydration is reduced and is more prone to fracture. From Figure 3, we observe that, as the dimensionless distance increases, the water saturation rate of the rock gradually decreases, and the FI of the formation rock gradually decreases to the initial value.

Case Study
The shale formation in the Da'anzhai Section in Central Sichuan has favorable conditions for lacustrine shale oil and gas accumulation; its high-quality shale is mainly developed in the second sub-member, with favorable conditions such as high organic carbon content, good physical properties, stable distribution, and good oil and gas properties. For example, NC2H is the first shale oil horizontal well in the Sichuan Basin. The depth of the Da'anzhai Section ranges from 2598 to 2703 m.
The method of calculating rock-mechanics parameters from well-logging data has been applied maturely in petroleum engineering. The calculation methods are as shown in Formula (23) [17] (Yuan 2013): where v p is the longitudinal wave velocity, m/s; v s is the shear wave velocity, m/s; γ max and γ min are Gamma values of pure sandstone and pure mudstone, respectively; I γ is argillaceous content index; ρ is rock density, g/cm 3 ; and G cur is Hilchie index, which is related to geological age. Generally, 37 is taken for Tertiary system and 2 is taken for old stratum. The in-situ stress calculation methods are as shown in Formula (24) [40] (Huang 1985): where ω 1 and ω 2 are the tectonic stress coefficient in the direction of maximum and minimum horizontal stress. In the same block, ω 1 and ω 2 are constant; TVD is the vertical depth at measuring, m. According to the logging data, the Young's modulus, tensile strength, shale content, Poisson's ratio and the in situ stress coefficient are shown in Figure 8. AC and DTS represent P-wave moveout and S-wave moveout, respectively.
According to the logging data, the Young's modulus of the Da'anzhai section of Well NC2H ranges from 18.23 to 23.70 GPa, the tensile strength ranges from 6.27 to 8.18 MPa, the shale content ranges from 0.39 to 0.65, and the Poisson's ratio varies between 0.18 and 0.24. The maximum horizontal principal stress coefficient (σ H ) ranges from 2.51 to 3.32, the minimum horizontal principal stress coefficient (σ h ) from 2.20 to 2.41, and the vertical in situ stress coefficient (σ v ) is about 2.55. The measured parameters are brought into the fracability evaluation model, and the variation law of the FI of the Da'anzhai Section and the depth, time, and dimensionless distance from the wellbore axis are calculated, as shown in Figure 9. According to the logging data, the Young's modulus of the Da'anzhai section of Well NC2H ranges from 18.23 to 23.70 GPa, the tensile strength ranges from 6.27 to 8.18 MPa, the shale content ranges from 0.39 to 0.65, and the Poisson's ratio varies between 0.18 and 0.24. The maximum horizontal principal stress coefficient (σH) ranges from 2.51 to 3.32, the minimum horizontal principal stress coefficient (σh) from 2.20 to 2.41, and the vertical in situ stress coefficient (σv) is about 2.55. The measured parameters are brought into the fracability evaluation model, and the variation law of the FI of the Da'anzhai Section and the depth, time, and dimensionless distance from the wellbore axis are calculated, as shown in Figure 9. In Figure 9, the redder the color, the higher the FI and the greater the fracturing feasibility; the bluer the color, the lower the fracturing index and the fracturing feasibility are. Figure 9 shows that, in the Da'anzhai section of Well NC2H, the fracability index at the well wall near 2600 m is 0.78, which is the maximum value in this interval, thus indicating that it is easier to conduct fracturing operations in that region. Simultaneously, at the same depth, the FI near the wellbore is higher because the rock near the wellbore demonstrates a high degree of hydration and is more prone to fracture. The change law of FI with time and dimensionless distance is similar to the change law of water saturation rate in Figure 3. Figure 9 shows that when H = 2600 m, the rock is most likely to fracture. Therefore, according to Formulas (9) and (10), the variations in rock fracture pressure with  According to the logging data, the Young's modulus of the Da'anzhai section of Well NC2H ranges from 18.23 to 23.70 GPa, the tensile strength ranges from 6.27 to 8.18 MPa, the shale content ranges from 0.39 to 0.65, and the Poisson's ratio varies between 0.18 and 0.24. The maximum horizontal principal stress coefficient (σH) ranges from 2.51 to 3.32, the minimum horizontal principal stress coefficient (σh) from 2.20 to 2.41, and the vertical in situ stress coefficient (σv) is about 2.55. The measured parameters are brought into the fracability evaluation model, and the variation law of the FI of the Da'anzhai Section and the depth, time, and dimensionless distance from the wellbore axis are calculated, as shown in Figure 9. In Figure 9, the redder the color, the higher the FI and the greater the fracturing feasibility; the bluer the color, the lower the fracturing index and the fracturing feasibility are. Figure 9 shows that, in the Da'anzhai section of Well NC2H, the fracability index at the well wall near 2600 m is 0.78, which is the maximum value in this interval, thus indicating that it is easier to conduct fracturing operations in that region. Simultaneously, at the same depth, the FI near the wellbore is higher because the rock near the wellbore demonstrates a high degree of hydration and is more prone to fracture. The change law of FI with time and dimensionless distance is similar to the change law of water saturation rate in Figure 3. Figure 9 shows that when H = 2600 m, the rock is most likely to fracture. Therefore, according to Formulas (9) and (10) In Figure 9, the redder the color, the higher the FI and the greater the fracturing feasibility; the bluer the color, the lower the fracturing index and the fracturing feasibility are. Figure 9 shows that, in the Da'anzhai section of Well NC2H, the fracability index at the well wall near 2600 m is 0.78, which is the maximum value in this interval, thus indicating that it is easier to conduct fracturing operations in that region. Simultaneously, at the same depth, the FI near the wellbore is higher because the rock near the wellbore demonstrates a high degree of hydration and is more prone to fracture. The change law of FI with time and dimensionless distance is similar to the change law of water saturation rate in Figure 3. Figure 9 shows that when H = 2600 m, the rock is most likely to fracture. Therefore, according to Formulas (9) and (10), the variations in rock fracture pressure with time and dimensionless distance at a depth of 2600 m were calculated. The results are shown in Figure 10.
Processes 2022, 10, x FOR PEER REVIEW 14 of 16 time and dimensionless distance at a depth of 2600 m were calculated. The results are shown in Figure 10. According to the data in Figure 10, the fracture pressure near the wellbore is 66.80 MPa, and the original fracture pressure of the formation is 72.31 MPa. The fracturing pressure is the lowest near the borehole wall and gradually increases to the original fracturing pressure of the formation as the distance increases. Combined with the previous variation law of water saturation rate, it shows that after the shale around the well is hydrated, the tensile strength of the rock decreases, and the fracture pressure of the formation decreases. The change in the fracture pressure is similar to the variation law of the water saturation rate of the rock.

Conclusions
In this study, the rock mechanical properties of the shale oil reservoirs in the Da'anzhai Section of the Sichuan Basin were discussed, and the fracability index calculation model was established based on three factors: brittleness index, horizontal differential stress, and fracture toughness. A method for evaluating the fracability of shale oil reservoirs, considering hydration, was established. The main conclusions are the following: 1. Through the triaxial mechanical test, we found that the average compressive strength of the immersed rock samples decreased by 37.99%, the average elastic modulus decreased by 53.36%, and the average Poisson's ratio increased by 68.75%. These results show that the mechanical strength and elastic parameters of the shale rock in the Da'anzhai Section are obviously affected by the immersion in the drilling fluid, which further indicates that the change in the water saturation rate of the shale will also change its formation fracability. 2. A fracability evaluation model that comprehensively considers the brittleness index, in situ stress, and fracture toughness was established, and the law of water phase intrusion into shale oil reservoirs was analyzed. The changes in the Young's modulus, Poisson's ratio, fracture toughness, and FI with time and dimensionless distance under hydration were discussed. The results show that the change in the formation FI is closely related to the change in the water saturation rate. The rock at the wellbore is mostly affected by hydration, while the FI is the highest. 3. Taking Well NC2H as an example, the distribution law of the FI of the Da'anzhai section of the well with the well depth, dimensionless distance, and time was calculated and analyzed. The results show that the well is most prone to fracture near H = 2600 m. Second, the variation law of fracture initiation pressure before and after hydration was compared and analyzed, and it was found that the formation fracture initiation pressure decreased after hydration, thus making it easier to conduct fracturing operations. According to the data in Figure 10, the fracture pressure near the wellbore is 66.80 MPa, and the original fracture pressure of the formation is 72.31 MPa. The fracturing pressure is the lowest near the borehole wall and gradually increases to the original fracturing pressure of the formation as the distance increases. Combined with the previous variation law of water saturation rate, it shows that after the shale around the well is hydrated, the tensile strength of the rock decreases, and the fracture pressure of the formation decreases. The change in the fracture pressure is similar to the variation law of the water saturation rate of the rock.

Conclusions
In this study, the rock mechanical properties of the shale oil reservoirs in the Da'anzhai Section of the Sichuan Basin were discussed, and the fracability index calculation model was established based on three factors: brittleness index, horizontal differential stress, and fracture toughness. A method for evaluating the fracability of shale oil reservoirs, considering hydration, was established. The main conclusions are the following:

1.
Through the triaxial mechanical test, we found that the average compressive strength of the immersed rock samples decreased by 37.99%, the average elastic modulus decreased by 53.36%, and the average Poisson's ratio increased by 68.75%. These results show that the mechanical strength and elastic parameters of the shale rock in the Da'anzhai Section are obviously affected by the immersion in the drilling fluid, which further indicates that the change in the water saturation rate of the shale will also change its formation fracability.

2.
A fracability evaluation model that comprehensively considers the brittleness index, in situ stress, and fracture toughness was established, and the law of water phase intrusion into shale oil reservoirs was analyzed. The changes in the Young's modulus, Poisson's ratio, fracture toughness, and FI with time and dimensionless distance under hydration were discussed. The results show that the change in the formation FI is closely related to the change in the water saturation rate. The rock at the wellbore is mostly affected by hydration, while the FI is the highest.

3.
Taking Well NC2H as an example, the distribution law of the FI of the Da'anzhai section of the well with the well depth, dimensionless distance, and time was calculated and analyzed. The results show that the well is most prone to fracture near H = 2600 m. Second, the variation law of fracture initiation pressure before and after hydration was compared and analyzed, and it was found that the formation fracture initiation pressure decreased after hydration, thus making it easier to conduct fracturing operations.

4.
The evaluation method of fracability considering hydration established in this paper cannot only be used to calculate the hydration caused by drilling fluid invading formation during drilling but also the hydration caused by the change in water saturation rate caused by other factors.