Theoretical Analysis of the Mechanism of Fracture Network Propagation with Stimulated Reservoir Volume (SRV) Fracturing in Tight Oil Reservoirs

Stimulated reservoir volume (SRV) fracturing in tight oil reservoirs often induces complex fracture-network growth, which has a fundamentally different formation mechanism from traditional planar bi-winged fracturing. To reveal the mechanism of fracture network propagation, this paper employs a modified displacement discontinuity method (DDM), mechanical mechanism analysis and initiation and propagation criteria for the theoretical model of fracture network propagation and its derivation. A reasonable solution of the theoretical model for a tight oil reservoir is obtained and verified by a numerical discrete method. Through theoretical calculation and computer programming, the variation rules of formation stress fields, hydraulic fracture propagation patterns (FPP) and branch fracture propagation angles and pressures are analyzed. The results show that during the process of fracture propagation, the initial orientation of the principal stress deflects, and the stress fields at the fracture tips change dramatically in the region surrounding the fracture. Whether the ideal fracture network can be produced depends on the geological conditions and on the engineering treatments. This study has both theoretical significance and practical application value by contributing to a better understanding of fracture network propagation mechanisms in unconventional oil/gas reservoirs and to the improvement of the science and design efficiency of reservoir fracturing.


Introduction
Stimulated reservoir volume (SRV) fracturing has become an important technique to improve well production and the recovery of unconventional reservoirs in oil/gas field development.
With the development of the SRV fracturing technique, studying the mechanism of fracture network propagation with SRV fracturing has been an academic focus, including the variation rules of stress fields, mechanical mechanisms of fracture propagations and the corresponding geological and engineering factors.
Fracture propagation is strongly influenced by the mechanical interaction between neighboring fractures throughout the fracture growth history. This mechanical interaction, called the stress shadowing effect, is quantified using a boundary-element numerical model originally developed for natural-fracture propagation simulation. Erdogan and Sih [1] calculated the stress distribution near the crack tip under two special conditions, which involving symmetric and skew-symmetric case, by using stress field equations of the composite type I-II fracture at fracture tips in polar coordinates. In a 2D, plane-strain, displacement discontinuity solution, Crouch and Starfield [2] described the normal and shear stresses acting on one fracture element caused by the opening and shearing displacement discontinuities from all fracture elements. To consider the 3D effect due to finite fracture height, Olson [3] introduced a 3D correction factor to the influence coefficients. Since then, many scholars have done a lot of work in this regard, whose including Olson [4], Weng [5] and Kresse [6] and [7], et al. Nagel et al. [8] presented a discrete element model numerical study of multi-well completions simulated in a fully hydro-mechanical coupled fashion, and considered the influence of changes in the in-situ pressure in order to more completely understand the mechanical interactions between propagating hydraulic fractures and natural fractures during multi-well completions.
The methods of microseismic monitoring, numerical simulation and rock physics experiments have typically been used to analyze the mechanisms of fracture propagation patterns (FPP) with SRV fracturing in unconventional reservoirs. Complex hydraulic fracture geometry is regarded as a common place occurrence in naturally fractured reservoirs largely based on the interpretation of microseismic data by Fisher et al. [9], Maxwell et al. [10], Daniels et al. [11] and Le Calvez et al. [12]. In recent years, some progress has been made in this study by use of the numerical simulation method with flexible ways, which including Planar Fracture Model (PFM) and Unconventional Fracture Model (UFM) involving Complex Fracture Network Model (CFNM) and Discrete Fracture Network Model (DFNM). Zhang et al. [13] presented a PFM incorporating a full solution of coupled elasticity and fluid equations, but it is limited to the 2D plane strain conditions. Olson et al. [14] and [15] presented and modified the CFNM that is capable of predicting hydraulic fracture propagation and interaction with pre-existing natural fractures based on the Extended Finite Element Method (XFEM), and it was developed by Weng et al. [5] and Kresse et al. [7]. However, great uncertainty exists with respect to the nature of the natural fracture system and importantly how the hydraulic fracture interacts with it. To provide a framework to address these uncertainties, DFNM was presented by Rogers et al. [16] and was improved based on the discrete element models by Nagel et al. [17]. Laboratory experiments are another way to simulate the process of hydraulic fracture propagation in perexisting fracture formation. Blanton [18], Renshaw and Pollard [19] and Gu et al. [20] studied and expanded the approach of hydraulic fracture crossing or arrest in frictional interfaces within Devonian shale and hydro-stone. Several events may occur during the time period of the hydraulic fracture propagating towards the natural fracture and beyond. Zhou et al. [21] delineated natural fractures by inserting paper sheets of varying strength into cement blocks to simulate per-existing fractures of various "shear strength". Olson et al. [22] performed hydrostone block experiments to examine the effect of cemented natural fractures on hydraulic fracture propagation. The results indicate that oblique embedded fractures are more likely to divert a fluid-driven hydraulic fracture and hydraulic fracture and natural fracture interaction took three usual patterns.
Griffith [23] first postulated the maximum stress theory, which can predict the direction of crack growth. After its initial publication, the theory was re-examined by Erdogan and Sih [1], who proposed that the maximum circumferential stress in the fracture tip should be regarded as a controllable parameter of the composite crack extension by experimenting with the extension of a central crack under symmetrical plane loading in brittle materials (Plexiglass).
Blanton [18] and [24] used blocks of Devonian Shale cast in hydro-stone to show that a hydraulic fracture crosses a pre-existing fracture only under high differential stresses and high a ngles of approach. He noted that in most cases, fractures were either diverted or arrested by pre-existing fractures. Warpinski and Teufel [25] recognized the differences between in-situ stresses as the governing factor of fracture cross-over or fracture diversion in mine-back experiments and presented the corresponding fracture propagation criterion. The resulting effect of the discontinuities depends mainly on ancillary parameters (permeability of the joints, frictional properties, in-situ stresses, etc.). Based on this former research, a simple criterion for predicting whether a fracture will propagate across a frictional interface oriented perpendicular to the approaching fracture has been derived by Renshaw and Pollard [19]. However, this criterion has many restrictions on intersections. Hopkins et al. [26] analyzed the performance of hydraulic fracturing in fractured Antrim shale and concluded that complex fracture paths are primarily controlled by the geometry of the natural fractures. This conclusion contrasted with Blanton's theory, which limits natural fracture dilation to locations with low differential stresses. Consequently, a simplified multiple hydraulic fracture propagation model was presented by Olson [14] to demonstrate several important in situ and operational procedure parameters that influence the resultant fracture network geometry from horizontal wells and in naturally fractured formations. The geometry is strongly influenced by horizontal plane in situ stress contrasts, net pressure magnitude, hydraulic fracture height and the presence of natural fractures. Subsequently, Olson and Dahi-Taleghani [4] defined a parameter called the relative net pressure and developed a simplified numerical model to account for the mechanical interaction between pressurized fractures. The results show that fracture pattern complexity is strongly controlled by the magnitude of the hydraulic fracture net pressure relative to the in situ horizontal differential stress and the geometry of the natural fractures.
Previous studies have investigated theoretical models based on the crack mechanics theory or rock stress experiments, but all the above approaches may be considered to be less rigorous because the branch-fracture propagation description is not comprehensive. The branch-fracture propagation description assumes the stress field and propagation angle to be constant. In fact, the interference of the combined stress field would result in a continuous, mutative propagation angle at the fracture tip. This study presents a theoretical model of fracture network propagation with SRV fracturing in tight oil reservoirs. Based on the modified DDM, mechanical mechanism analysis and initiation and propagation criteria, a reasonable solution of the theoretical model is obtained and verified with the numerical discrete method. Using theoretical arithmetic and computer programming, the variation in the formation of the stress field, FPP and branch fracture propagation angle and pressure are analyzed.

Geological Conditions of Fracture Network Formation in Tight Oil Reservoirs
SRV fracturing is essential in the development of tight oil reservoirs, and the formation of fracture networks is key in increasing production. Several important geological factors affect the formation of fracture networks: rock mineral composition, rock mechanical properties, natural fractures and horizontal stress field.
1. Mineral Composition. Mineral composition is key to rock brittleness [27]. The brittleness of a mineral increases when it contains silicon or calcium, making the mineral more easily fractured and facilitating the activation of natural fractures. Therefore, an induced fracture network is easily formed after fracturing. However, high clay concentrations may create a softer and more ductile formation, complicating the formation of stable flow paths and complex fracture networks.
2. Rock Mechanical Properties. To a large extent, rock mechanical properties are affected by the mineral composition (or rock brittleness). The concept of rock brittleness combines both Poisson's ratio and Young's modulus. Poisson's ratio is used to reflect the ability of the rock to fail under stress; as values of Poisson's ratio decrease, the brittleness of the rock increases. Young's modulus reflects the ability of the rock to maintain a fracture [28]. A higher Young's modulus indicates more brittle rock.
3. Natural Fractures. The orientation of the natural fracture will affect the propagation of the hydraulic fracture and the formation of fracture network [29]. The results show that when the angle between the natural fracture and hydraulic fracture is less than 30 degrees, the probability that the natural fracture will be activated and the original propagation path will be changed increases, causing the easier formation of a fracture network. However, once the angle is larger than 60 degrees, the formation of a fracture network becomes impossible. In the case that the angle is between 30 and 60 degrees, the formation of a fracture network depends on the horizontal stress difference.

Horizontal Stress
Field. Under the condition of a low horizontal stress difference, the hydraulic fracture would propagate along natural fractures, and a radial fracture network would be induced. Hence, the stress anisotropy in fractured formations with low horizontal stress difference is weak, and as fractures propagate along different directions, the resulting net pressure difference is less [30] and [31]. Therefore, hydraulic fractures propagate more easily along natural fractures in random directions to form a fracture network.

Combined Stress Field Distribution
In view of the problems of stress concentration in the fracture tips, this section uses a modified displacement discontinuity method (DDM) by introducing the shape function as a Lagrange interpolation cardinal function to improve the precision of discontinuous displacement for the discrete state. The computational model of the combined stress field is constructed by the stress shadow effect and stress boundary condition. This model considers the influence of the fracture length on the stress field distribution during the fracturing treatments. 3.1.1 DDM Based on the Lagrange Interpolation. DDM is based on the analytical solution to the problem of constant discontinuity in displacement over a finite line segment in the x, y plane of an infinite elastic solid. The basic solutions proposed by DDM to the problem of conventional boundary and crack boundary are discontinuous displacements [2].
Considering a 2-D model, when the surfaces of several arbitrary fracture elements bear a stress field at a discrete state, slippage will occur between the upper and lower surfaces. This slippage is displacement discontinuity, and the relative magnitude of the slippage is a discontinuous displacement. We divide the single fracture into several equal elements, each of which has a length of 2a. SON is the local coordinate system. D s and D n are the tangential and normal discontinuous displacements, respectively (Fig 1). When |s| a, the equation of a discontinuous displacement can be defined as follows: D s ¼ u s ðs; 0 À Þ À u s ðs; 0 þ Þ D n ¼ u n ðn; 0 À Þ À u n ðn; 0 þ Þ ð1Þ ( where u s and u n are the tangential and normal displacement components, respectively; '+' indicates the upper surface and '-' indicates the lower surface of the fracture element. To improve the accuracy of the discontinuous displacement between one fracture element and the other two adjacent elements, the shape function N m (ξ)(m = 1 * 3) is introduced as the Lagrange-interpolation cardinal function. Consequently, the discontinuous displacement of arbitrary point can be rewritten as follows: Based on the theory of elastic mechanics, assuming that the fracture element surface produces symmetrical tangential and normal discontinuous displacements, the relational expression among the displacements (u s and u n ), stress fields (σ ss , σ nn and τ sn ) and discontinuous displacements (D s and D n ) of the fracture element in the local coordinate system are given by the following: where F n ðn ¼ 2 $ 5Þ is the derivative function of f m (s, n) that is only related to the fracture element coordinates; Àf m;nn ¼ @ 2 f m ðs;nÞ @s 2 , For the function F m ðI 0 ; I 1 ; , the analytical formula can be given as follows: I 0 ¼ nðy 1 À y 2 Þ À ðs À aÞln r 1 þ ðs þ aÞln r 2 À 2a where y 1 ¼ arctan n sÀa and y 2 ¼ arctan n sþa are the angles between the S-axis and two lines connecting any point (s, n) to (-a, 0) and (a, 0), respectively; r 1 ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi ðs À aÞ 2 þ n 2 q and r 2 ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi ðs þ aÞ 2 þ n 2 q are the distances of any point (s, n) to (-a, 0) and (a, 0), respectively.

Combined Stress Field and Stress Boundary Conditions.
In the local coordinate system SO'N, the stress field (normal stress and shear stress) of the fracture element i under an in-situ stress field can be defined as follows: where ðs i nn Þ 0 is the normal stress of fracture element i in the N direction; ðt i sn Þ 0 is the shear stress of fracture element i; ðs i xx Þ 0 and ðs i yy Þ 0 are the original normal stresses of fracture element i in the X and Y directions, respectively; and ðt i xy Þ 0 is the original shear stress of fracture element i.
For a fracture element, assuming the surface displays a symmetrical net pressure, the normal and shear stresses on the surface of the element are generated from the superposition of the stress shadows produced by other elements (as shown in Fig 2). Combining the equations of the additional stress field [2], the stress field (normal and shear stress) of fracture element i becomes the following: where C ij ss and C ij sn are the elastic coefficients of the normal stress due to the effects of shear and the open discontinuous displacements (D j s and D j n ) of element j on fracture element i in the X direction, respectively; C ij C ij ns and C ij nn are the elastic coefficients of the normal stress due to the effects of shear and the open discontinuous displacements (D j s and D j n ) of element j on fracture element i in the Y direction, respectively; cos 2b F 7 Þ; β is the angle between the N direction of element i and the X axis.
The above equations are the stress boundary conditions for the discrete fracture element. The equations on the left are the known stress boundary conditions, and the right side are the 2N variables containing the fracture element discontinuous displacements (D j s and D j n ) which can be solved through the establishment of 2N equations.
Both the propagation of hydraulic fractures and their extension upon encountering natural fractures are influenced by the combined stress field during the SRV fracturing of naturally fractured reservoirs. The combined stress field is formed by the superposition of two stress fields, which are the in-situ stress field caused by the long-term geological conditions and the additional stress field induced by multi-fracture propagation. The equations of the combined stress field are as follows: where σ xx and σ yy are the normal stresses in the X and Y directions, respectively, under the combined stress field; τ yy is the shear stress in the combined stress field; ðs i xx Þ 0 and ðs i yy Þ 0 are the additional normal stresses of fracture element i in the X and Y directions, respectively; ðt i xy Þ 0 is the additional shear stress of fracture element i; C ij sns and C ij snn are the elastic coefficients of the shear stress from the effects of the shear and open discontinuous displacements (D j s and D j n ) of fracture element j on fracture element I, respectively, C ij sns ¼ 2G½sin2b

Numerical Solution by
Programming. In the process of iteratively solving the pressure and geometric parameters for each fracture element, the additional stress field (normal and shear stress) produced by stress shadow effects must be recalculated at every time step. The stress field must be superimposed at the final time step to determine the new composite stress field distribution in the global coordinate system XOY.  Because natural fractures develop in tight oil reservoirs, SRV fracturing often induces complex fracture network growth. However, to form the complex network, the natural fractures must initiate and produce open rupture or shear slippage, which can effectively link natural fractures to form branch fractures. When a hydraulic fracture propagates along a natural fracture tip after intersecting with the natural fracture, branching of the hydraulic fracture to form a complex fracture network is possible. However, the fluid pressure at the intersection point must exceed the fluid pressure drop from the intersection point of the natural fracture tip and satisfy the initiation and propagation criteria at the natural fracture tip [30]. Based on the above analysis of the mechanical relationship between the fluid pressure at the intersection point of a hydraulic and a natural fracture and other geological factors (normal stress, shear stress and rock tensile strength), the hydraulic fracture and natural fracture interaction occurred in the following three patterns (as shown in Fig 4). FPP 1. The hydraulic fracture does not pass through the natural fracture but induces the diversion of fracture propagation. The propagation of the hydraulic fracture diverts when the fluid pressure of the intersection point cannot cause the hydraulic fracture to pass through the natural fracture. Rather, the hydraulic fracture induces natural fracture slippage, initiates at the tips, and diverts propagation along the direction of maximum principal stress. The mechanical condition is given by the following: where σ n and σ t are the normal stresses perpendicular and parallel, respectively, to the natural According to the theory of two-dimensional linear elastic fracture mechanics, the following equations can be obtained: s n ¼ 1 2 ðs xx þ s yy Þ À 1 2 ðs xx À s yy Þcos2b, s t ¼ 1 2 ðs xx þ s yy Þþ where k nf is the natural fracture permeability, expressed in mD; μ is the fluid viscosity in the reservoir, expressed in mPa.s; ϕ nf is the dimensionless natural fracture porosity; C t is the natural fracture comprehensive compression coefficient, expressed in 1/MPa; p 0 is the original fluid pressure in the reservoir, expressed in MPa; L nf is the natural fracture length, expressed in m; and p i = p net(i) + σ yy , where p net(i) is the net pressure at the intersection point, expressed in MPa. Thus, the formula for the net pressure at the intersection point becomes the following: FPP 2. he hydraulic fracture passes through natural fracture, inducing a partial opening of the natural fracture. The hydraulic fracture passes directly through the natural fracture and propagates along the initial direction (the maximum principal stress direction). According to the Blanton criteria [24], the fluid pressure at the intersection point must satisfy the following formula: Therefore, the following formula for the net pressure at the intersection point is obtained: FPP 3. The hydraulic fracture passes through the natural fracture, inducing the diversion of natural fracture propagation. Simultaneously, the net pressure at the intersection point cannot only exceed the rock tensile strength and normal stress to pass through the natural fracture but also must induce natural fracture slipping.

Fracture Initiation and Propagation Criteria
In this section, the mathematical model of the branch fracture initiation and propagation criteria is constructed based on the maximum circumferential stress theory and the stress intensity factor theory of the composite type I-II fracture.
3.3.1 Maximum Circumferential Stress Theory. Based on the maximum circumferential stress criterion, we can assume that the crack extends along the direction of the circumferential stress, reaching a maximum value; the crack initiates in this direction when the circumferential stress reaches a critical value.
(1) The stress field equations for a composite type I-II fracture at a fracture tip in polar coordinates (as shown in Fig 5) can be written as follows: 8 > > > > > > > > < > > > > > > > > : (2) The propagation angle equations of a composite type I-II fracture (for which the circumferential stress is equal to the limiting value) can be written as follows: The equation can be rewritten as follows: ffiffiffiffiffiffiffi 2pr p cos y 2 ½K I sin y þ K II ð3cos y À 1Þ ¼ 0 ð13Þ

Mechanism of Fracture Network Propagation with SRV Fracturing
Because cos y 2 À Á 6 ¼ 0, the corresponding simplified equation can be given as follows: where M ¼ K I K II . Considering that the value of θ 0 is constant when the shear stress changes, Eq (14) can be rewritten as cos . This equation shows that the propagation angle does not depend on the material properties because the parameters in this equation do not contain the material constants E and ν.
(3) Maximum Circumferential Stress In the function ðs yy Þ max ¼ 1 ffiffiffiffi ffi 2pr p ½K I g I yy ðyÞ À K II g II yy ðyÞ, g I yy ðyÞ ¼ cos 3 y 2 , g II yy ðyÞ ¼ À 3 2 cos y 2 sin y, and their curves are shown in Fig 6. For the composite type I-II fracture in Fig 6, when M is positive and negative, the conditions at the fracture tips achieve maximum circumferential stress values of -70.53° θ 0 and 0 θ 70.53°, respectively.
(4) Fracture Initiation Criteria where (σ θθ ) cr is the critical value of the maximum circumferential stress and depends on the fracture toughness of the type I fracture, which is also a constant related only to the material properties and independent of the composite fracture state. The fracture criterion must describe the type I fracture, and the material properties must be invariant. For the pure type I fracture, θ 0 equals zero. The value of K I equals K IC when initiating the fracture. Therefore, the critical stress equation of fracture initiation becomes the following: Consequently, the composite type I-II fracture initiation criteria becomes the following: where K eq and K IC are the equivalent fracture strength and the fracture toughness at the fracture tips, respectively.

Stress Intensity Factor (SIF).
Considering a vertical buried natural fracture in infinite elastic rock. We assume that the long axis of the elliptical cross section is 2a, the natural fracture is in a state of two-axial compression stress, and the surface produces a shear stress τ xy under the combined stress field (maximum and minimum principal stress: σ xx and σ yy ). We construct a local rectangular coordinate system, SON, by using the center of the fracture surface as the origin and the S axis as the long axis. In the global coordinate system, XOY, the angle between the normal direction of the crack surface and the ox axis is α (as shown in Fig 7).
According to the 2-D linear elastic theory of fracture mechanics, natural fractures are an independent system unit. Therefore, when hydraulic fractures interact with natural fractures in the direction perpendicular to the normal stress caused by the combined stress field, the pressure component of the fracturing fluid can be expressed as the following: where p net(i) is the net pressure inside the fracture at the intersection point between the natural fracture and hydraulic fracture, expressed in MPa. Therefore, under the multi-fracture combined stress field, the type I SIF at the branch fracture tips that is produced by the normal stress can be defined as follows: The shear stress parallel to the natural fracture surface under the stress field can be obtained as follows: Consequently, under a multi-fracture combined stress field, the type II SIF at the natural fracture tips produced by the shear stress can be defined as follows:

Results, Analysis and Discussion
We use the basic parameters in Table 1 to analyze the mechanism of fracture network propagation with SRV fracturing in tight oil reservoirs, including the variation rules of stress field formation, FPP, branch-fracture propagation angle and pressure.

Model Comparison and Validation
To compare the superiority of different model solution methods, a line fracture of length 6a subjected to a pressure p was analyzed. The analytical solution of the normal discontinuous displacement can be given by the following: We divide the line fracture into three equal elements, each of which has a length of 2a. Considering the solutions of analytical method, DDM and modified DDM, the calculated results  including the normal discontinuous displacement and normal stress of fracture in the N direction can be obtained ( Table 2). As shown in Table 2, for the normal discontinuous displacement and normal stress at elements adjacent to the fracture center and tips, the modified DDM's solution is in excellent agreement with the analytical solutions and more accurate than the DDM's solution. The relative errors of modified DDM's solutions are always less than 5% for all the elements. But the relative errors of DDM's solutions are 17.33% and 31.67% for the center element and both sides of the elements, respectively. Therefore, the simple DDM's solution cannot give reliable results for field points close to the fracture tip. The more elaborate DDM based on the Lagrange interpolation (modified DDM) will be necessary in this case.

Variation Rule of Ground Stress Field
According to the formula for the combined stress field, the stress distribution curve on the X and Y axis and the 2-D ground stress field distribution over different fracture lengths can be obtained, the results are shown in Fig 8 and Fig 9.  Fig 8A-8D and Fig 9A-9F compare the differences of the stress distribution curves on the X and Y axis and the stress field distribution over different fracture lengths during a signal fracture propagation.
1. The phenomenon of local stress concentration at fracture tips is induced by the stress near the fracture tips, which changes dramatically in the X-direction. The normal stress in the Xdirection increases at the fracture tips and declines rapidly far away from the fracture tip; the normal stress in the X-direction is always greater than the Y-direction.
2. The initial orientation of the maximum horizontal principal stress deflects in the vicinity of the fracture in the Y-direction. The normal stress in the X-direction is less than in the Y-direction for the deflection area and is consistent with the initial orientation out of the deflection area.
3. In the process of fracture propagation, the induced shear stress is small and forms the approximate shape the numeral "8" in the 2-D stress field distribution.
During fracture propagation, the stress influence range increases gradually, the phenomena of local stress concentration and stress deflection quickly become increasingly evident, and the stress field changes intensify. In this table, M ¼ À D n E pað1Àn 2 Þ ; RE is the relative error; D n is the normal discontinuous displacement; E is the young's modulus; p is the net pressure; v is the poisson's ratio; a is the half length of element; σ nn is the normal stress of fracture in the N direction. doi:10.1371/journal.pone.0125319.t002

Mechanical Analysis of the FPP
FPP 3 forms a fracture network system that combines the main fracture and branch fractures whereas FPP 1 can only form the fracture network that has a specific bandwidth. Therefore, the mechanical conditions of FPP 3 are the most favorable for forming a complex fracture network, followed by FPP 1; FPP 2 only forms a single fracture. Considering the mechanical formula for the FPP and the factors of the contact angle, horizontal stress difference and net pressure inside the fracture, three FPPs can be obtained under different mechanical control conditions, as shown in Fig 10. Fig 10A-10F compare the differences of FPP under different control conditions. The contact angle, stress difference and net pressure inside the fracture and other factors combined determine the FPP.
1. The stress difference has a greater effect on FPP. When the stress difference is 0 MPa, only the FPP 1 and FPP 3 conditions may arise; whereas when stress differences exist, all three FPPs may appear. For greater stress differences, greater mechanical values are required to produce FPP 3, which is not conducive to the formation of a fracture network.
2. The contact angle has some influence on FPP. A small contact angle with increases in pressure during fracture treatment will cause the intersection between the hydraulic fracture and natural fracture to first satisfy the mechanical conditions of FPP 1. Then, this condition will divert the natural fracture propagation to cause FPP 3 by passing through natural fractures, forming a more complex fracture network. A greater contact angle corresponds to a higher critical net pressure required to form the fracture network, i.e., resistance to the formation of the complex fracture network.
3. The rock tensile strength influences the hydraulic fractures passing through the natural fracture and directly affects the boundary between FPP 1 and FPP 3. Greater rock tensile strength corresponds to a greater value of the critical net pressure required to generate FPP 3, which is required for SRV fracture treatment.

Variation Rule of Branch Fracture Propagation Angle
In light of the formula for the compound fracture propagation angle, the variation curve of the branch fracture propagation angle with contact angle, the horizontal stress difference and net pressure inside the fracture can be calculated, and the result is shown in Fig 11. 1. The contact angle is an uncontrollable parameter. When the hydraulic fracture interacts with the natural fracture along the direction of maximum principal stress in a tight reservoir, the branch fracture propagation angle increases with the increasing contact angle. When the contact angle is 0°, the rock tensile failure will cause a pure type I (open type) fracture, in which the hydraulic fracture passes through the natural fracture without deflection. When the contact angle is 90°, shear and slip failure occur, causing pure type II (slip type) fractures; the deflection angle of the branch fracture then reaches a maximum value of 70.53°.
2. Greater net pressure favors the formation of a complex fracture network. When the stress difference is constant and the net pressure inside the fracture increases, the branch fracture propagation deflection angle decreases. Additionally, the radius of curvature of the deflecting propagation trajectory increases. Therefore, the branch fracture can encounter more natural fractures in the vertical direction, increasing its bandwidth. This process is conducive to the production of a more complex fracture network. During hydraulic fracturing, the net pressure is dynamic with fracture propagation and is a controllable factor; thus, the net pressure can be increased by increasing the flow rate of the fracturing fluid in a tight oil reservoir.
3. Smaller stress differences favor the formation of the complex fracture network. When the net pressure inside the fracture is constant and the stress difference increases, on the one hand, the branch fracture deflecting propagation angle increases and the radius of curvature of the deflecting propagation trajectory decreases. Additionally, the branch fracture converges towards the direction of principal stress over a small range because of the higher pressure drop within the fracture, allowing the branch fracture to continue to propagate more difficultly. On the other hand, when the stress difference are too low, the hydraulic fracture propagates more easily along the natural fracture surface that has a weak strength. This situation is not conducive to increasing of band length of the fracture network because it reduces the chance to encounter more natural fractures in the horizontal direction.
Therefore, the SRV fracture treatment requires not only the hydraulic fracture to propagate farther along the direction of maximum principal stress to encounter more natural fractures but also a smaller propagation angle to produce a larger fracture-network band width and to encounter more natural fractures. Combining both of these conditions, we can achieve the formation of a complex-fracture network system.

Variation Rule of Branch Fracture Propagation Pressure
Based on the fracture initiation and propagation criteria of mixed-mode I and II, the variation rule of branch fracture propagation pressure with contact angle, natural fracture half-length, horizontal stress difference and rock fracture toughness can be calculated, as seen in Fig 12.  Fig 12A-12F compares the differences in branch fracture propagation pressure under different control conditions, of which the following are several of the notable effects: 1. The contact angle has a large effect on the propagation pressure required to divert the propagation of the natural fracture. When the contact angle is 0°, pure type I open fractures are produced. These fractures are produced simply by overcoming the minimum principal stress and rock fracture toughness. The path of the hydraulic fracture coincides with the natural fracture and passes through it. When the contact angle increases further, the effective net pressure at the intersection of the hydraulic fracture and natural fracture decreases, and the required propagation pressure required, hindering the formation of a fracture network. However, when the contact angle approaches 90°, pure type II fractures occur, in which steering the natural fracture requires maximum net pressure. Therefore, producing a fracture network is difficult.
2. The natural fracture half-length has some influence on the propagation pressure required to induce a branch fracture. On the one hand, with increases in the natural fracture length, the generated equivalent-stress intensity factor increases, reducing the magnitude of the net pressure that can induce the extension of a branch fracture; this reduction in the magnitude of net pressure occurs when the natural fracture is shorter. On the other hand, longer natural fractures correspond to a greater loss of net pressure at the tip when the hydraulic fractures encounter the natural fractures, resulting in the need for higher pressure to form a fracture network.
3. Stress differences have a greater effect on the propagation pressure required for the deviation of the natural fracture. When the stress difference increases, the net pressured required to steer the natural fractures to extend along the branch increases markedly. It becomes more difficult to expand or initiate a branch; this difficulty occurs because of the increases in the branch propagation angle. The increased angle reduces the radius of curvature to a smaller trajectory that reaches the main branch in the direction of the maximum principle horizontal stress. This intersection decreases the net pressure rapidly, resulting in the need for a higher pressure to induce a branch fracture.
4. The rock fracture toughness has a smaller effect on the propagation pressure required for the deviation of the natural fracture. With increases in the rock fracture toughness, the natural fracture propagation resistance increases. Therefore, natural fractures require greater pressure to produce a net branch fracture, resulting in increased difficulty to form a fracture network.
The simulation results show that the natural fracture not only steering extends with the insitu differential stress and rock fracture toughness but is also influenced by the length and orientation of the natural fractures. Therefore, during the SRV fracturing in the tight oil reservoir, more mature and dense natural fractures are suggested to be selected in the block, and the influence of the natural fracture length and orientation of the fracture network morphology should be considered.

Influence of SRV on Oil Well Performance
Varied lithologic compaction, poor physical property and high flow resistance are inherent characteristic for tight reservoir. Hydraulic fracturing technology are widely used to change the flowing regime near the wellbore, and increase the channel of oil flowing to improve the well performance. To reveal the influence of SRV on oil well performance, considering the situation of without fracturing (SRV=0), conventional fracturing (Bi-wing & Symmetric fracture) and SRV fracturing, the production performance including daily production and cumulative production curve of oil well can be obtained by numerical reservoir simulator, as shown in Fig 13. Results show that without fracturing, oil production is extremely low in tight oil reservoir and it's hard to achieve economic rates, so hydraulic fracturing is a necessary process to accelerate production. Conventional fracturing could create bi-wing & symmetric fracture, providing significantly more fracture-surface contact with the reservoir, thus accerates oil production obvisously. However, SRV fracturing produces higher rates with the following mechanism.
SRV fracturing could communicate hydraulic fracture with natural fracture and create 3D fracture network with certain length, width and height, that will increase the stimulated reservoir volume and enhance the overall system permeability. In tight oil reservoir, diffusivity-related pore-pressure changes cannot move very far from the fracture, so little fluid is produced outside the extent of the SRV. Within the SRV, SRV maximize the contact area between fracture-surface and matrix, so oil could flow from the matrix to the created fracture network through the shortest path easily. SRV could maximize producing ratio and recovery percent of reservoir, having significant effect on oilfield production.
So it's necessary to study the mechanism of fracture network propagation with SRV fracturing in tight oil reservoirs. The relevant results studied including geological conditions and engineering treatments have both theoretical significance and practical application value by contributing to a better understanding of fracture network propagation mechanisms in tight oil reservoirs and to the improvement of the science and design efficiency of reservoir fracturing.
Conclusions SRV fracturing has emerged as an important technology to exploit oil/gas from unconventional reservoirs, in which fluids are injected at high rate and pressure to crack the reservoir rock in order to create pathways with higher permeability to accelerate production. This paper analyzes the mechanism of fracture network propagation with SRV fracturing in tight oil reservoirs from both geological and engineering perspectives, the following conclusions can be drawn: 1. During the process of fracture propagation, the stress field produces local stress concentration at the fracture tips, and the initial orientation of maximum horizontal principal stress deflects into the area around the fracture.
2. Propagation pattern 3 (in which a hydraulic fracture passes through a natural fracture, thereby inducing the diversion of the natural fracture propagation) is the most conducive to forming a network fracture system by combining a principle fracture and a branch fracture.
3. A specific set of reservoir geological attributes (higher brittle mineral contents and stronger elastic characteristic of the rock (Poisson's ratio and Young's modulus), lower tensile strength and fracture toughness of rock, mature natural fractures (natural fracture density), and lower stress anisotropy (horizontal stress difference)) are more conducive to the propagation of hydraulic fractures into a complex fracture network. According to the engineering conditions of fracturing operations, higher net treatment pressure (propagation pressure), reasonable fracturing azimuth direction (contact angle) and larger fracturing scale (hydraulic fracture length) will assist the formation of complex fracture networks with SRV fracturing in tight oil reservoirs.

Author Contributions
Conceived and designed the experiments: YS LR. Performed the experiments: YS CX. Analyzed the data: YS LR WW. Contributed reagents/materials/analysis tools: LR FM WW. Wrote the paper: LR FM CX.