Research Article  Open Access
Junchao Wang, Haitao Li, Yongqing Wang, Ying Li, Beibei Jiang, Wei Luo, "A New Model to Predict Productivity of MultipleFractured Horizontal Well in Naturally Fractured Reservoirs", Mathematical Problems in Engineering, vol. 2015, Article ID 892594, 9 pages, 2015. https://doi.org/10.1155/2015/892594
A New Model to Predict Productivity of MultipleFractured Horizontal Well in Naturally Fractured Reservoirs
Abstract
In order to predict productivity of multiplefractured horizontal well in fractured reservoir, flow models of reservoir and hydraulic fractures based on the volumetric source idealization are developed. The models are solved by utilizing Laplace transformation and orthogonal transformation, and flow rate of the well is calculated by coupling the two models. Compared to traditional point source functions, volumetric source function has many advantages in properties of function and programming calculation. The productivity predicting model is verified via an analytical ternaryporosity model. Moreover, a practical example of fractured horizontal well is studied to analyze the productivity and its influent factors. The result shows that flow rate of each fracture is different and inner fracture contributes least to productivity. Meanwhile, there are optimizing ranges for number, length, and conductivity of hydraulic fractures. In lowpermeability reservoir, increasing surface area in contact with reservoir by increasing number and length of hydraulic fractures is the most effective method to improve the productivity.
1. Introduction
Hydraulic fracturing is widely applied to improve the productivity of horizontal wells. Analytical methods to predict the productivity of fractured horizontal well have become invalid tools to optimize the hydraulic fracturing and determine productivity of horizontal well. There are two main methods to accomplish the target, which are analytical flow models [1–6] and source function (Green function) models [7–16].
Analytical flow models such as triplelinear model [1, 2] and ternaryporosity model [3] were built to simulate productivity of fractured horizontal well. But these models are limited in application [17, 18]; the triplelinear model is not accurate if the regions beyond the well tips control the well performance; the triporosity model could not simulate the properties of natural fractures very well.
In order to obtain more practical solutions, point source functions [19–21] for fractured reservoirs are developed to simulate pressure behavior and productivity. Application of source functions (Green’s functions) is the main technique since Gringarten and Ramey Jr. [22] introduced methods to solve transient flow problems for the first time. They also studied the application of Newman’s principle to transform 3dimensional system into the product of three 1dimensional systems. As for simplification, subsequent works have been focusing on point source functions which neglected the flow in source and its volume. Ogunsanya et al. [23, 24] pointed out that point source functions have computational problems and the no volume assumption may not be adequate in some specific situations, and then they revisited the concept of volumetric source and presented a new “solid bar source” model. Amini et al. [11–14] developed distributed volumetric source (DVS) model of homogeneous reservoirs and validated the applications in vertical well, horizontal well, and fractured vertical well with wellknown existing solutions. Amini and Valkó [14] also demonstrated the application of DVS model in fractured horizontal well and emphasized advantages of the method for volumetric source. The solution of volumetric source model does not suffer from the problem of inherent singularity, and pressure can be calculated at any point. Meanwhile, it is more adjustable than point source solutions; the solution is always in the same form, no matter what the types of the source are.
This study is focused on the derivation and application of volumetric source model which represents pressure behavior or productivity is generalized for fractured reservoir; a new steady volumetric source model is applied to calculate pressure drop in hydraulic fractures for finiteconductivity assumption. Laplace transformation and orthogonal transformation are adopted to derive the solutions of two models. Productivity of fractured horizontal well is calculated by coupling the two models. An example of fractured horizontal well is studied; factors that influence the well productivity are analyzed.
2. Volumetric Source Model
2.1. Assumption
The porous media are assumed to be dualpermeability reservoir whose shape is a box with closed boundaries, and the size of the reservoir is . A volumetric source with strength produces oil in the reservoir, for which the central coordinate of the source is , and the size is . Figure 1 shows the schematic of the system.
2.2. Mathematical Model
In dualporosity idealizations, Warren and Root introduced the dualporosity model in terms of bulk properties. If represents an intrinsic property of medium (matrix) or (natural fractures), then bulk property of medium is where is the ratio of volume of medium , .
Characteristics of the matrix and natural fractures system are incorporated by the storativity and flow capacity ratio is defined as follows: The pressure drop in the matrix or fracture system from initial pressure is Dimensionless directions and time are
Volumetric source model is where is Heaviside function:
2.3. Solution of the Model
Take Laplace transformation in model (5): Volumetric source model after transformation is where
Corresponding to model (9), define the orthogonal transformation as Inverse transformation formula is where characteristic function of the orthogonal transformation is characteristic value is and norm of characteristic function is
After orthogonal transformation, the solution of the model is Applying the inverse transformation formula (12), volumetric source function in Laplace Space is where presents pressure drop at arbitrary point when volumetric source is created with unit strength
3. Flow Model of Hydraulic Fracture
3.1. Assumption of Flow in Hydraulic Fracture
The length of hydraulic fracture is , height is , and width is . The flow in hydraulic fracture only occurs on the plane of the fracture. Oil from reservoir flows into the well through the hydraulic fracture immediately, and the pressure of hydraulic fracture recovers to reach balance at each time step.
Divide hydraulic fracture into segments ( segments in direction and segments in direction, both and are even which is shown in Figure 2). The th segment is described as , in which central coordinate is , total rate is , and volume is (). The central coordinate of horizontal well which intersected by the hydraulic fracture is , and the volume is (, ).
3.2. Pressure Drop in Hydraulic Fracture
For each time step, pressure drop in hydraulic fracture is where is the average hydraulic fracture pressure at time .
Assume that only the th segment flows to the well; flow model of hydraulic fracture is where
Corresponding to model (20), the orthogonal transformation is defined as inverse transformation formula is where characteristic function of the orthogonal transformation is characteristic value is and norm of characteristic function is
Through application of the inverse transformation formula (23), pressure drop in hydraulic fracture is where
Pressure difference from arbitrary point of hydraulic fracture to the well is where According to the superposition theory, when all the segments of hydraulic fracture flow to the well, pressure difference from arbitrary point to the well is Take Laplace transformation to (31); the pressure difference in Laplace Space is
4. Productivity of Fractured Horizontal Well
The th hydraulic fracture is divided into segments, and the size is . Then () volumetric source exists in the reservoir, matrix of each source’s rate in Laplace Space is , 3D coordinate is , and 2D coordinate is .
According to the pressure drop caused by flowing in reservoir and hydraulic fracture, the problem can be written in the matrix formation as where The matrix of each source’s rate can be calculated, and the total flow rate of the well in Laplace Space is According to Stehfest algorithm, the flow rate of the well in real time space can be calculated numerically.
5. Influential Factors Analysis
Volumetric source function in this paper is generalized from the case in homogeneous reservoirs developed by Valkó and Amini [11], and validation of volumetric source function has also been tested. In the paper, the volumetric source model will be verified by analytical flow model [6]. Table 1 shows the basic parameter values of this simulation.

Figure 3 shows the schematic of the simulated example, five evenly spaced hydraulic fractures are distributed along the horizontal well, and each hydraulic fracture is divided into volumetric sources. Productiontime relationship simulated by two methods is shown in Figure 4, which shows great consistency except for the very beginning of production. As the analytical flow model was proposed based on linear flow assumption, the flow resistance is underestimated at the beginning of production, which leads to result of analytical flow model which is higher than the result of volumetric source models.
Figure 5 implies that the contribution of each hydraulic fracture is different. Hydraulic fractures at the end of horizontal well dominate large portion of reservoir and have less interference with other fractures. Therefore, they have greater contribution than fractures in middle part. Ratio of flow rate shown in Figure 5 is defined as
5.1. Factors of Hydraulic Fractures
In this subsection, influences of factors of hydraulic fractures on productivity are analyzed. In order to reveal the effects of hydraulic fracture and reservoir factors, incremental rate () and ratio of cumulative production () are adopted: where represents factor value, such as value of number or length or conductivity; represents the minimal factor’s value simulated; represents cumulative production of horizontal well.
With increasing number, the flow rate increases rapidly at the beginning, but interference between hydraulic fractures shows up earlier, and flow rate decreases faster with time. Figure 6 presents the ratios of cumulative production and increment rate for different number of hydraulic fractures. As illustrated, increment rate is less than 0.1 as the number of hydraulic fracture exceeds 45; cumulative production is more than 3 times larger than that of a single hydraulic fracture.
Figure 7 reveals effect of length on productivity of horizontal well. One of the major advantages of hydraulic fracturing is increasing the surface area in contact reservoir. For a certain number, it is accomplished by increasing the length of hydraulic fracture. Result shows that increment rate is less than 0.1 if length of hydraulic fracture approaches 250 m.
Another mechanism for hydraulic fracturing to improve productivity of horizontal well is creating highconductivity flow path. In this case, conductivity should remain at 20 cm · μm^{2}–30 cm · μm^{2} (see Figure 8).
Comparing the results of Figures 6–8, it is obvious that the increasing number improves productivity most significantly, while increasing conductivity improves the least. In unconventional reservoirs, flow capacity of reservoirs may be not sufficient to supply the hydraulic fractures to its full flow capacity. Therefore, higher conductivity of hydraulic fractures is undesirable once basic flow capacity of reservoir is achieved.
The simulated example reveals a clear method to design hydraulic fracturing by applying volumetric source model. For a real horizontal well, parameters of hydraulic fractures could be optimally designed by following the method.
5.2. Permeability of Reservoir
Flow capacity of reservoir is the primary factor that affects designing of hydraulic fracturing, which is controlled by natural fractures in the dualporosity idealization. Figures 9–11 show the effects of hydraulic fractures on ratios of cumulative production when permeability of reservoir is 0.001 μm^{2}–0.1 μm^{2}. When permeability of reservoir reduces the range of lowpermeability, oil flows with more difficulty, and interference between hydraulic fractures will show up much later. Therefore, increasing the surface area in contacted reservoir seems more important, and the highconductivity is less vital for hydraulic fractures, which can easily satisfy the flow capacity of reservoir. This outcome emphasizes the importance of improving contacted area when designing hydraulic fracturing.
6. Conclusions
(1)New volumetric source models of flow in reservoir and hydraulic fractures are developed to simulate the productivity of fractured horizontal well in naturally fractured reservoir, which is verified by analytical flow model.(2)Compared with traditional point source functions, the volumetric source function brings many advantages. Firstly, it shows no inherent singularity, and pressure can be calculated at any point in reservoir. Secondly, it is more adjustable; the solution is always in the same form, regardless of the types of the source. Last, it is more convenient for programming calculation.(3)A practical example is simulated to show the essence of productivity of fractured horizontal wells. Result shows that flow rate of each fracture is different and the inner fracture contributes least to productivity. Increasing number improves productivity most significantly, while increasing length is the second important factor, and increasing conductivity is the least.(4)The defined parameters and could be employed to optimize number, length, and conductivity of hydraulic fractures. More and longer hydraulic fractures are required for lower permeability reservoirs, but larger conductivity is unnecessary.
Nomenclature
(m):  Length of reservoir 
(m):  Width of reservoir 
(m):  Height of reservoir 
(m):  Position of the center of the source in direction 
(m):  Position of the center of the source in direction 
(m):  Position of the center of the source in direction 
(m):  Source width in direction 
(m):  Source width in direction 
(m):  Source width in direction 
:  Strength of the source 
(dimensionless):  The ratio of volume 
(m^{−2}):  Shape factor 
(D):  Permeability 
(MPa^{−1}):  Storage 
:  Viscosity 
(m):  Reference length 
(MPa):  Initial pressure 
(MPa):  Pressure 
(dimensionless):  Flow capacity ratio 
(dimensionless):  Storativity. 
:  Matrix 
:  Natural fractures 
:  Dimensionless 
:  Bulk property. 
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This study is supported by the National Science and Technology Major Project of China (Grant no. 2011ZX05022006004HZ).
References
 M. Brown, E. Ozkan, R. Raghavan, and H. Kazemi, “Practical solutions for pressure transient responses of fractured horizontal wells in unconventional reservoirs,” in Proceedings of the SPE Annual Technical Conference and Exhibition, New Orleans, La, USA, 2009, SPE 125043. View at: Publisher Site  Google Scholar
 M. Brown, E. Ozkan, R. Raghavan, and H. Kazemi, “Comparison of fracturedhorizontalwell performance in tight sand and shale reservoirs,” SPE Reservoir Evaluation & Engineering, vol. 14, no. 2, pp. 248–259, 2014. View at: Publisher Site  Google Scholar
 A. Hasan, A tripleporosity model for fractured horizontal wells [M.S. thesis], Texas A&M University, 2010.
 X. Wang, W. Luo, X. Hou, and J. Wang, “Transient pressure analysis of multiplefractured horizontal wells in boxed reservoirs,” Petroleum Exploration and Development, vol. 41, no. 1, pp. 74–78, 2014. View at: Publisher Site  Google Scholar
 L. I. Longlong, Y. Jun, Y. Li, M. Wu, Q. Zeng, and R. Lu, “Productivity calculation and distribution of staged multicluster fractured horizontal wells,” Petroleum Exploration and Development, vol. 41, no. 4, pp. 504–508, 2014. View at: Publisher Site  Google Scholar
 L. Haitao, W. Junchao, W. Yongqing, J. Beibei, and L. Wei, “Generalized productivity model for designing hydraulic fractures in horizontal wells located in naturally fissured lowpermeability gas reservoirs,” Chemistry and Technology of Fuels and Oils, vol. 50, no. 4, pp. 299–314, 2014. View at: Publisher Site  Google Scholar
 R. Raghavan and C. ChihCheng, “An analysis of horizontal wells intercepted by multiple fractors,” SPE Journal, vol. 2, no. 3, pp. 235–245, 1997. View at: Publisher Site  Google Scholar
 C.C. Chen and R. Raghavan, “A multiplyfractured horizontal well in a rectangular drainage region,” SPE Journal, vol. 2, no. 4, pp. 455–462, 1997. View at: Publisher Site  Google Scholar
 A. Zerzar and Y. Bettam, “Interpretation of multiple hydraulically fractured horizontal well in closed systems,” in Proceedings of the Canadian International Petroleum Conference, PETSOC2004027, June 2004. View at: Publisher Site  Google Scholar
 J. Guo, F. Zeng, J. Zhao, and Y. Xu, “A new model to predict fractured horizontal well production,” in Proceedings of the Canadian International Petroleum Conference, Calgary, Canada, June 2006. View at: Publisher Site  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 Hydraulic Fracturing Technology Conference, SPE 106279, pp. 426–439, January 2007. View at: Google Scholar
 D. Zhu, F. Magalhaes, and P. P. Valko, “Predicting productivity of multiplefractured horizontal gas wells,” in Proceedings of the SPE Hydraulic Fracturing Technology Conference, SPE 106280, College Station, Tex, USA, January 2007. View at: Publisher Site  Google Scholar
 F. Magalhaes, D. Zhu, S. Amini, and P. P. Valko, “Optimization of fracturedwell performance of horizontal gas wells,” in Proceedings of the SPE 2nd International Oil Conference and Exhibition, SPE 108779, pp. 503–509, Veracruz, Mexico, June 2007. View at: Publisher Site  Google Scholar
 S. Amini and P. P. Valkó, “Using distributed volumetric sources to predict production from multiplefractured horizontal wells under nondarcyflow conditions,” SPE Journal, vol. 15, no. 1, pp. 105–115, 2010. View at: Publisher Site  Google Scholar
 P. Q. Lian, L. S. Cheng, and J. Y. Cui, “A new computation model of fractured horizontal well coupling with reservoir,” International Journal for Numerical Methods in Fluids, vol. 67, no. 8, pp. 1047–1056, 2011. View at: Publisher Site  Google Scholar
 J. Lin and D. Zhu, “Predicting well performance in complex fracture systems by slab source method,” in Proceedings of the Hydraulic Fracturing Technology Conference, SPE 151960, The Woodlands, Tex, USA, February 2012. View at: Publisher Site  Google Scholar
 S. K. Siddiqui, A. Ali, and H. Dehghanpour, “New advances in production data analysis of hydraulically fractured tight reservoirs,” in Proceedings of the SPE Canadian Unconventional Resources Conference, SPE 162830, Calgary, Canada, OctoberNovember 2012. View at: Publisher Site  Google Scholar
 A. Alkouh, D. Schechter, R. A. Wattenbarger et al., “Practical use of simulators for characterization of shale reservoirs,” in Proceedings of the SPE Canadian Unconventional Resources Conference, SPE162645MS, OctoberNovember 2012. View at: Publisher Site  Google Scholar
 E. Ozkan and R. Raghavan, “New solutions for welltestanalysis problems. Part 1. Analytical considerations,” SPE Formation Evaluation, vol. 6, no. 3, pp. 359–368, 1991. View at: Publisher Site  Google Scholar
 E. Ozkan and R. Raghavan, “New solutions for welltestanalysis problems. Part 2. Computational considerations and applications,” SPE Formation Evaluation, vol. 6, no. 3, pp. 369–378, 1991. View at: Publisher Site  Google Scholar
 E. Ozkan and R. Raghavan, “New solutions for welltestanalysis problems: part IIIadditional algorithms,” in Proceedings of the Annual Technical Conference and Exhibition, SPE 28424, New Orleans, La, USA, September 1994. View at: Publisher Site  Google Scholar
 A. C. Gringarten and H. J. Ramey Jr., “The use of source and Green’s functions in solving unsteadyflow problems in reservoir,” Society of Petroleum Engineers Journal, vol. 13, no. 5, pp. 285–296, 1973. View at: Publisher Site  Google Scholar
 B. O. Ogunsanya, T. P. Oetama, J. F. Lea, L. R. Heinze, and P. S. Adisoemarta, “A coupled model for analyzing transient pressure behavior of horizontal drainholes,” in Proceedings of the SPE Production and Operations Symposium, pp. 619–632, April 2005. View at: Publisher Site  Google Scholar
 B. O. Ogunsanya, T. P. Oetama, J. F. Lea, L. R. Heinze, and P. S. Adisoemarta, “A robust type curve solution for analyzing pressuretransient behaviors of both vertical and horizontal fractures systems,” in Proceedings of the Annual SPE International Technical Conference and Exhibition, SPE 105979, Abuja, Nigeria, August 2006. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Junchao 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.