Inﬂuence of Water–Oil Saturation on the Fracture Process Zone: A Modiﬁed Dugdale–Barenblatt Model

: The wetting and nonwetting ﬂuid saturations in porous reservoirs always change during long-term injection and production. The fracture process zone (FPZ) is a prominent feature in the rock fracture process. If the FPZ properties are inﬂuenced by pore ﬂuids, the process of hydraulic fracturing will change greatly. The existing models do not consider the role of pore ﬂuid when characterizing the FPZ. In this paper, a modiﬁed Dugdale–Barenblatt (D–B) model with capillary pressure is proposed. The model reﬂects the fact that the FPZ length decreases nonlinearly with the increase in capillary pressure, and it reveals the mechanism of capillary pressure on the equivalent fracture cohesion in the FPZ, which affects the FPZ length. Three-point bending tests were carried out on sandstone under various ﬂuid saturations through digital image correlation (DIC), acoustic emission (AE), and scanning electron microscope (SEM). It was found that the FPZ length of the water–oil-saturated samples was 30–50% smaller than that of water-saturated/oil-saturated samples due to the capillary pressure effect, and the modiﬁed D–B model was well consistent with the experiments. The AE behaviors of different saturated samples were not the same: The cumulative AE signals changed abruptly at 90% of the peak load for the water–oil-saturated samples and at 50% of the peak load for water-saturated samples. This demonstrated that the effect of capillary pressure was more obvious than the weakening effect of microstructural damages. The signiﬁcant inﬂuence of capillary pressure on FPZ requires continuous recognition in hydraulic fracturing design.


Introduction
As global energy provided by conventional reservoirs has become more difficult, the development of oil and gas in unconventional reservoirs is increasingly necessary [1]. The wetting and nonwetting fluid saturations in reservoirs always change as a result of long-term injection and production [2]. The fractures produced by hydraulic fracturing and the existing fractures in depleted reservoirs are affected by the water-oil saturation. During hydraulic fracturing, fracture initiation and propagation are closely related to the nonlinear elastic response that occurs in the fracture process zone (FPZ) [3]. Therefore, investigating the effect of water-oil saturation on the FPZ is of great importance for production improvement.
The effect of pore fluids on mechanical properties needs to be taken into consideration during the rock fracture process. Many scholars have performed a series of experiments to research the influence of water content on the mechanical properties of rocks. The results showed that the static Young's modulus [3,4], dynamic Young's modulus [5,6], compressive stress [3,7], tensile stress [4], and fracture toughness [8][9][10] decreased drastically with the increase in water content. The effect of pore fluids on the mechanical properties is a complicated process [5], which includes physical-chemical and mechanical effects (such as capillary pressure). Researchers [11][12][13][14][15][16][17][18][19] have focused on the physical-chemical effects of water on rock strength, such as clay swelling and chemical reactions. Besides, capillary pressure (the pressure difference of wetting and nonwetting fluids) has been shown to increase strength [5] and fracture toughness [9]. When the content of the wetting and nonwetting fluids are different, there are differences in the effects of the physical-mechanical and capillary pressure, which results in different trends in the mechanical properties. The fracture process zone is a prominent feature of the rock fracture process. As the mechanical properties of rocks can be influenced by pore fluids, it is not difficult to speculate that the fracture process zone of rock is closely related to pore fluid saturation.
In addition to studying the effects of pore fluids on mechanical parameters and the fracture process, researchers have proposed various models to describe the FPZ. In 1960, Dugdale [20] found that a plastic zone at the fracture tip was concentrated at an angle of 45 • to the plane of the mild sheet, and the ribbon yielding model was proposed. In 1962, Barenblatt [21] found that a small nonlinear stress cohesion zone at the fracture tip eliminated stress singularity. Researchers have since combined these two models into the classic Dugdale-Barenblatt (D-B) model, and many modified D-B models have been proposed. Classical FPZ models include the virtual fracture model [22], the passivation model [23], and the equivalent model [24]. All these models assume that strain softening occurs at the fracture tip. The FPZ length can be calculated based on these models. Labuz et al. [25] presented an experimental analysis of fracture growth in dry charcoal granite that incorporated the concept of the process zone, and they proposed a model to calculate the FPZ length. Although many models have been provided, little research has concentrated on the influence of capillary pressure on the FPZ length.
The primary objective of this study was to investigate the influence of capillary pressure on the FPZ length. For this, the authors propose a modified classical D-B model considering capillary pressure. Three-point bending experiments were conducted with digital image correlation (DIC) and acoustic emission (AE) to verify the proposed model. The microstructural variations of different saturated samples were compared with scanning electron microscope (SEM) results.

Modified Fracture Process Zone Model
As capillary pressure can increase strength and fracture toughness [5,9] and as the magnitude of capillary pressure is similar to tensile strength [26], the capillary pressure in the fracture tip is a non-neglectable factor during the fracture process. It can be inferred that capillary pressure will have an effect on the FPZ length. However, the influence of capillary pressure on the FPZ length is not considered in classical FPZ models. In this paper, the authors therefore propose a modified Dugdale-Barenblatt model that considers capillary pressure.

Liquid Bridges between Idealized Grains
When the pores of porous rock are filled with wetting and nonwetting fluids, the liquid bridge between the grains is shown in Figure 1. The liquid bridge force between grains comprises two parts: surface tension from the wetting phase and the capillary pressure formed by the two fluids. The direction of the liquid bridge force points to the center of the two grains.
The liquid bridge force is expressed as in Equation (1) [27]: where d (µm) is the half-length between the two grains centers; R 1 (µm) is the grain radius; ϕ (degree) is the half-filling angle; F cap (N) is the liquid force caused by the capillary pressure; γ (mN/m) is the surface tension of the two phases; F surf (N) is the liquid bridge force caused by the surface tension of the wetting phase fluid; and F (N) is the resultant liquid bridge force.

Modified FPZ Model
The geometry of a fracture presented in Figure 1 consists of a prefabricated fracture a0 (mm), a real fracture a1 (mm), and a microfracture a2 (mm) as an FPZ [28].
For the modified FPZ model, the authors assumed that the dimension of the fracture process zone was small in comparison with the length of the whole fracture. The grain size of the samples was the same. There was no stress singularity at the tip of the fracture.
The superposition of the stress intensity factors caused by the external load, yield stress, and liquid bridge force at the fracture tip was equal to 0, as shown in Equation (2).
where (Pa·m 0.5 ), (Pa·m 0.5 ), (Pa·m 0.5 ), (mm), (mm) are the stress intensity factors produced by the external load, the yield stress, the liquid bridge force at the fracture tip, the prefabricated fracture length, and the prefabricated fracture length, respectively. The stress intensity factor produced by the liquid bridge force at the fracture tip can be shown as in Equation (3). The FPZ length with pore fluids can be simplified into Equation (5), and the stress produced by the liquid bridge can be expressed as in Equation (6).
where (mm) is the FPZ length, and (MPa) is the stress produced by the liquid bridge. As shown in Equation (5), the formulation to determine the FPZ length will change when the capillary pressure exists.

Estimating the Influence of Capillary Pressure on an Idealized Sample
In order to quantify the influence of capillary pressure on FPZ length, the authors assumed that the grain radius , the contact angle , the surface tension of the water and oil, and the peak load were 25 μm, 13° [29], 36 mN/m [30], and 0.8 kN, respectively. For simplification, the relationship between and d was assumed to be linear. It was assumed that when the largest filling angle was 60°, the corresponding length between the grains was 20 μm; when the smallest filling angle was 30°,

Modified FPZ Model
The geometry of a fracture presented in Figure 1 consists of a prefabricated fracture a 0 (mm), a real fracture a 1 (mm), and a microfracture a 2 (mm) as an FPZ [28].
For the modified FPZ model, the authors assumed that the dimension of the fracture process zone was small in comparison with the length of the whole fracture. The grain size of the samples was the same. There was no stress singularity at the tip of the fracture.
The superposition of the stress intensity factors caused by the external load, yield stress, and liquid bridge force at the fracture tip was equal to 0, as shown in Equation (2).
where K σ (Pa·m 0.5 ), K y (Pa·m 0.5 ), K F (Pa·m 0.5 ), a (mm), c (mm) are the stress intensity factors produced by the external load, the yield stress, the liquid bridge force at the fracture tip, the prefabricated fracture length, and the prefabricated fracture length, respectively. The stress intensity factor produced by the liquid bridge force at the fracture tip can be shown as in Equation (3). The FPZ length with pore fluids can be simplified into Equation (5), and the stress produced by the liquid bridge can be expressed as in Equation (6).
where l p (mm) is the FPZ length, and σ F (MPa) is the stress produced by the liquid bridge. As shown in Equation (5), the formulation to determine the FPZ length will change when the capillary pressure exists.

Estimating the Influence of Capillary Pressure on an Idealized Sample
In order to quantify the influence of capillary pressure on FPZ length, the authors assumed that the grain radius R 1 , the contact angle θ, the surface tension γ of the water and oil, and the peak load were 25 µm, 13 • [29], 36 mN/m [30], and 0.8 kN, respectively. For simplification, the relationship between ϕ and d was assumed to be linear. It was assumed that when the largest filling angle was 60 • , the corresponding length between the grains was 20 µm; when the smallest filling angle was 30 • , the corresponding length between the grains was 0.1 µm. This linear relationship can be expressed as in Equation (7).
where d (µm) is the half-length between the two grains centers, and ϕ (degree) is the half-filling angle.
As the half-filling angle is not constant in reality, two extreme cases for simplification were assumed, i.e., that all the filling angles were either 30 • or 60 • . Based on this assumption, the authors obtained the corresponding FPZ length, and the true length was between the maximum and the minimum. According to Equation (5), the decrease in FPZ length could be up to 30% due to the influence of capillary pressure because the FPZ length decreased by 30% and 5% compared with the dry sample for filling angles of 30 • and 60 • , respectively.

Sample Preparation
Due to the relative uniformity in diameter of the particles and the homogeneity of the mechanical properties, sandstone was selected as the samples, with an average grain diameter in the range of 25-100 µm, uniaxial compressive strength of 70 MPa, tensile strength of 7 MPa, elastic modulus of 8 GPa, porosity of 15%, and Poisson's ratio of 0.2. The X-ray diffraction (XRD) analysis found the average mineral contents, which are listed in Table 1. I/S, S, and K represent illite/smectite, smectite, and kaolinite, respectively. Table 1. Characterization of the minerals in the samples. Following the International Society for Rock Mechanics and Rock Engineering (ISRM) test standard [31], the samples were cut into semicircular plates with prefabricated fractures in the center. The dimensions of the samples are shown in Table 2, where a (mm) represents the length of the prefabricated fracture, R (mm) represents the radius of the samples, and B (mm) represents the thickness of samples. After processing the samples, the fluid saturating treatments [19] were conducted. Firstly, the dry sample was placed in a sealed container, which was then evacuated. Secondly, the saturated fluid was injected into the container under the negative pressure produced by the vacuum. Finally, the sample was saturated in the fluids for two months so that the oil/water molecules moved into the small pores. The saturating fluids were water and kerosene. For sandstone, water was the wetting fluid, while kerosene was the nonwetting fluid. For DIC, artificial speckles needed to be made on the surface of the sample. After the saturating treatments, one surface was sealed with paraffin and the other was painted with a thin layer of matte white paint to create a white base film, which was facing the camera. Then, a matte black paint was applied to make speckles. For high-quality specks, matte white and black paint should be sprayed alternately.

Experimental Equipment and Methods
The experimental devices are illustrated in Figure 2, where DIC and AE were used to record the fracture process. Internal damages to the rock during the fracture process are due to the change of the microstructure [32][33][34] and are detected by AE [35]. The DIC technology can be used to obtain the deformation of the grains on the surface [36][37][38][39].
To calculate the FPZ length from the displacement [36], DIC technology was deployed during the fracture process. Speckles were located and recorded under loading through a charge-coupled device (CCD) high-speed camera with an effective pixel of 2062 × 2062, a 35 mm fixed focus lens, and a high-speed capacity at a rate of two frames/s. The width and height of the camera window view focused on the 50 mm on both sides of the prefabricated fractures and 40 mm at the top of the prefabricated fracture, respectively. A stabilized white light source was used for a better light effect. As shown in Figure 2a, there were four sensors arranged on the light source side and the back side. The AE signal acquisition device was a SAEU2S-centralized AE system, and the threshold was 40 dB.
The load was applied with a hydraulic servo control system at the rate of 0.1 mm/min. Simultaneous recording of the AE signals (AE energy, amplitude, hit count, and the arrival time of signals) and high-speed camera were carried out until the whole fracture process was completed. For DIC, artificial speckles needed to be made on the surface of the sample. After the saturating treatments, one surface was sealed with paraffin and the other was painted with a thin layer of matte white paint to create a white base film, which was facing the camera. Then, a matte black paint was applied to make speckles. For high-quality specks, matte white and black paint should be sprayed alternately.

Experimental Equipment and Methods
The experimental devices are illustrated in Figure 2, where DIC and AE were used to record the fracture process. Internal damages to the rock during the fracture process are due to the change of the microstructure [32][33][34] and are detected by AE [35]. The DIC technology can be used to obtain the deformation of the grains on the surface [36][37][38][39].
To calculate the FPZ length from the displacement [36], DIC technology was deployed during the fracture process. Speckles were located and recorded under loading through a charge-coupled device (CCD) high-speed camera with an effective pixel of 2062 2062, a 35 mm fixed focus lens, and a high-speed capacity at a rate of two frames/s. The width and height of the camera window view focused on the 50 mm on both sides of the prefabricated fractures and 40 mm at the top of the prefabricated fracture, respectively. A stabilized white light source was used for a better light effect. As shown in Figure 2a, there were four sensors arranged on the light source side and the back side. The AE signal acquisition device was a SAEU2S-centralized AE system, and the threshold was 40 dB.
The load was applied with a hydraulic servo control system at the rate of 0.1 mm/min. Simultaneous recording of the AE signals (AE energy, amplitude, hit count, and the arrival time of signals) and high-speed camera were carried out until the whole fracture process was completed.

Determination of FPZ Length
Under an external load, the initiation and propagation of a microfracture will lead to an increase in the fracture tip displacement. The fracture extends unstably only when the crack tip opening displacement (CTOD) reaches the critical value. For tight rocks, the fracture propagates unstably at the peak load, and the zone in front of prefabricated fracture can be defined as the FPZ. Figure 3 shows that the incremental horizontal displacement contours merge at one point, which is defined as the FPZ tip [36]. At the peak load, the FPZ length is defined to be the distance between the merged point and the tip of the prefabricated fracture. The exact definition of the merged point is closely correlated with the size of the FPZ. There are different material properties inside and outside the FPZ. The deformation characteristics in those regions are different, and the position where the slope of incremental displacement abruptly changes is defined to be the merged point [36]. As shown in Figure 4, the horizontal displacements at the peak load of y = 0, y = 2.47 mm, y = 4.93 mm, and y = 5.32 mm were selected for the dry sample. It was not difficult to find that the displacement profiles of y = 0, y = 2.47 mm, y = 4.93 mm had the same characteristics; the slope could clearly be divided into three sections. However, the slope of the incremental horizontal displacement profile at y = 5.32 mm was nearly the same, which indicated that the slope abruptly changed at y = 4.93 mm. The corresponding point where the slope changed was defined as the merged point.

Determination of FPZ Length
Under an external load, the initiation and propagation of a microfracture will lead to an increase in the fracture tip displacement. The fracture extends unstably only when the crack tip opening displacement (CTOD) reaches the critical value. For tight rocks, the fracture propagates unstably at the peak load, and the zone in front of prefabricated fracture can be defined as the FPZ. Figure 3 shows that the incremental horizontal displacement contours merge at one point, which is defined as the FPZ tip [36]. At the peak load, the FPZ length is defined to be the distance between the merged point and the tip of the prefabricated fracture. The exact definition of the merged point is closely correlated with the size of the FPZ. There are different material properties inside and outside the FPZ. The deformation characteristics in those regions are different, and the position where the slope of incremental displacement abruptly changes is defined to be the merged point [36]. As shown in Figure 4, the horizontal displacements at the peak load of y = 0, y = 2.47 mm, y = 4.93 mm, and y = 5.32 mm were selected for the dry sample. It was not difficult to find that the displacement profiles of y = 0, y = 2.47 mm, y = 4.93 mm had the same characteristics; the slope could clearly be divided into three sections. However, the slope of the incremental horizontal displacement profile at y = 5.32 mm was nearly the same, which indicated that the slope abruptly changed at y = 4.93 mm. The corresponding point where the slope changed was defined as the merged point.

Results and Discussion
In the following section, the influence of capillary pressure on the FPZ length is discussed based on the measurement of the FPZ length according to DIC and the combined damage data from the AE.

Determination of FPZ Length
Under an external load, the initiation and propagation of a microfracture will lead to an increase in the fracture tip displacement. The fracture extends unstably only when the crack tip opening displacement (CTOD) reaches the critical value. For tight rocks, the fracture propagates unstably at the peak load, and the zone in front of prefabricated fracture can be defined as the FPZ. Figure 3 shows that the incremental horizontal displacement contours merge at one point, which is defined as the FPZ tip [36]. At the peak load, the FPZ length is defined to be the distance between the merged point and the tip of the prefabricated fracture. The exact definition of the merged point is closely correlated with the size of the FPZ. There are different material properties inside and outside the FPZ. The deformation characteristics in those regions are different, and the position where the slope of incremental displacement abruptly changes is defined to be the merged point [36]. As shown in Figure 4, the horizontal displacements at the peak load of y = 0, y = 2.47 mm, y = 4.93 mm, and y = 5.32 mm were selected for the dry sample. It was not difficult to find that the displacement profiles of y = 0, y = 2.47 mm, y = 4.93 mm had the same characteristics; the slope could clearly be divided into three sections. However, the slope of the incremental horizontal displacement profile at y = 5.32 mm was nearly the same, which indicated that the slope abruptly changed at y = 4.93 mm. The corresponding point where the slope changed was defined as the merged point.

Results and Discussion
In the following section, the influence of capillary pressure on the FPZ length is discussed based on the measurement of the FPZ length according to DIC and the combined damage data from the AE.

Results and Discussion
In the following section, the influence of capillary pressure on the FPZ length is discussed based on the measurement of the FPZ length according to DIC and the combined damage data from the AE.

Capillary Pressure Decreases FPZ Length
As stated in the previous section, FPZ length can be represented as the distance-along with the prefabricated fracture-between the initial position and the merged point, whose coordinates are specified with DIC at peak loads. For the dry, water-saturated, oil-saturated, and the water-oil-saturated samples, Figures 3 and 5-7 show the contours of the incremental horizontal displacements, respectively. Based on the method described in Section 3.3, the FPZ lengths for the samples with various saturations were obtained, as seen in Figure 8.
Capillary pressure can reduce the FPZ length. The average FPZ length of the water-oil-saturated samples was 3.32 mm-the shortest of all the samples-and was about 15-25 times the grain radius. The lengths of the FPZ for the oil-saturated, dry, and water-saturated samples were about 33%, 63%, and 110% larger than that of the water-oil-saturated sample, respectively. The capillary pressure only existed in the water-oil-saturated samples, thus it can be concluded that capillary pressure decreased the FPZ length.
Water-oil saturation can affect the pore structure of the sample and the cohesive strength of the grains. The changes in the microstructure of different saturated samples were observed by SEM in order to describe the influence of different types of pore fluids; the results are described in the following sections. Combined with the analysis of the forces between the grains, the influence of capillary pressure on FPZ length is then discussed.

Capillary Pressure Decreases FPZ Length
As stated in the previous section, FPZ length can be represented as the distance-along with the prefabricated fracture-between the initial position and the merged point, whose coordinates are specified with DIC at peak loads. For the dry, water-saturated, oil-saturated, and the water-oilsaturated samples, Figures 3 and 5-7 show the contours of the incremental horizontal displacements, respectively. Based on the method described in Section 3.3, the FPZ lengths for the samples with various saturations were obtained, as seen in Figure 8.
Capillary pressure can reduce the FPZ length. The average FPZ length of the water-oil-saturated samples was 3.32 mm-the shortest of all the samples-and was about 15-25 times the grain radius. The lengths of the FPZ for the oil-saturated, dry, and water-saturated samples were about 33%, 63%, and 110% larger than that of the water-oil-saturated sample, respectively. The capillary pressure only existed in the water-oil-saturated samples, thus it can be concluded that capillary pressure decreased the FPZ length.
Water-oil saturation can affect the pore structure of the sample and the cohesive strength of the grains. The changes in the microstructure of different saturated samples were observed by SEM in order to describe the influence of different types of pore fluids; the results are described in the following sections. Combined with the analysis of the forces between the grains, the influence of capillary pressure on FPZ length is then discussed.

Capillary Pressure Decreases FPZ Length
As stated in the previous section, FPZ length can be represented as the distance-along with the prefabricated fracture-between the initial position and the merged point, whose coordinates are specified with DIC at peak loads. For the dry, water-saturated, oil-saturated, and the water-oilsaturated samples, Figures 3 and 5-7 show the contours of the incremental horizontal displacements, respectively. Based on the method described in Section 3.3, the FPZ lengths for the samples with various saturations were obtained, as seen in Figure 8.
Capillary pressure can reduce the FPZ length. The average FPZ length of the water-oil-saturated samples was 3.32 mm-the shortest of all the samples-and was about 15-25 times the grain radius. The lengths of the FPZ for the oil-saturated, dry, and water-saturated samples were about 33%, 63%, and 110% larger than that of the water-oil-saturated sample, respectively. The capillary pressure only existed in the water-oil-saturated samples, thus it can be concluded that capillary pressure decreased the FPZ length.
Water-oil saturation can affect the pore structure of the sample and the cohesive strength of the grains. The changes in the microstructure of different saturated samples were observed by SEM in order to describe the influence of different types of pore fluids; the results are described in the following sections. Combined with the analysis of the forces between the grains, the influence of capillary pressure on FPZ length is then discussed.

Capillary Pressure Affects FPZ More than Microstructural Damages
The FPZ length of the water-saturated samples was the longest due to the damaging effect of the water-rock reaction on the microstructure. Some swelling, dissolution, and cementitious mineral loss always occur as a result of the existence of pore fluids [10], which weakens the ability to sustain the load. Figure 9a-c show the microstructures of the dry sample, and there were some flaky and floccular minerals among the grains. The composition of the constituent elements of these minerals, as found by the SEM analysis, are shown in Figure 10. According to Figure 10, the significant presence of Al was an indication of the presence of aluminosilicate, which can react with water as shown in Equations (8) and (9) [40]. The reaction led to mineral loss among the grains, which weakened the cementation; the microstructure of the sample then became the type shown in Figure 9d, and there were many microfractures among the grains. Figure 9e,f show the microstructure of the oil-saturated and water-oil-saturated samples, respectively. In Figure 9, it can be seen that the water-saturated sample had the most microfractures, and the water-oil-saturated sample was slightly affected; however, the microstructure of the oil-saturated sample was almost unaffected. This showed that the change in the microstructure was one of the reasons for the longer length of the FPZ in the watersaturated sample.
K0.9Al2.9Si3.1O10(OH)2 + nH2O → K0.9Al2.9Si3.1O10(OH)2nH2O (9) The attraction between the grains at the fracture tip due to the liquid bridge force from the capillary pressure was stronger than the weakening effect of the microstructural damages. As sandstone is a porous medium, pore fluids can produce a liquid bridge, which produces a liquid

Capillary Pressure Affects FPZ More than Microstructural Damages
The FPZ length of the water-saturated samples was the longest due to the damaging effect of the water-rock reaction on the microstructure. Some swelling, dissolution, and cementitious mineral loss always occur as a result of the existence of pore fluids [10], which weakens the ability to sustain the load. Figure 9a-c show the microstructures of the dry sample, and there were some flaky and floccular minerals among the grains. The composition of the constituent elements of these minerals, as found by the SEM analysis, are shown in Figure 10. According to Figure 10, the significant presence of Al was an indication of the presence of aluminosilicate, which can react with water as shown in Equations (8) and (9) [40]. The reaction led to mineral loss among the grains, which weakened the cementation; the microstructure of the sample then became the type shown in Figure 9d, and there were many microfractures among the grains. Figure 9e,f show the microstructure of the oil-saturated and water-oil-saturated samples, respectively. In Figure 9, it can be seen that the water-saturated sample had the most microfractures, and the water-oil-saturated sample was slightly affected; however, the microstructure of the oil-saturated sample was almost unaffected. This showed that the change in the microstructure was one of the reasons for the longer length of the FPZ in the watersaturated sample.
K0.9Al2.9Si3.1O10(OH)2 + nH2O → K0.9Al2.9Si3.1O10(OH)2nH2O (9) The attraction between the grains at the fracture tip due to the liquid bridge force from the capillary pressure was stronger than the weakening effect of the microstructural damages. As sandstone is a porous medium, pore fluids can produce a liquid bridge, which produces a liquid

Capillary Pressure Affects FPZ More than Microstructural Damages
The FPZ length of the water-saturated samples was the longest due to the damaging effect of the water-rock reaction on the microstructure. Some swelling, dissolution, and cementitious mineral loss always occur as a result of the existence of pore fluids [10], which weakens the ability to sustain the load. Figure 9a-c show the microstructures of the dry sample, and there were some flaky and floccular minerals among the grains. The composition of the constituent elements of these minerals, as found by the SEM analysis, are shown in Figure 10. According to Figure 10, the significant presence of Al was an indication of the presence of aluminosilicate, which can react with water as shown in Equations (8) and (9) [40]. The reaction led to mineral loss among the grains, which weakened the cementation; the microstructure of the sample then became the type shown in Figure 9d, and there were many microfractures among the grains. Figure 9e,f show the microstructure of the oil-saturated and water-oil-saturated samples, respectively. In Figure 9, it can be seen that the water-saturated sample had the most microfractures, and the water-oil-saturated sample was slightly affected; however, the microstructure of the oil-saturated sample was almost unaffected. This showed that the change in the microstructure was one of the reasons for the longer length of the FPZ in the water-saturated sample.
The attraction between the grains at the fracture tip due to the liquid bridge force from the capillary pressure was stronger than the weakening effect of the microstructural damages. As sandstone is a porous medium, pore fluids can produce a liquid bridge, which produces a liquid bridge force that resists particle separation. The liquid bridge force (especially capillary pressure) can strengthen the cohesion between particles. Although the damage to the water-oil-saturated samples was more serious than that of the oil-saturated samples, the length of the FPZ of the water-oil-saturated sample was smaller than that of the oil-saturated sample shown in Figure 8, which indicated that hydration reaction played a less important role than capillary pressure.
Energies 2018, 11, x FOR PEER REVIEW 9 of 14 bridge force that resists particle separation. The liquid bridge force (especially capillary pressure) can strengthen the cohesion between particles. Although the damage to the water-oil-saturated samples was more serious than that of the oil-saturated samples, the length of the FPZ of the water-oilsaturated sample was smaller than that of the oil-saturated sample shown in Figure 8, which indicated that hydration reaction played a less important role than capillary pressure.

Acoustic Emission Characteristics Affected by Water-Oil Saturation
Different saturated samples showed different AE behavior, which could be described from two aspects: the cumulative AE hit count and the AE cumulative energy [41,42].
The damage caused by water was the greatest of all the saturation scenarios. The cumulative AE hit count describes the damage magnitude of the sample with different water-oil saturations. If the damage is more serious, the significant cumulative hit count occurs earlier. The blue line in Figures  11 and 12 show typical curves for the cumulative AE hit count for the water-saturated sample (W-2) and the water-oil-saturated sample (W-O-1). The cumulative AE signals changed abruptly at 90% of the peak load for the water-oil-saturated samples; this change occurred at 50% of the peak load for the water-saturated samples.
The cohesion of the water-saturated samples was the weakest of all the samples. The AE cumulative energy can represent the strength of the sample with different water/oil contents. When the cohesion is weak, the energy needed is low. The green line in Figures 11 and 12 shows the cumulative energy over time for the water-saturated sample (W-2) and the water-oil-saturated

Acoustic Emission Characteristics Affected by Water-Oil Saturation
Different saturated samples showed different AE behavior, which could be described from two aspects: the cumulative AE hit count and the AE cumulative energy [41,42].
The damage caused by water was the greatest of all the saturation scenarios. The cumulative AE hit count describes the damage magnitude of the sample with different water-oil saturations. If the damage is more serious, the significant cumulative hit count occurs earlier. The blue line in Figures 11  and 12 show typical curves for the cumulative AE hit count for the water-saturated sample (W-2) and the water-oil-saturated sample (W-O-1). The cumulative AE signals changed abruptly at 90% of the peak load for the water-oil-saturated samples; this change occurred at 50% of the peak load for the water-saturated samples.
The cohesion of the water-saturated samples was the weakest of all the samples. The AE cumulative energy can represent the strength of the sample with different water/oil contents. When the cohesion is weak, the energy needed is low. The green line in Figures 11 and 12 shows the cumulative energy over time for the water-saturated sample (W-2) and the water-oil-saturated sample (W-O-1). It was evident that the energy released from the water-oil-saturated sample was about 5 times that of the water-saturated sample. It was shown that the cohesion between the grains in the water-saturated sample was weaker than that of the water-oil-saturated samples. When the cohesion was weaker, the FPZ length was larger, which was well consistent with the modified model. sample (W-O-1). It was evident that the energy released from the water-oil-saturated sample was about 5 times that of the water-saturated sample. It was shown that the cohesion between the grains in the water-saturated sample was weaker than that of the water-oil-saturated samples. When the cohesion was weaker, the FPZ length was larger, which was well consistent with the modified model.

Impact on Engineering Design
Drawing on the above discussion, it is not difficult to find that capillary pressure reduced the FPZ length. This may be of benefit to modify the design of hydraulic fracturing.
During the process of designing hydraulic fractures, engineers always set the volume of the injected liquid to control the length of the hydraulic fracture. However, capillary pressure is always neglected. Due to capillary pressure, the fracture toughness is higher than that obtained in the lab, and the actual length of the fracture may be shorter than the designed one. The design length of the fracture should be longer than the designed one. This may be one of the reasons for the low production.  sample (W-O-1). It was evident that the energy released from the water-oil-saturated sample was about 5 times that of the water-saturated sample. It was shown that the cohesion between the grains in the water-saturated sample was weaker than that of the water-oil-saturated samples. When the cohesion was weaker, the FPZ length was larger, which was well consistent with the modified model.

Impact on Engineering Design
Drawing on the above discussion, it is not difficult to find that capillary pressure reduced the FPZ length. This may be of benefit to modify the design of hydraulic fracturing.
During the process of designing hydraulic fractures, engineers always set the volume of the injected liquid to control the length of the hydraulic fracture. However, capillary pressure is always neglected. Due to capillary pressure, the fracture toughness is higher than that obtained in the lab, and the actual length of the fracture may be shorter than the designed one. The design length of the fracture should be longer than the designed one. This may be one of the reasons for the low production.

Impact on Engineering Design
Drawing on the above discussion, it is not difficult to find that capillary pressure reduced the FPZ length. This may be of benefit to modify the design of hydraulic fracturing.
During the process of designing hydraulic fractures, engineers always set the volume of the injected liquid to control the length of the hydraulic fracture. However, capillary pressure is always neglected. Due to capillary pressure, the fracture toughness is higher than that obtained in the lab, and the actual length of the fracture may be shorter than the designed one. The design length of the fracture should be longer than the designed one. This may be one of the reasons for the low production.

Conclusions
In this paper, the influence of water-oil saturation on the FPZ was studied by means of capillary pressure and the weakening effect of microstructural damages. More attention was paid to the influence of capillary pressure. The FPZ is a prominent feature in the rock fracture process. However, existing models do not consider the role of pore fluid when characterizing the FPZ. In this study, a modified Dugdale-Barenblatt (D-B) model with capillary pressure is therefore proposed. The model reflects the fact that the FPZ length decreases nonlinearly with the increase in capillary pressure, and it reveals the mechanism of capillary pressure on the equivalent fracture cohesion in the FPZ, which affects the FPZ length. Three-point bending tests were carried out on sandstone under various fluid saturations. It was found that the FPZ length of the water-oil-saturated samples was 30-50% smaller than that of the water-saturated/oil-saturated samples due to the capillary pressure effect; the modified D-B model was well consistent with the experiments. The AE behaviors of the different saturated samples were not the same: The cumulative AE signals changed abruptly at 90% of the peak load for water-oil-saturated samples; this change occurred at 50% of the peak load for the water-saturated samples. This demonstrated that the effect of capillary pressure was more obvious than the weakening effect of the microstructural damage. This research will be beneficial for understanding the fracture process in reservoirs, where the saturations of wetting and nonwetting fluids always change.