Pressure Transient Performance for a Horizontal Well Intercepted by Multiple Reorientation Fractures in a Tight Reservoir

: A fractured horizontal well is an e ﬀ ective technology to obtain hydrocarbons from tight reservoirs. In this study, a new semi-analytical model for a horizontal well intercepted by multiple ﬁnite-conductivity reorientation fractures was developed in an anisotropic rectangular tight reservoir. Firstly, to establish the ﬂow equation of the reorientation fracture, all reorientation fractures were discretized by combining the nodal analysis technique and the fracture-wing method. Secondly, through coupling the reservoir solution and reorientation fracture solution, a semi-analytical solution for multiple reorientation fractures along a horizontal well was derived in the Laplace domain, and its accuracy was also veriﬁed. Thirdly, typical ﬂow regimes were identiﬁed on the transient-pressure curves. Finally, dimensionless pressure and pressure derivative curves were obtained to analyze the e ﬀ ect of key parameters on the ﬂow behavior, including fracture angle, permeability anisotropy, fracture conductivity, fracture spacing, fracture number, and fracture conﬁguration. Results show that, for an anisotropic rectangular tight reservoir, horizontal wells should be deployed parallel to the direction of principal permeability and fracture reorientation should be controlled to extend along the direction of minimum permeability. Meanwhile, the optimal fracture number should be considered for economic production and the fracture spacing should be optimized to reduce the ﬂow interferences between reorientation fractures.


Introduction
As one of the most effective stimulation techniques to increase the recovery of hydrocarbons, fractured horizontal wells have been widely adopted to develop tight reservoirs. In some tight reservoirs, reorientation fractures may be formed when an increase in pore pressure in a hydraulically-fractured area or the depletion of the reservoir pressure occurs [1,2]. For reorientation fractures, previous studies mainly focused on field case studies and fracture propagation [2,3]. Compared with hydraulically-fractured vertical wells, fractured horizontal wells can significantly enhance well productivity by increasing the reservoir contact area. For reservoir engineers, pressure transient analysis is one of the most effective techniques to diagnose the flow behavior of hydraulically-fractured (1) The tight reservoir is rectangular in shape. Single-phase Darcy flow occurs in the anisotropic tight reservoir. The permeability of the reservoir in the xand y-direction is k x and k y , respectively. (2) The horizontal well is fully intercepted by arbitrary reorientation fractures with constant height and width. The horizontal wellbore is deployed parallel to the x-axis, and all fractures' height is equal to the reservoir thickness. (3) Fluid in the horizontal wellbore and reservoir is constant viscosity and slightly compressible, and flow rates from each reorientation fracture contribute to the horizontal well total rate, although they may change over time. The horizontal wellbore is infinite-conductivity and no pressure loss occurs along the wellbore. (4) This reservoir is fully penetrated by all reorientation fractures. For the i-th reorientation fracture, its principal fracture angle is θ i,1 , and its reoriented fractures angles are θ i,2 and θ i,3 , respectively. All reorientation fractures have a finite conductivity and their tips are assumed to be impermeable boundaries. (5) Gravity effect is negligible, as well as the influence of the temperature on different reservoir parameters.
Energies 2019, 12, x FOR PEER REVIEW 3 of 23 in Section 4. In Section 5, the flow characteristics and production rate distribution along each reorientation fracture will be discussed in detail. In addition, the influence of key parameters on pressure transient will also be investigated in Section 5, such as fracture number, fracture spacing and fracture configuration.

Physical Model
In this section, the physical model is introduced for pressure transient analysis of a horizontal well with multiple reorientation fractures in an anisotropic rectangular tight reservoir. Figure 1 presents the schematic of a horizontal well with two reorientation fractures, and every hydraulic fracture reoriented twice, which is significantly different from the physical model of Wu et al. [22]. Three parts, including principle fracture, irregular curve fracture and reoriented fracture, constitute the reorientation fracture. The basic assumptions are made as follows: (1) The tight reservoir is rectangular in shape. Single-phase Darcy flow occurs in the anisotropic tight reservoir. The permeability of the reservoir in the x-and y-direction is kx and ky, respectively. (2) The horizontal well is fully intercepted by arbitrary reorientation fractures with constant height and width. The horizontal wellbore is deployed parallel to the x-axis, and all fractures' height is equal to the reservoir thickness. (3) Fluid in the horizontal wellbore and reservoir is constant viscosity and slightly compressible, and flow rates from each reorientation fracture contribute to the horizontal well total rate, although they may change over time. The horizontal wellbore is infinite-conductivity and no pressure loss occurs along the wellbore. (4) This reservoir is fully penetrated by all reorientation fractures. For the i-th reorientation fracture, its principal fracture angle is θi,1, and its reoriented fractures angles are θi,2 and θi, 3, respectively. All reorientation fractures have a finite conductivity and their tips are assumed to be impermeable boundaries. (5) Gravity effect is negligible, as well as the influence of the temperature on different reservoir parameters.

Dimensionless Definitions
To simplify the governing equation and definite conditions, the following dimensionless parameters are defined for the mathematical model.
The dimensionless pressure in the tight reservoir and all reorientation fractures, and the dimensionless time, are expressed as follows: In the fracture model, the dimensionless fracture conductivity, F cD , the dimensionless flow rate, q fD , and the flow rate at the center of each fracture segment, q fwD , are defined as follows: Other dimensionless definitions in the reservoir model or the reorientation fracture model are given as follows: where β is the anisotropic factor and L R is the reference length and k is the equivalent system permeability, m i is the total discrete number of all irregular curve fractures and L i is the length of the equivalent i-th planar fracture for the irregular curve fracture.

Reservoir Flow Model
On the basis of the Laplace transform, a fracture segment solution can be obtained for an anisotropic rectangular tight reservoir as follows [22]: where (x Dmi , y Dmi ) is a movable infinitesimal unit along the fracture segment and (x D , y D ) is any point in the fracture segment. However, Equation (6) has an implicit assumption that the minute fracture segment should be parallel to the coordinate axis. In order to get the pressure solution of a minute fracture segment with a reorientation angle θ (Figure 2) between the fracture and the horizontal wellbore, the following equations are used: x Dmi = x Di + χ cos θ i y Dmi = y Di + χ sin θ i Energies 2019, 12, 4232 5 of 23 where (x Di , y Di ) is the midpoint of the fracture segment, θ i represents the angle of i-th fracture segment, and l Di represents the dimensionless length of i-th reorientation fracture segment.
permeability, mi is the total discrete number of all irregular curve fractures and Li is the length of the equivalent i-th planar fracture for the irregular curve fracture.

Reservoir Flow Model
On the basis of the Laplace transform, a fracture segment solution can be obtained for an anisotropic rectangular tight reservoir as follows [22]:  Dmi  D  D i  D  eD  eD  l   n  eD  D  Dmi  Dmi  D  Di  D  n  eD  n  n eD  eD  eD  l   n  eD   y  y  y  p  q  d l  x  y   y  y  y  n x  n x  q  d l  x  y where (xDmi, yDmi) is a movable infinitesimal unit along the fracture segment and (xD, yD) is any point in the fracture segment. As stated previously, the reorientation fracture cannot be described by the ideal planar fracture. In order to stimulate the flow characteristics of a horizontal well intercepted by multiple reorientation fractures, we used the approximate approach similar to the boundary element method to handle the irregular shaped boundaries [23]. As illustrated in Figure 3, the irregular curve fracture is replaced by a series of planar fracture segments. Thus, the reorientation fracture system is considered to be composed of several planar fracture segments. The number of planar fracture segments used to describe the irregular curve fracture largely depends on the seepage of the irregular curve fracture. However, Equation (6) has an implicit assumption that the minute fracture segment should be parallel to the coordinate axis. In order to get the pressure solution of a minute fracture segment with a reorientation angle θ (Figure 2) between the fracture and the horizontal wellbore, the following equations are used: where (xDi, yDi) is the midpoint of the fracture segment, θi represents the angle of i-th fracture segment, and △lDi represents the dimensionless length of i-th reorientation fracture segment. As stated previously, the reorientation fracture cannot be described by the ideal planar fracture. In order to stimulate the flow characteristics of a horizontal well intercepted by multiple reorientation fractures, we used the approximate approach similar to the boundary element method to handle the irregular shaped boundaries [23]. As illustrated in Figure 3, the irregular curve fracture is replaced by a series of planar fracture segments. Thus, the reorientation fracture system is considered to be composed of several planar fracture segments. The number of planar fracture segments used to describe the irregular curve fracture largely depends on the seepage of the irregular curve fracture.

Reorientation Fracture Flow Model
In this work, we assumed that the flow in the reorientation fracture is one-dimensional. Taking the i-th fracture segment as an example, we assume that the fluid flows from the (i+1)-th fracture segment to the i-th fracture segment ( Figure 4). Based on the results of Zhou et al. [19], the pressure difference between the i-th fracture segment and (i+1)-th fracture segment can be written as:

Reorientation Fracture Flow Model
In this work, we assumed that the flow in the reorientation fracture is one-dimensional. Taking the i-th fracture segment as an example, we assume that the fluid flows from the (i+1)-th fracture Energies 2019, 12, 4232 6 of 23 segment to the i-th fracture segment ( Figure 4). Based on the results of Zhou et al. [19], the pressure difference between the i-th fracture segment and (i+1)-th fracture segment can be written as: In order to obtain the wellbore pressure, we take the integration method from the first fracture segment near the horizontal wellbore to the i-th fracture segment and derive Equation (10): In addition, according to the principle of mass balance, we can obtain the following equations for each fracture segment: The inflow and outflow of each fracture node should satisfy the mass balance: where NI is total fracture segments of a wing of the i-th fracture.
Substituting Equation (11) and Equation (12) into Equation (10) yields: 1 Based on above dimensionless definitions, Equation (13) can be rewritten as follows: ( ) To the best of our knowledge, Equation (14) is a new fracture equation and can be employed to describe the flow inside the reorientation fracture.
In order to obtain the fracture solution in Laplace domain, we should apply the Laplace transform to Equation (14) based on tD and thus we can have: ( )  In order to obtain the wellbore pressure, we take the integration method from the first fracture segment near the horizontal wellbore to the i-th fracture segment and derive Equation (10): In addition, according to the principle of mass balance, we can obtain the following equations for each fracture segment: The inflow and outflow of each fracture node should satisfy the mass balance: where N I is total fracture segments of a wing of the i-th fracture. Substituting Equation (11) and Equation (12) into Equation (10) yields: Based on above dimensionless definitions, Equation (13) can be rewritten as follows: Energies 2019, 12, 4232 7 of 23 To the best of our knowledge, Equation (14) is a new fracture equation and can be employed to describe the flow inside the reorientation fracture.
In order to obtain the fracture solution in Laplace domain, we should apply the Laplace transform to Equation (14) based on t D and thus we can have: For the uniform dimensionless length of l Di for each fracture segment, Equation (15) will have the same form as that reported by Cinco-Ley et al., [24]. The simplification of Equation (15) is presented in Appendix A.
The flow equation of the reorientation fracture (Equation (15)) can be given in the matrix form: whereq fD represents the vector of the flow-rate in the reorientation fracture,p wD is the vector of the horizontal wellbore pressure andp fD is the vector of the pressure in the reorientation fracture, C, D, and E are the corresponding coefficient matrices obtained from Equation (15).

Semi-Analytical Solutions
In this paper, we get the semi-analytical pressure solution by discretizing each reorientation fracture into a large number of planar fracture segments, similar to the method illustrated in Figure 4. On the basis of the superposition principle, we can easily get the dimensionless pressure solution for the i-th fracture segment in the Laplace domain, which is presented as a discretized form: where N is the reorientation fracture number and: Considering the continuity conditions of pressure and flux along the fracture surface, the following equations can be obtained in the Laplace domain: Equations (17), (19), and (20) can be further written in the matrix form as follows: where A and B are the corresponding coefficient matrix obtained from Equation (17). Combining Equations (16) and (21) reveals the following simplified matrix: where G and H are the corresponding coefficient matrices obtained from Equation (15) and Equation (17). The sum of the flow rates of the reorientation fracture segments contributes to the total flow rate in the horizontal wellbore: Combining Equations (22) and (23), the following equation in the matrix form is obtained: where M T is the vector of dimensionless length for each reorientation fracture segment.
Equation (24) is the new semi-analytical solution for multiple reorientation fractures along a horizontal well in a tight reservoir in Laplace domain. According to Equation (24), the reference pressure of all fracture segments except the pressure of the horizontal wellbore does not need to be calculated and, thus, the calculation speed can be significantly increased, which is different from the work of Wu et al. [22]. By solving Equation (25) with the Gaussian elimination method, the pressure of the horizontal wellbore and the flow rate of each reorientation fracture segment can be obtained in Laplace space, and then can be further inverted into the real space by using the Stehfest numerical inverse method [25].

Verification and Comparison
The major advantage of the novel semi-analytical model is that we take the fracture reorientation and permeability anisotropy into consideration. To the best of our knowledge, there are no reports on a horizontal well perforated by a single reorientation fracture or multiple reorientation fractures in anisotropic rectangular reservoirs. Therefore, we conduct two special cases presented by Chen and Raghavan to verify the accuracy of our semi-analytical solution [8]. They calculated the pressure transient of a multiple fractured horizontal well in an isotropic rectangular reservoir using the reformulation of the point-source solution presented by Ozkan and Raghavan [4]. The comparison results are shown in Figure 5, where the semi-analytical solution shows a good fit to the results of Chen and Raghavan at every production time [8].

Transient Flow Characteristics
The transient flow behavior of this model is illustrated by transient-pressure type curves, which can be employed for pressure and rate transient analysis. In addition, some reservoir parameters can be obtained by type curves matching, such as the average reservoir permeability, fracture number, and fracture reorientation.
For a given set of parameters, the dimensionless pressure and its derivative of two reorientation fractures along the horizontal well in an anisotropic tight reservoir are displayed in Figure 6 to analyze its typical flow regimes. Figure  The major advantage of the novel semi-analytical model is that we take the fracture reorientation and permeability anisotropy into consideration. To the best of our knowledge, there are no reports on a horizontal well perforated by a single reorientation fracture or multiple reorientation fractures in anisotropic rectangular reservoirs. Therefore, we conduct two special cases presented by Chen and Raghavan to verify the accuracy of our semi-analytical solution [8]. They calculated the pressure transient of a multiple fractured horizontal well in an isotropic rectangular reservoir using the reformulation of the point-source solution presented by Ozkan and Raghavan [4]. The comparison results are shown in Figure 5, where the semi-analytical solution shows a good fit to the results of Chen and Raghavan at every production time [8].

Transient Flow Characteristics
The transient flow behavior of this model is illustrated by transient-pressure type curves, which can be employed for pressure and rate transient analysis. In addition, some reservoir Dimensionless pressure derivative,dp wD Figure 5. Comparison of pressure derivative for a horizontal well with three vertical fractures with the results of previous literature [8].
For a given set of parameters, the dimensionless pressure and its derivative of two reorientation fractures along the horizontal well in an anisotropic tight reservoir are displayed in Figure 6 to analyze its typical flow regimes. Figure 7 illustrates the fluid flow in the reservoir and reorientation fractures with different flow regimes. Without considering the transition flow regimes, we can easily observe six typical flow regimes in Figure 6: bilinear flow regime (BF), first linear flow regime (FLF), first radial flow regime (FRF), second radial flow regime (SRF), pseudo-radial flow regime (PRF), and pseudo-steady-state flow regime (PSSF).
Bilinear flow: In this regime, two types of linear flow, including the linear flow from the reservoir to reorientation fractures and from reorientation fractures to the horizontal wellbore, occur at the same time (Figure 7a). This flow behavior makes the slope of the dimensionless pressure derivative curve constant. In contrast to conventional multi-stage fractured horizontal wells, the streamline from the reservoir to the reorientation fracture may be not perpendicular to the fracture due to permeability anisotropy and fracture reorientation existing simultaneously. Consequently, the slope of the straight line at this regime is characterized by 1/4-1/3 in the log-log plot, not 1/4 as presented by Cinco-Ley et al. [24]. In addition, this flow behavior only occurs when the reorientation fracture conductivity is considered.
First linear flow: For the conventional finite-conductivity fractured horizontal well model, the slope of the pressure derivative curve is constant to 1/2 in this regime. Due to the fracture reorientation and permeability anisotropy, the slope in this regime is also larger than that of the conventional finite-conductivity fractured horizontal wells. During this regime, the fluid from the reservoir will linearly flow into the individual reorientation fracture and each reorientation fracture works independently on other reorientation fractures.
First radial flow: In the first radial flow regime, one of the most important transient flow characteristics is that the pressure derivative is expressed as a fixed value of 1/(2N). Consequently, the slope of the pressure derivative curve for various fractured horizontal well model may be different, which mainly depends on the fracture number. In this regime, the fluid around reorientation fractures will individually flow into the reorientation fractures. The duration of this regime is highly related to adjacent fracture spacing and reorientation fracture length. Generally, a large reorientation fracture length or a wide fracture spacing can definitely lead to a long duration of the first radial flow regime.   and reservoir. After this regime, the transient behavior in the reservoir surrounded by the outer fractures reaches a pseudo-radial flow regime. Pseudo radial flow: In this regime the dimensionless pressure derivative has a constant value of 0.5, and the derivative curve is parallel to the time axis ( Figure 6). In the meantime, the fluid flow among the reservoir, the reorientation fractures, and the horizontal wellbore reache a dynamic balance and, thus, the dimensionless pressure derivative stabilizes at 0.5. Pseudo-steady-state flow: A pseudo-steady-state flow regime occurs in the late-time period for all closed reservoirs. The typical characteristic of this regime is that the imensionless pressure and corresponding derivative curves increase rapidly and finally normalize to a straight line with a constant slope of 1. Generally, the seepage area in this flow regime always presents a circle ( Figure  7f).

Production Rate Distribution
Figures 8 and 9 present the production rate distribution in each reorientation fracture for a horizontal well with three reorientation fractures and the data used in Figures 8 and 9 are presented in Table. 1. In this section, we assumed that the inner reorientation fracture is fixed in the center of Bilinear flow: In this regime, two types of linear flow, including the linear flow from the reservoir to reorientation fractures and from reorientation fractures to the horizontal wellbore, occur at the same time (Figure 7a). This flow behavior makes the slope of the dimensionless pressure derivative curve constant. In contrast to conventional multi-stage fractured horizontal wells, the streamline from the reservoir to the reorientation fracture may be not perpendicular to the fracture due to permeability anisotropy and fracture reorientation existing simultaneously. Consequently, the slope of the straight line at this regime is characterized by 1/4-1/3 in the log-log plot, not 1/4 as presented by Cinco-Ley et al. [24]. In addition, this flow behavior only occurs when the reorientation fracture conductivity is considered.
First linear flow: For the conventional finite-conductivity fractured horizontal well model, the slope of the pressure derivative curve is constant to 1/2 in this regime. Due to the fracture reorientation and permeability anisotropy, the slope in this regime is also larger than that of the conventional finite-conductivity fractured horizontal wells. During this regime, the fluid from the reservoir will linearly flow into the individual reorientation fracture and each reorientation fracture works independently on other reorientation fractures.
First radial flow: In the first radial flow regime, one of the most important transient flow characteristics is that the pressure derivative is expressed as a fixed value of 1/(2N). Consequently, the slope of the pressure derivative curve for various fractured horizontal well model may be different, which mainly depends on the fracture number. In this regime, the fluid around reorientation fractures will individually flow into the reorientation fractures. The duration of this regime is highly related to Energies 2019, 12, 4232 11 of 23 adjacent fracture spacing and reorientation fracture length. Generally, a large reorientation fracture length or a wide fracture spacing can definitely lead to a long duration of the first radial flow regime.
Second radial flow: After the first radial flow, a second radial flow occurs in the tight reservoir. As demonstrated in Figure 6, the dimensionless pressure and its derivative curves are parallel to each other and the slope of these curves is a constant value of 0.36 in this regime. During this period, the interference of reorientation fractures starts to affect the flow between reorientation fractures and reservoir. After this regime, the transient behavior in the reservoir surrounded by the outer fractures reaches a pseudo-radial flow regime.
Pseudo radial flow: In this regime the dimensionless pressure derivative has a constant value of 0.5, and the derivative curve is parallel to the time axis ( Figure 6). In the meantime, the fluid flow among the reservoir, the reorientation fractures, and the horizontal wellbore reache a dynamic balance and, thus, the dimensionless pressure derivative stabilizes at 0.5.
Pseudo-steady-state flow: A pseudo-steady-state flow regime occurs in the late-time period for all closed reservoirs. The typical characteristic of this regime is that the imensionless pressure and corresponding derivative curves increase rapidly and finally normalize to a straight line with a constant slope of 1. Generally, the seepage area in this flow regime always presents a circle (Figure 7f).

Production Rate Distribution
Figures 8 and 9 present the production rate distribution in each reorientation fracture for a horizontal well with three reorientation fractures and the data used in Figures 8 and 9 are presented in Table 1. In this section, we assumed that the inner reorientation fracture is fixed in the center of the tight reservoirs and the distance between the inner reorientation fracture and the outer reorientation fracture is assumed to be constant at FS = 20. Meanwhile, the horizontal wellbore extends parallel to the x-axis. reorientation fracture is assumed to be constant at FS = 20. Meanwhile, the horizontal wellbore extends parallel to the x-axis. The production rate distribution of each reorientation fracture along the horizontal well in an isotropic tight reservoir is illustrated in Figure 8 when kx = ky. In the early flow regime (the black line in Figure 8), the production rate of each fracture segment is symmetrically distributed with respect to the wellbore node. As the pressure wave expands, the production rate in the reoriented fracture significantly increases along the extension of the reoriented fracture (the green line in Figure 8). For a long production time (Figure 8, tD = 10 4 ), the production rate distribution changes significantly. Since the outermost fractures (RF1 and RF3) have a larger drainage area, they have a higher production rate. In addition, owing to the flow interference, the inner fracture (RF2) has a lower production rate. Furthermore, the production rate distribution of the outermost fractures (RF1 and RF3) (the red line in Figure 8) is asymmetrical. The reoriented fractures away from the horizontal wellbore have a larger production rate due to the larger seepage area caused by fracture reorientation.     Figure 9 presents the production rate distribution for each reorientation fracture in a horizontal well in an anisotropic tight reservoir (kx = 4ky, FN = 3). The production rate distribution within each reorientation fracture is similar to that in Figure 8. However, for an anisotropic reservoir, fracture reorientation has a more significant effect on the production rate in each reorientation fracture. In the reorientation section, the production rate distribution curves vary sharply (Figure 9), which suggests that the effect of fracture reorientation in anisotropic tight reservoir on type curves is significant and, thus, cannot be neglected.

Parameter Influence on Transient Pressure Behavior
In this section, we analyze the sensitivity of the transient pressure of multiple reorientation fractures along a horizontal well in an anisotropic tight reservoir. We consider the case of a horizontal well with two reorientation fractures to analyze the effects of some key parameters, such as principal fracture angle (PFA), reoriented fracture angle (RFA), permeability anisotropic factor (PAF), and adjacent reorientation fractures spacing (FS), on the dimensionless pressure and its derivative curves. We also analyze the effect of reorientation fracture number and complex fracture configuration on type curves. Figure 10 presents the effect of PFA on type curves for a horizontal well with two reorientation fractures. The dimensionless fracture spacing is set at 20. As displayed in Figure 10, PFA has a weak impact on type curves in the early-time period, i.e., the bilinear flow regime, first linear flow regime, and first radial flow regime. For an anisotropic tight reservoir, the dimensionless pressure and its derivative decrease as PFA increases, indicating that a large PFA is beneficial to improve the productivity of fractured horizontal wells. In addition, as PFA increases, the first radial flow regime occurs later. The effect of PFA on other flow regimes can be neglected.   Figure 9. Rate distribution of a horizontal well with three reorientation fractures (k x = 4k y ). Figures 8 and 9 for a horizontal well with three reorientation fractures.

Basic Model Parameters Value
Drainage area (dimensionless) 200 × 200 Fracture conductivity (dimensionless) 20 Principal fracture angle (deg) 20 Inner fracture position (dimensionless) x D = y D = 100 The production rate distribution of each reorientation fracture along the horizontal well in an isotropic tight reservoir is illustrated in Figure 8 when k x = k y . In the early flow regime (the black line in Figure 8), the production rate of each fracture segment is symmetrically distributed with respect to the wellbore node. As the pressure wave expands, the production rate in the reoriented fracture significantly increases along the extension of the reoriented fracture (the green line in Figure 8). For a long production time (Figure 8, t D = 10 4 ), the production rate distribution changes significantly. Since the outermost fractures (RF1 and RF3) have a larger drainage area, they have a higher production rate. In addition, owing to the flow interference, the inner fracture (RF2) has a lower production rate. Furthermore, the production rate distribution of the outermost fractures (RF1 and RF3) (the red line in Figure 8) is asymmetrical. The reoriented fractures away from the horizontal wellbore have a larger production rate due to the larger seepage area caused by fracture reorientation. Figure 9 presents the production rate distribution for each reorientation fracture in a horizontal well in an anisotropic tight reservoir (k x = 4k y , FN = 3). The production rate distribution within each reorientation fracture is similar to that in Figure 8. However, for an anisotropic reservoir, fracture reorientation has a more significant effect on the production rate in each reorientation fracture. In the reorientation section, the production rate distribution curves vary sharply (Figure 9), which suggests that the effect of fracture reorientation in anisotropic tight reservoir on type curves is significant and, thus, cannot be neglected.

Parameter Influence on Transient Pressure Behavior
In this section, we analyze the sensitivity of the transient pressure of multiple reorientation fractures along a horizontal well in an anisotropic tight reservoir. We consider the case of a horizontal well with two reorientation fractures to analyze the effects of some key parameters, such as principal  Figure 10 presents the effect of PFA on type curves for a horizontal well with two reorientation fractures. The dimensionless fracture spacing is set at 20. As displayed in Figure 10, PFA has a weak impact on type curves in the early-time period, i.e., the bilinear flow regime, first linear flow regime, and first radial flow regime. For an anisotropic tight reservoir, the dimensionless pressure and its derivative decrease as PFA increases, indicating that a large PFA is beneficial to improve the productivity of fractured horizontal wells. In addition, as PFA increases, the first radial flow regime occurs later. The effect of PFA on other flow regimes can be neglected.  Figure 11 shows the effect of RFA on the dimensionless pressure and pressure derivative curves. In order to completely reveal the effect of RFA on type curves, we compare the dimensionless pressure (PFA = RFA = π/3) with the dimensionless pressure (PFA RFA) and then find the point 'MRPC' (maximum relative pressure change) ( Figure 12). The relative pressure change (RPC) is defined as follows: In Figure 12, t DMRPC represents the corresponding dimensionless time for the MRPC. Before t DMRPC , as the RFA increases, the RPC decreases which means that RFA should be larger to maintain the constant rate under the same pressure depletion rate. However, after t DMRPC , the effect of RFA on RPC weakens, which indicates that the effect of RFA on type curves tends to disappear. When RFA > PFA, RPC gets smaller and even becomes negative, which suggests that the larger the RFA (when RFA > PFA), the higher the well productivity.
curves. In order to completely reveal the effect of RFA on type curves, we compare the dimensionless pressure (PFA = RFA = π/3) with the dimensionless pressure (PFA ≠ RFA) and then find the point 'MRPC' (maximum relative pressure change) ( Figure 12). The relative pressure change (RPC) is defined as follows: In Figure 12, tDMRPC represents the corresponding dimensionless time for the MRPC. Before tDMRPC, as the RFA increases, the RPC decreases which means that RFA should be larger to maintain the constant rate under the same pressure depletion rate. However, after tDMRPC, the effect of RFA on Figure 11. Effect of the reoriented fracture angle on pressure and the pressure derivative. RPC weakens, which indicates that the effect of RFA on type curves tends to disappear. When RFA > PFA, RPC gets smaller and even becomes negative, which suggests that the larger the RFA (when RFA > PFA), the higher the well productivity.

Effect of Permeability Anisotropy on Type Curves
The impact of permeability anisotropy on dimensionless pressure and pressure derivative is illustrated in Figure 13 (β > 1) and Figure 14 (β < 1). In Figures 13 and 14, permeability anisotropy affects all the flow regimes except the pseudo-steady-state flow regime.
The influence of permeability anisotropy (β ≥ 1) on type curves is displayed in Figure 13. Here, the horizontal wellbore is parallel to the principal permeability axis. As permeability anisotropic factor (β) increases, the dimensionless pressure and its derivative decrease in the early linear flow regimes, indicating that the pressure difference of the horizontal wellbore can be smaller to maintain the constant rate. However, before the pseudo-steady-state flow regime, the pressure and its derivative increase when the permeability anisotropic factor (β) increases, for the reason that the first radial flow regime tends to disappear when permeability anisotropic factor (PAF) increases to some extent.

Effect of Permeability Anisotropy on Type Curves
The impact of permeability anisotropy on dimensionless pressure and pressure derivative is illustrated in Figure 13 (β > 1) and Figure 14 (β < 1). In Figures 13 and 14, permeability anisotropy affects all the flow regimes except the pseudo-steady-state flow regime.
The influence of permeability anisotropy (β ≥ 1) on type curves is displayed in Figure 13. Here, the horizontal wellbore is parallel to the principal permeability axis. As permeability anisotropic factor (β) increases, the dimensionless pressure and its derivative decrease in the early linear flow regimes, indicating that the pressure difference of the horizontal wellbore can be smaller to maintain the constant rate. However, before the pseudo-steady-state flow regime, the pressure and its derivative increase when the permeability anisotropic factor (β) increases, for the reason that the first radial flow regime tends to disappear when permeability anisotropic factor (PAF) increases to some extent. Similarly, the effect of permeability anisotropic factor (β ≤ 1) on type curves is shown in Figure  14. At this time, the horizontal wellbore is perpendicular to the principal permeability axis. Compared with Figure 13, we can find that the second radial flow regime tends to disappear when the horizontal wellbore is perpendicular to the principal permeability axis. Meanwhile, the compound linear flow regime tends to appear when the first radial flow regime ends, for the reason that the flow in the principal permeability direction dominates and, thus, the second radial flow regime is converted into the compound linear flow regime.  Figure 15 depicts the pressure and pressure derivative curves for a horizontal well intercepted by two reorientation fractures with various fracture conductivities: π, 10π, 100π, and 1000, respectively. Dimensionless pressure, p wD Dimensionless pressure derivative, dp wD Figure 13. Effect of permeability anisotropy on pressure and the pressure derivative (β ≥ 1). Figure 13. Effect of permeability anisotropy on pressure and the pressure derivative (β ≥ 1).

Effect of Fracture Conductivity on Type Curves
Similarly, the effect of permeability anisotropic factor (β ≤ 1) on type curves is shown in Figure  14. At this time, the horizontal wellbore is perpendicular to the principal permeability axis. Compared with Figure 13, we can find that the second radial flow regime tends to disappear when the horizontal wellbore is perpendicular to the principal permeability axis. Meanwhile, the compound linear flow regime tends to appear when the first radial flow regime ends, for the reason that the flow in the principal permeability direction dominates and, thus, the second radial flow regime is converted into the compound linear flow regime.  Figure 15 depicts the pressure and pressure derivative curves for a horizontal well intercepted by two reorientation fractures with various fracture conductivities: π, 10π, 100π, and 1000, respectively. Dimensionless pressure, p wD Dimensionless pressure derivative, dp wD Similarly, the effect of permeability anisotropic factor (β ≤ 1) on type curves is shown in Figure 14. At this time, the horizontal wellbore is perpendicular to the principal permeability axis. Compared with Figure 13, we can find that the second radial flow regime tends to disappear when the horizontal wellbore is perpendicular to the principal permeability axis. Meanwhile, the compound linear flow regime tends to appear when the first radial flow regime ends, for the reason that the flow in the principal permeability direction dominates and, thus, the second radial flow regime is converted into the compound linear flow regime. Figure 15 depicts the pressure and pressure derivative curves for a horizontal well intercepted by two reorientation fractures with various fracture conductivities: π, 10π, 100π, and 1000, respectively. Here PFA and RFA are kept constant and equal to 60° and 30°, respectively. As seen in Figure  15, the fracture conductivity affects the transient pressure significantly at the early-time period, mainly focusing on bilinear flow regime and first linear flow regime. The lower the fracture conductivity is, the earlier the bilinear flow regime appears. Moreover, when the fracture conductivity is low enough, the first linear flow regime becomes shorter and even tends to disappear. Meanwhile, the bilinear flow regime disappears and consequently the duration of the first linear flow lasts longer when fracture conductivity tends to be infinite (FcD > 100π).

Effect of Dimensionless Fracture Spacing (FS) on Type Curves
Fracture spacing (FS) determines the intensity of flow interferences happened between the adjacent fractures [10], and can mask some typical flow behaviors. For example, the first radial flow regime disappears when FS is smaller than the fracture length. When FS is large enough, the flow interference occurs later, or each fracture flows independently and does not interfere with each other.
The influence of FS on type curves is presented in Figure 16 and it can be concluded that the difference among the various type curves is mainly concentrated on the duration of first radial flow regime and pseudo-radial flow regime. Smaller fracture spacing generates larger pressure depletion to maintain the constant production rate due to the strong flow interference between the adjacent fractures. Moreover, the larger the fracture spacing is, the later the pseudo radial flow regime appears and also the shorter this regime's duration is. For multiple fractured horizontal wells, the optimal fracture spacing mainly depends on the reorientation fracture number and the drainage area. x wD = y wD = 100 x eD = y eD = 200 β = 2 FN=2 FS = 20 θ 1,1 = θ 2,1 = π /3 θ 1,2 = θ 2,2 = π /6 θ 1,3 = θ 2,3 = π /6 Dimensionless pressure, p wD Dimensionless pressure derivative, dp wD Figure 15. Effect of fracture conductivity on pressure and the pressure derivative.
Here PFA and RFA are kept constant and equal to 60 • and 30 • , respectively. As seen in Figure 15, the fracture conductivity affects the transient pressure significantly at the early-time period, mainly focusing on bilinear flow regime and first linear flow regime. The lower the fracture conductivity is, the earlier the bilinear flow regime appears. Moreover, when the fracture conductivity is low enough, the first linear flow regime becomes shorter and even tends to disappear. Meanwhile, the bilinear flow regime disappears and consequently the duration of the first linear flow lasts longer when fracture conductivity tends to be infinite (F cD > 100π).

Effect of Dimensionless Fracture Spacing (FS) on Type Curves
Fracture spacing (FS) determines the intensity of flow interferences happened between the adjacent fractures [10], and can mask some typical flow behaviors. For example, the first radial flow regime disappears when FS is smaller than the fracture length. When FS is large enough, the flow interference occurs later, or each fracture flows independently and does not interfere with each other.
The influence of FS on type curves is presented in Figure 16 and it can be concluded that the difference among the various type curves is mainly concentrated on the duration of first radial flow regime and pseudo-radial flow regime. Smaller fracture spacing generates larger pressure depletion to maintain the constant production rate due to the strong flow interference between the adjacent fractures. Moreover, the larger the fracture spacing is, the later the pseudo radial flow regime appears and also the shorter this regime's duration is. For multiple fractured horizontal wells, the optimal fracture spacing mainly depends on the reorientation fracture number and the drainage area. Energies 2019, 12, x FOR PEER REVIEW 17 of 23 Figure 16. Effect of dimensionless fracture spacing on pressure and the pressure derivative.

Effect of Fracture Number (FN) on Type Curves
Fracture spacing between the adjacent fractures is assumed to be constant (FS = 20). The fracture number can control the area between fractures and the reservoir. Figure 17 presents the influence of fracture number (FN) on pressure transient. The increasing FN means an increasing fractured area when the adjacent fracture spacing is fixed. In general, the pressure drop decreases with the increase of the FN. As seen in Figure 17, bilinear flow, first linear flow, first radial flow and second radial flow are all affected by reorientation fracture number. A large fracture number corresponds to a small pressure and derivative, suggesting that a smaller pressure drop is required to remain the same production rate. Consequently, the production is definitely enhanced by the increase in the fracture number, reflecting the advantage of multi-stage fractured horizontal wells. However, the degree of wD & dp wD FN= 1 FN= 3 FN= 5 FN= 7 x wD = y wD = 100 x eD = y eD = 200 β = 2 FS=20 F cD = 20 θ 1,1 = θ 2,1 = π /3 θ 1,2 = θ 2,2 = π /4 θ 1,3 = θ 2,3 = π /6 Dimensionless pressure, p wD Dimensionless pressure derivative, dp wD  As seen in Figure 17, bilinear flow, first linear flow, first radial flow and second radial flow are all affected by reorientation fracture number. A large fracture number corresponds to a small pressure and derivative, suggesting that a smaller pressure drop is required to remain the same production rate. Consequently, the production is definitely enhanced by the increase in the fracture number, reflecting the advantage of multi-stage fractured horizontal wells. However, the degree of wD & dp wD FN= 1 FN= 3 FN= 5 FN= 7 x wD = y wD = 100 x eD = y eD = 200 β = 2 FS=20 F cD = 20 θ 1,1 = θ 2,1 = π /3 θ 1,2 = θ 2,2 = π /4 θ 1,3 = θ 2,3 = π /6 Dimensionless pressure, p wD Dimensionless pressure derivative, dp wD Figure 17. Effect of fracture number on pressure and the pressure derivative.
As seen in Figure 17, bilinear flow, first linear flow, first radial flow and second radial flow are all affected by reorientation fracture number. A large fracture number corresponds to a small pressure and derivative, suggesting that a smaller pressure drop is required to remain the same production rate. Consequently, the production is definitely enhanced by the increase in the fracture number, reflecting the advantage of multi-stage fractured horizontal wells. However, the degree of the decrease of the pressure drop weakens when FN increases to some degree, indicating that the optimal reorientation fracture number needs to be demonstrated to maintain economic production.  For large fracture spacing, the influence of fracture configuration on type curves is largely focused on the transition flow regime (10 -5 ≤ tD ≤ 10 -1 ). Due to large fracture spacing, the difference between the effects of the staggered form of the reorientation fractures and the asymmetrical reorientation fractures on type curves can be neglected because each reorientation fracture behaves individually and does not interfere with each other. However, symmetrical reorientation fractures have higher production rates than the asymmetrical reorientation fractures when FS is large. For large fracture spacing, the influence of fracture configuration on type curves is largely focused on the transition flow regime (10 -5 ≤ t D ≤ 10 -1 ). Due to large fracture spacing, the difference between the effects of the staggered form of the reorientation fractures and the asymmetrical reorientation fractures on type curves can be neglected because each reorientation fracture behaves individually and does not interfere with each other. However, symmetrical reorientation fractures have higher production rates than the asymmetrical reorientation fractures when FS is large.

Effect of Complex Fracture Configuration on Type Curves
To determine the impact of fracture configuration on type curves under small fracture spacing (FS = 4), three typical cases are considered ( Figure 19). Compared with the cases under large fracture spacing (Figure 18), high production rates are obtained when the reorientation fractures are staggered because the mutual dislocation of the reorientation fractures can reduce the flow interference between the adjacent reorientation fractures.
Further, the influence of different fracture configurations under the same fracture spacing (FS = 20) and fracture number (FN = 3) on type curves is illustrated in Figure 20. Case A depicts that the length of the outermost fractures are larger than that of the inner fractures, Case B represents that the length of the outermost fractures are smaller than that of the inner fractures, and Case C shows that each fracture in the system has the same length. Therefore, the equal fracture length model has a low pressure depletion under the same conditions, which suggests that during the hydraulic fracturing treatment the length of each fracture should be equal to maintain economic production ( Figure 20).
For large fracture spacing, the influence of fracture configuration on type curves is largely focused on the transition flow regime (10 -5 ≤ tD ≤ 10 -1 ). Due to large fracture spacing, the difference between the effects of the staggered form of the reorientation fractures and the asymmetrical reorientation fractures on type curves can be neglected because each reorientation fracture behaves individually and does not interfere with each other. However, symmetrical reorientation fractures have higher production rates than the asymmetrical reorientation fractures when FS is large. To determine the impact of fracture configuration on type curves under small fracture spacing (FS = 4), three typical cases are considered ( Figure 19). Compared with the cases under large fracture spacing (Figure 18), high production rates are obtained when the reorientation fractures are staggered because the mutual dislocation of the reorientation fractures can reduce the flow interference between the adjacent reorientation fractures.  Figure 20. Case A depicts that the length of the outermost fractures are larger than that of the inner fractures, Case B represents that the length of the outermost fractures are smaller than that of the inner fractures, and Case C shows that each fracture in the system has the same length. Therefore, the equal fracture length model has a low pressure depletion under the same conditions, which suggests that during the hydraulic fracturing treatment the length of each fracture should be equal to maintain economic production ( Figure 20).

Conclusions
On the basis of the new fracture equation, we present a semi-analytical solution to investigate the pressure behavior of multiple reorientation fractures along the horizontal well in an anisotropic tight reservoir. The main conclusions include the following:

Conclusions
On the basis of the new fracture equation, we present a semi-analytical solution to investigate the pressure behavior of multiple reorientation fractures along the horizontal well in an anisotropic tight reservoir. The main conclusions include the following: (1) For a horizontal well perforated by multiple reorientation fractures with finite-conductivity in an anisotropic rectangular reservoir, six typical flow regimes can be found on type curves, including bilinear flow regime, first linear flow regime, first radial flow regime, second radial flow regime, pseudo-radial flow regime, and pseudo-steady-state flow regime. Meanwhile, the occurrence and duration of these typical regimes are determined by some significant parameters, such as the permeability anisotropic factor, fracture reorientation, fracture spacing, fracture number, and fracture configuration. (2) Fracture reorientation has an obvious effect on the production rate distribution along the fracture extension, particularly in the anisotropic tight reservoir. The production rate distribution curves for a horizontal well with three reorientation fractures in an anisotropic tight reservoir oscillate, and the outermost reorientation fracture has a larger production rate than the inner reorientation fracture due to a larger contact area with the reservoir. (3) Fracture reorientation is one of the highlights of this work. For an anisotropic tight reservoir, the horizontal well should be deployed parallel to the direction of the principal permeability axis, and all reoriented fractures should be perpendicular to the direction of the horizontal well. (4) The influence of permeability anisotropy on type curves is significant in all typical flow regimes except the late-time period. When horizontal wellbore is parallel to the principal permeability axis, first radial flow regime tends to disappear. However, when the horizontal wellbore is perpendicular to the principal permeability axis, the pseudo-radial flow regime changes into a compound linear flow regime. (5) During the fracturing treatment, it is necessary to make the adjacent fractures stagger over each other and make the fracture spacing large to reduce the flow interference, and the length of each fracture in multiple fractured horizontal wells should remain the same to maintain economic production.
Author Contributions: All authors have contributed to this work. G.X. proposed the mathematical model and wrote the main manuscript text. As the corresponding author, M.W. made substantial contributions to the conception/design of the work. S.W., H.W., B.W., and J.C. contributed to the modeling, programming, results analysis, and discussion. J.W. mainly focuses on the modeling, programming of the manuscript. Acknowledgments: I am greatly indebted to my supervisor, Shuhong Wu, for her valuable instructions and suggestions on this work as well as her careful reading of the manuscript. In addition, I deeply appreciate the contribution to this work made in various ways by other co-authors.

Conflicts of Interest:
The authors declare that they have no competing interests.

Appendix A Simplified Algorithm for Discrete Equations of a Reorientation Fracture
The discretized form of the reorientation fracture is: For i > j, Equation (A1) can be rewritten as: For i =j, Equation (A1) is given as: Otherwise, Equation (A1) can be expressed as: p wD −p f Di = 2π F cD l Di ∆l D jq f D j (A4)