Research Article  Open Access
A Semianalytical Solution for Multifractured Horizontal Wells in BoxShaped Reservoirs
Abstract
This paper presented a 3D point sink model through using Dirac function. Then, 3D point sink solution in boxed reservoirs was obtained through using Laplace transform and Fourier transform methods. Based on the flux and pressure equivalent conditions in Laplace space, a semianalytical solution for multifractured horizontal wells was also proposed for the first time. The production rate distribution was discussed in detail for multifractured horizontal wells. The calculative results show the outermost fractures had higher production ratio due to larger drainage area and the inner fractures were lower due to the strong interface between fractures. Type curves were established to analyze the flow characteristics, which would be divided into six stages, for example, bilinear flow region, the first linear flow region, the first radial flow region, the second linear flow region, the second radial flow region, and the boundary dominated flow region, respectively. Finally, effects of some sensitive parameters on type curves were also analyzed in detail.
1. Introduction
The development of tight reservoirs has been paid more and more attention in China. Stimulation for a horizontal well in tight reservoirs may further enhance well productivity. Unlike a vertical well, a horizontal well may be more than one point along the well length. During the last two decades, horizontal wells have become a common applied completion in the petroleum industry. With a large reservoir contact area, horizontal wells can greatly improve well productivity; however, it is most advantageous to drill horizontal well in thin and tight reservoirs with vertical fractures [1, 2].
The main difficulty of transientpressure analysis for multifractured horizontal wells (MFHWs) is how to solve those problems of multiplefracture interference, finite conductivity transverse fractures coupled with the formation flow. Some typical research papers are worth reviewing [2–9]. Since 1972, some attempts have been made to simulate the pressure transient behavior for either horizontal or vertical wells, with or without hydraulic fractures. Larsen and Hegre [3] rigorously presented a transient pressure solution of horizontal wells with circular finiteconductivity fractures in the threedimensional unbounded formation. The results just showed the earlystage and middlestage features of fracture system because the effects of outer boundary were not considered into their model. Guo et al. [2] presented the pressure transient behavior for a horizontal well with multiple randomly distributed vertical fractures in the infinite reservoir and bounded reservoir. Wan and Aziz [6] described a new semianalytical solution for horizontal wells with multiple hydraulic fractures. The fractures can be rotated at any horizontal angle to the well and need not fully penetrate the formation in the vertical direction. AlKobaisi et al. [7] presented a hybrid numerical analytical model with a finiteconductivity vertical fracture intercepted by a horizontal well, which dynamically couples a numerical fracture model with an analytical reservoir model. Their approach allows us to include finer details of the fractures characteristics while keeping the computational work manageable. Valkó and Amini [8] proposed a DVS (distributed volumetric sources) method to predict gas production for a horizontal well with multiple transverse fractures in a bounded reservoir. But it was only an approximate approach and the fracture conductivity was not considered in their model. AI Rbeawi and Tiab [9] introduced a new technique for interpreting the pressure transient behavior for a horizontal well with multiple infinite conductivity fractures which could be longitudinal or transverse, vertical or inclined, or symmetrical or asymmetrical. Recently, some meshless methods have been presented [10–13] and were possibly applied to solve the fracture problems.
The main objective of this paper is to develop a computational model to investigate the transientpressure behavior for multifractured horizontal wells (MFHWs) in boxshaped reservoirs. A semianalytical solution for MFHWs in boxed reservoirs is obtained through using Laplace transform and Fourier transform. The fluid flow in the fracture system and reservoir system is computed separately, and the flux and pressure equivalent conditions in Laplace space are applied in the fracture wall to couple the fluids flow in both systems. The result is validated accurately through comparing with previous results in the literature.
2. Physical Model
Figures 1 and 2 show the physical model assumed and the assumptions of this model are as follows.(1)A singlephase flow was assumed in the reservoir. There is a fractured horizontal well in the boxed, closed, and homogeneous formation whose length is , width is , and height is (see Figure 1).(2)The well produces with constant viscosity and slightly compressible fluid at the total flow rate of , but the flow rates from each fracture may change over time.(3)The oil production is assumed to be an isothermal process. Darcy flow is assumed in the fracture. The linear flow occurs in the fracture (see Figure 2).(4)Two wings of each fracture may be of unequal length, and the tips of the fractures are assumed as noflow boundaries (see Figure 2).(5)The horizontal wellbore is assumed to be infinitely conductive and the fractures are assumed to be finitely conductive.(6)The effect of conductivity can be described by the common vertically fractured wells model.(7)At the starting time of production, the pressure is uniformly distributed and is equal to the initial pressure ().
3. Mathematical Model
3.1. Reservoir Model
To obtain the fracture solution in the boxed reservoirs, we must firstly obtain the point sink solution. The partial differential equation which is closed in all directions in the Cartesian coordinate system is given by where is porosity, is total compressibility, is the time, is reservoir pressure, is reservoir permeability, is point flux, is Dirac function, (, , ) is the point position, and and are the boundary length and width. All the boundaries of the boxshaped drainage volume are closed; therefore, boundary conditions can be given by And the initial condition can be given by In order to simplify (1)–(3), we take the following dimensionless transforms: where is the characteristic length, is the point convergence intensity at the point sink (, , ), and the total productivity of the well is ; there holds Then, through using dimensionless transforms, (1) becomes Equation (2) becomes And initial condition becomes The Laplace transform is based on and functions as follows: Applying the Laplace transforms to (6)–(8), we have Outer boundary conditions in Laplace space are Fourier cosine transform based on can be defined as The characteristic equation is Through solving (13), characteristic number is obtained by Imposing the Fourier cosine transform on the variable in (10), we have Similarly, imposing the Fourier cosine transform on the variable in (15), we have Similarly, imposing the Fourier cosine transform on the variable in (16), we have where
Imposing the Fourier inverse transform on the variable in (17a) and (17b), we obtain Imposing the Fourier inverse transform on the variable in (18), we have Imposing the Fourier inverse transform on the variable in (19), we have where we use the following formulas: Therefore, we have and the cosine product formula can be given by Through using (22)(23), the following formulas can be obtained: Substituting (21)–(24) into (20), we can have Further, (25) could be simplified as the following equation: Equation (26) is the point sink solution in boxshaped reservoirs.
3.2. Fracture Model
Fluid flow inside th fracture in the system may be defined by the following set of equations: where where is the well position in the fracture, is the fracture permeability, is well rate in the fracture, is the fracture length, is the source position, is the reservoir thick, and is the fracture flux. In order to simplify (27), we take the following dimensionless transforms: Therefore, the following equations can be obtained by The initial condition can be given by The boundary conditions are given by Imposing Laplace transforms on in (30)–(32), we have The boundary conditions become
3.3. Semianalytical Solution
The reservoir and fractureflow equations are described in Sections 3.1 and 3.2. To build the reservoirflow model, Fourier transformation on space and Laplace transformation on time have been applied and point sink solution has been also obtained in Section 3.1. Thus, line sink solution can be derived by integrating the solution given in (26) with respect to over 0 to . (Note that only the integral of cos () needs to be determined.) The result can be given by where
Fracture number for MFHW is ; the pressure drop formulas at each fracture segment for MFHW can be obtained through using superposition principle and point sink integral method, where is the angle between th fracture and axis, is defined in (36), and is source position in the fracture. Now, we divide per fracture into segments; thus, pressure drop formulas for th segment of th fracture can be given by To obtain the solution for the fractureflow models, Laplace transformations on both space and time have been carried out. The pressuredrop evaluation for a fracture is presented in the following equation. The pressure drop at th segment on the th fracture is given by where where is the Heaviside unit step function.
For each fracture, the total flux equation can be given by and the pressure for each fracture is equal to the pressure of horizontal well; it must have In addition, total rate of horizontal wells should be equal to sum of flux for each fracture; thus, Coupling the fracture and reservoirflow models at the fracture faces, we obtain a system of equations with unknowns. is the total number of fracture segments, and is the total number of fractures. This system of equations could be rearranged and written in the matrix format. Once the matrix is solved for unknowns, the rate at each fracture segment and the pressure drop at each fracture will be obtained. These rates and pressure drops are used to calculate the pressure drop at the wellbore for MFHW.
4. Results and Discussion
4.1. Validation of Solutions in This Paper
Riley et al. [14] presented an analytical solution for a single fracture case. To validate the solutions proposed in this paper, some data were utilized (see Table 1) for comparing the solutions with those from Riley et al. [14] (see Figure 3). For the convenience of the comparison, only a fracture is considered into the MFHWs. As shown in Figure 3, there is a good agreement between the solutions obtained in this work and the results presented by Riley. This might indicate that previous models are special cases of our MFHWs model and also implies that our model is validated accurately.

4.2. The Production Rate Distribution
The production rate distribution is presented in Figure 4 and the data is shown in Table 2. In this case, five transversal fractures have identical length and uniform distribution along a horizontal well located in a boxshaped reservoir. As shown in Figure 4, the outermost fractures have higher production ratio due to larger drainage area and the inner fractures are lower due to strong interference between fractures. For a specific fracture, the flux distribution along the discrete elements has the same property as the whole fracture system.

4.3. Flow Characteristics Analysis
The transient transport characteristics are graphically showed by type curves, which can be used to analyzing transient pressure and rate decline so as to recognize the flow characteristics of fluids in the reservoir. In addition, by type curves matching some reservoir property parameters, such as permeability, skin factor, gas in place, fracture halflength, and gas reservoir drainage area, can be obtained [15, 16].
The pressure and its derivative curves are presented in Figure 5, which shows flow characteristics for a fractured horizontal well with 5 identical lengths and uniformdistribution fractures, and the used data is shown in Table 2. As shown in Figure 5, the flow can be divided into six stages, which could be corresponding with Figure 6 and described as follows.
(a)
(b)
(c)
(d)
(e)
(f)
Stage 1. In this stage, the segment has a straight line with 1/4 slope, reflecting the bilinear flow region (see Figure 6(a)). In this region, fluids flow from fracture to wellbore and from reservoirs to fracture at the same time. This region only occurs when the fracture conductivity is relatively small.
Stage 2. The segment has a 1/2 slope straight line on both pressure and pressure derivative curves, namely, the first linear flow region. In this region, flow occurs linearly and directly from formation to individual fractures and each fracture behaves independently of the other fractures (see Figure 6(b)). The linear flow is the optimal flow mode because it can reduce the seepage resistance.
Stage 3. If the pressure derivative curve shows the 1/() constant, the first radial flow region will be observed. In this region, radial flow occurs directly (see Figure 6(c)) from formation to individual fractures and the continuous time of radial flow for each fracture depends on the fracture length and fracture spacing. A big fracture length or big spacing between the fractures will lead to a short time for radial flow.
Stage 4. The segment has a 0.36 slope straight line on both pressure and pressure derivative curves, namely, the second linear flow region. The reservoir linear flow occurs in this region (Figure 6(d)). Fluids flow will start to occur when the pressure wave spreads farther. The linear flow starts from reservoir to the vertical plane which contains the horizontal wellbore.
Stage 5. The segment has a straight line with 0.5 constant, namely, the second radial flow region (Figure 6(e)). During this flow period, flow across the outer most producing elements becomes dominant.
Stage 6. The segment has a unit slope straight line on both pressure and pressure derivative curves, namely, boundary dominated flow region. This flow period occurs when the reservoir boundary is reached (Figure 6(f)). Due to the fracturedhorizontal wells centered in the reservoir, the fluids flow reached the six boundaries nearly at the same time. If the fracturedhorizontal wells were not centered in the reservoir, the characteristics of boundary dominated flow may be not a unit slope straight line on both pressure and pressure derivative curves.
The proposed flow regions above may not exist in a single test. Depending on the specific properties of fracture and reservoir, some of the flow regions may be absent. For the moderate ratios of fracture length and spacing, the boundary dominated flow region may be nonexistent or replaced by a transitional flow region. It is also possible that radial flow region may not exist because of well interference or boundary effects. Admittedly, affected by fracture conductivity and relative position, some characteristics of flow pattern may not clearly emerge at times. In the real case, engineers should analyze the results of flow regions according to the concrete case.
4.4. Effects of Sensitive Parameters on Type Curves
Based on the presented model, the main factors affecting transientpressure characteristics of fractured horizontal wells in coal reservoirs are presented in Figures 7–9. Those parameters are including fracture number, fracture conductivity , and reservoir width . The basic used parameters are listed in Table 2.
The effect of fracture number is shown in Figure 7 and fracture number can be set to 3 and 7, respectively. It appears that the fracture number mainly affects the flow characteristics of bilinear flow region, the first linear flow region, and the first radial flow region, which could be seen on the pressure derivative curve. It can be seen from pressure derivative curve that the time of the first radial flow region with big fracture number is longer than that with small fracture number and the time of the second linear flow region with small fracture number is longer than that with big fracture number. It is known that the production mainly benefits from the second linear flow region, and the longer the second linear flow region is, the better the production is. Therefore, through the analysis above, increasing the fracture number is helpful to production, but it is not the more the fracture, the better the development. It can be analyzed from the pressure curve that pressure curve with big fracture number is always higher than the one with small fracture number, which means the increasing fractures number will cause small pressure depletion.
The effect of fracture conductivity is shown in Figure 8 and fracture conductivity can be set to 6 and 7, respectively. The derivative curve shows that the effect of the fracture conductivity on flow characteristics occurs obviously in the bilinear flow region. As shown in the figure, the bilinear flow characteristic will gradually disappear as the conductivity increases; finally, linear flow will appear when . It can be also found from pressure curve that the pressure with small conductivity is higher than the one with big conductivity, which means fracture with high conductivity will cause big pressure depletion at early time.
The effect of reservoir width is shown in Figure 9 and reservoir width can be set to 1000 and 4000, respectively. Reservoir width mainly affects the boundary dominated flow region. It can be found from the pressure curve that pressure with big reservoir width is higher than the one with small reservoir width, which means that a small reservoir will cause big pressure depletion in the boundary dominated flow region. We also see from the pressure derivative curve that a small reservoir width value will shorten the time of the second radial flow and reach the boundary dominated flow region early.
5. Conclusions
Through using Fourier cosine transform and Laplace transform method, a semianalytical solution for multifractured horizontal wells in boxshaped reservoirs is presented to investigate the pressure transient behavior. The accuracy of the semianalytical solutions is verified by graphically comparing them with the published analytical solutions which are special cases of this paper. According to the presented solutions, the production rate distribution is also proposed for MFHWs, which show the outermost fractures have higher production ratio due to larger drainage area and the inner fractures are lower due to the strong interface between fractures. Type curves are established to analyze the flow characteristics, which can be divided into six stages. Effects of main sensitive parameters on type curves are analyzed in detail including fracture number, fracture conductivity, and reservoir width. The presented model could be applicable to develop well test soft for MFHWs, which will be our future research topic.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This paper was supported by the Major National Science and Technology Projects (Grants nos. 2011ZX05009004 and 2011ZX05013002).
References
 B. E. Roberts, H. van Engen, and C. P. J. W. van Kruysdijk, “Productivity of multiply fractured horizontal wells in tight gas reservoirs,” in Proceedings of the Offshore Europe Conference, pp. 133–141, September 1991, SPE paper no. 23113. View at: Google Scholar
 G. Guo, R. D. Evans, and M. M. Chang, “Pressuretransient behavior for a horizontal well intersecting multiple random discrete fractures,” in Proceedings of the SPE Annual Technical Conference & Exhibition, pp. 277–290, September 1994, SPE paper no. 28390. View at: Google Scholar
 L. Larsen and T. M. Hegre, “Pressure transient analysis of multifractured horizontal wells,” in Proceedings of the SPE Annual Technical Conference & Exhibition, pp. 265–276, September 1994, SPE paper no. 28389. View at: Google Scholar
 R. N. Horne and K. O. Temeng, “Relative productivities and pressure transient modeling of horizontal wells with multiple fractures,” in Proceedings of the 9th Middle East Oil Show & Conference, pp. 563–574, March 1995, SPE paper no. 29891. View at: Google Scholar
 C. C. Chen and R. Raghavan, “A multiplyfractured horizontal well in a rectangular drainage region,” in Proceedings of the 2nd International Conference on Horizontal Well Technology, pp. 311–319, November 1996, SPE paper no. 37072. View at: Google Scholar
 J. Wan and K. Aziz, “Semianalytical well model of horizontal wells with multiple hydraulic fractures,” SPE Journal, vol. 7, no. 4, pp. 437–445, 2002. View at: Google Scholar
 M. AlKobaisi, E. Ozkan, and H. Kazemi, “A hybrid numericalanalytical model of finiteconductivity vertical fractures intercepted by a horizontal well,” in Proceedings of the SPE International Petroleum Conference in Mexico, pp. 599–609, November 2004, SPE paper no. 92040. View at: Google Scholar
 P. P. Valkó and S. Amini, “The method of distributed volumetric sources for calculating the transient and pseudosteadystate productivity of complex wellfracture configurations,” in Proceedings of the SPE Hydraulic Fracturing Technology Conference, pp. 426–439, January 2007, SPE paper no. 106279. View at: Google Scholar
 S. AI Rbeawi and D. Tiab, “Transient pressure analysis of a horizontal well with multiple inclined hydraulic fractures using typecurve matching,” in Proceedings of the SPE International Symposium and Exhibition on Formation Damage Control, pp. 1–20, 2012, SPE paper no. 149902. View at: Publisher Site  Google Scholar
 P. Zhu, L. W. Zhang, and K. M. Liew, “Geometrically nonlinear thermomechanical analysis of moderately thick functionally graded plates using a local PetrovGalerkin approach with moving Kriging interpolation,” Composite Structures, vol. 107, pp. 298–314, 2014. View at: Google Scholar
 L. W. Zhang, P. Zhu, and K. M. Liew, “Thermal buckling of functionally graded plates using a local Kriging meshless method,” Composite Structures, vol. 108, pp. 472–492, 2014. View at: Google Scholar
 L. W. Zhang, Z. X. Lei, K. M. Liew, and J. L. Yu, “Static and dynamic of carbon nanotube reinforced functionally graded cylindrical panels,” Composite Structures, vol. 111, pp. 205–212, 2014. View at: Google Scholar
 K. M. Liew, Z. X. Lei, J. L. Yu, and L. W. Zhang, “Postbuckling of carbon nanotubereinforced functionally graded cylindrical panels under axial compression using a meshless approach,” Computer Methods in Applied Mechanics and Engineering, vol. 268, pp. 1–17, 2014. View at: Google Scholar
 M. F. Riley, W. E. Brigham, and R. N. Horne, “Analytic solutions for elliptical finiteconductivity fractures,” in Proceedings of the SPE Annual Technical Conference and Exhibition, pp. 35–48, October 1991, SPE paper no. 22656. View at: Google Scholar
 L. Wang, X. D. Wang, J. Q. Li, and J. Wang, “Simulation of pressure transient behavior for asymmetrically fractured wells in coal reservoirs,” Transport in Porous Media, vol. 97, no. 3, pp. 352–372, 2013. View at: Publisher Site  Google Scholar
 R.S. Nie, Y.F. Meng, J.C. Guo, and Y.L. Jia, “Modeling transient flow behavior of a horizontal well in a coal seam,” International Journal of Coal Geology, vol. 92, pp. 54–68, 2012. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Lei Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.