Research Article  Open Access
Productivity Formulae of an InfiniteConductivity Hydraulically Fractured Well Producing at Constant Wellbore Pressure Based on Numerical Solutions of a Weakly Singular Integral Equation of the First Kind
Abstract
In order to increase productivity, it is important to study the performance of a hydraulically fractured well producing at constant wellbore pressure. This paper constructs a new productivity formula, which is obtained by solving a weakly singular integral equation of the first kind, for an infiniteconductivity hydraulically fractured well producing at constant pressure. And the two key components of this paper are a weakly singular integral equation of the first kind and a steadystate productivity formula. A new midrectangle algorithm and a Galerkin method are presented in order to solve the weakly singular integral equation. The numerical results of these two methods are in accordance with each other. And then the solutions of the weakly singular integral equation are utilized for the productivity formula of hydraulic fractured wells producing at constant pressure, which provide fast analytical tools to evaluate production performance of infiniteconductivity fractured wells. The paper also shows equipotential threads, which are generated from the numerical results, with different fluid potential values. These threads can be approximately taken as a family of ellipses whose focuses are the two endpoints of the fracture, which is in accordance with the regular assumption in Kuchuk and Brigham, 1979.
1. Introduction
The main objective of hydraulic fracturing for well stimulation is to increase well productivity by creating a highly conductive path some distance away from the skin zone into the formation. The fracture creates more surface area to the wellbore without drilling another well. Since more reservoir area is in direct communication with the wellbore, a greater volume of fluid can be produced into the wellbore per unit time, resulting in an increased production rate. This basic objective has not changed since hydraulic fracturing was introduced in the early 1950s. At the beginning, hydraulic fracturing was considered a good mean to increase productivity of the wells completed in lowpermeability reservoirs. Now, it has become an integral part of most well completions.
For hydraulically fractured wells, unsteadystate pressure transient testings are useful tools for evaluating insitu reservoir and wellbore parameters that describe the production characteristics. The use of transient well testing for determining reservoir parameters and productivity of fractured wells has become common, and determination of transient pressure behavior and productivity for fractured wells has aroused considerable interest over the past decades.
Numerous analytic solutions have been presented for the pressure behavior of a fractured well producing at constant flow rate. An extensive literature survey on fractured wells producing at constant flow rates can be found. Guppy et al. present for the first time a technique that analyzes buildup and drawdown data from wells produced at constant pressure with turbulent flow in the fracture [1]. CincoLey and SamaniegoV studied the earlytime pressure data for a well intercepted by a finiteconductivity vertical fracture [2]. Ozkan et al. described the characteristics of a well producing at a constant pressure in a natural fractured reservoir [3]. Guppy and CincoLey presented semianalytic solutions for unsteadystate flow behavior of a well intersecting a vertical fracture and analytic solutions for defining certain portions of the early time data for various types of fracture conductivity [4]. Nashawi discussed a semianalytical equation which incorporates the effects of nonDarcy flow in the fracture [5]. Kuchuk and Brigham presented analytical solutions to elliptical flow problems that are applicable to infiniteconductivity vertically fractured wells, elliptically shaped reservoirs, and anisotropic reservoirs producing at constant rate or pressure [6].
Although most well test analysis methods for hydraulically fractured wells assume constant rate production, it must be pointed out that the constant rate production is difficult to maintain, and constant wellbore pressure production conditions are not uncommon. Examples of conditions under which constant pressure is maintained at a well include production into a constant pressure separator or pipeline, open flow to the atmosphere, or production from a low permeability reservoir. It is important to consider wells producing at constant pressure rather than constant rate during large portions of the production life of tight reservoirs [4].
This paper proposes a new productivity model of an infiniteconductive hydraulically fractured well producing at constant pressure. The weakly singular integral equation [7–9] arising from this model has attracted much attention for its numerical solutions [10–12].
In order to solve the problem of constant wellbore pressure production, the numerical solution to a weakly singular integral equation is required, and a new midrectangular algorithm is used to solve the weakly singular integral equation; the corresponding numerical results are compared with those by the Galerkin method.
Furthermore, the solutions of the weakly singular integral equation are used in the productivity formulae of hydraulically fractured wells producing at constant pressure. These formulae provide fast analytical tools to evaluate production performance of infiniteconductivity fractured wells. And equipotential threads with different fluid potential values, which can be approximately taken as a family of ellipses whose focuses are the two endpoints of the fracture, are also shown.
The rest of this paper is organized as follows. In Section 2, the new model is derived. In Section 3, the Galerkin method and the corresponding numerical results are presented. In Section 4, the midrectangle algorithm is constructed, and the corresponding numerical experiments are presented. Then conclusions are reached in Section 5.
2. InfiniteConductivity Hydraulically Fractured Well Model Producing at Constant Wellbore Pressure
In this section, we will construct a new productivity model of an infiniteconductivity hydraulically fractured well producing at constant wellbore pressure.
Figure 1 represents an ideal vertical fracture. The usual assumptions apply; that is, the porous medium is isotropic, horizontal, homogeneous, and uniform in thickness and has constant permeability. Also, the fracture fully penetrates the vertical extent of the formation and is the same length on both sides of the well.
The porous media domain is as follows: where is cylinder radius and is the cylindrical body. Since both the vertical well and the hydraulic fracture are fully penetrating, we may use twodimensional model to study the pressure behavior. Assume the fracture length is , the fracture is taken as a line sink in the twodimensional space, and the coordinates of the two ends are and . Supposing that point is on the fracture, in order to obtain the pressure at point caused by the point , we have to obtain the basic solution of the following partial differential equation in : where is the pressure at point . is the formation permeability, is viscosity, is the flow rate at point , are Dirac functions and that is, the hydraulic fracture is between .
Let where is the pressure at infinite distance away from the fracture. Because the hydraulically fractured well is infinite conductivity, the pressure of the fracture denoted by is uniform. That is
Thus, if is known, the total productivity is given by
In order to simplify the previously mentioned equations, we take the following dimensionless transforms:
Obviously, we obtain
Since then
In the same manner, we have So That is, where
Note that after taking the dimensionless transforms, producing section of the fractured well is changed to , and the initial conditions are changed as follows:
It is well known that the fundamental solution of twodimensional Laplace equation is as follows [13, 14]:
Thus, if is known, the solution to (2.13) is given by
If the hydraulically fractured well is infinite conductivity, then we have where is a constant.
From (2.17), we can obtain the following equation: and it can be simplified as follows: Let
Consequently, we obtain the following first kind weakly singular integral equation:
Once is obtained from the previous equation, we can use it to compute
Using we can obtain
The total productivity of the hydraulically fractured well is as follows:
Remark 2.1. Equation (2.22) is the first kind weakly singular integral equation, which can be solved with the numerical methods introduced in Sections 3 and 4. And the total productivity can be obtained by (2.26).
Remark 2.2. From Tables 1 and 2 in Section 3 or Tables 3 and 4 in Section 4, we can find that is a constant. So the total productivity formula is as follows: where is the length of the hydraulically fractured well.




3. Galerkin Method for SteadyState Productivity Computation of the Hydraulically Fractured Well
In this section, we will present the Galerkin method, which has a mature analysis framework, to solve the weakly singular integral equation (2.22). Once of (2.22) is known, the approximation of (2.26) is also presented.
Let
Clearly,
Let
Define the energy inner product and the energy norm
The weak formulation is to find such that
And the Galerkin formulation is to find such that
Then the matrix representation of the approximation is where
Once the solution of (3.9) is obtained, we can utilize it to obtain the total productivity of the hydraulically fractured well by the following formula:
In fact which completes the proof of Formula (3.11).
In the following, we will present the error estimate of the Galerkin approximation. Since , then (3.6) and (3.8) lead to
Hence where is the orthogonal projection with respect to the energy inner product. Consequently, there holds Using (3.15), the following theorem, which is recalled from [7, 8], is valid.
Theorem 3.1 (see [7, 8]). If is sufficiently smooth, g is a smooth function on , and is a subspace which is construct of piecewise smooth constant functions, then there is a constant such that
Now we turn to the implementation issue of (3.9). Clearly,
Using (3.17), we can obtain
In the numerical experiment, we choose and in the previous algorithm in order to obtain the numerical solutions of (2.22). Table 1 shows the Galerkin approximations with , from which we can see that the solution is symmetric, that is, . And Table 2 shows the total productivity of fractured well with different , which is convergent when . More observations and conclusions from the numerical results will be presented at the end of Section 4.
4. Midrectangle Algorithm for SteadyState Productivity Computation of the Hydraulically Fractured Well
In this section, we will construct and analyze a midrectangle algorithm, which is much more convenient to implement but still has similar accuracy to that of the Galerkin method.
We first prove the following theorem, which will lead to the algorithm and provide the error estimate of the method.
Theorem 4.1. For any , there holds
Proof. Let
Obviously,
By Taylor expansion, we can obtain
That is,
Hence
Here
where we have used the result of . Now we consider . Let . There exists such a constant that
Plugging and into (4.6), the proof is completed.
Let
Clearly,
Let be the approximation of , and let be the singular points in (2.22). Then Hence we obtain the following theorem.
Theorem 4.2. , , there holds
From Theorem 4.2, we can construct the following midrectangle algorithm for (2.22):
The error estimate of this method is given by Theorem 4.2. And the matrix representation of (4.13) is where
After in (4.14) is obtained, the total productivity of the hydraulically fractured well can be approximated as follows:
Table 3 shows the midrectangle approximations with . And Table 4 shows the total productivity of fractured well with different . The numerical results perform similarly to those of the Galerkin method presented in Tables 1 and 2.
From Tables 1 and 3, we can find that when a hydraulically fractured well is producing at constant wellbore pressure, the point convergence intensity at endpoint is bigger than that of midpoint, and the point convergence intensity is at monotone decreasing trend from endpoint to midpoint.
From Tables 2 and 4, we can see that the numerical results converge to
Furthermore, this constant is independent of well length, permeability, viscosity, and other fluid properties. Hence this result can be used to calculate flow rates for many realistic cases as follows:
Here is the length of the hydraulically fractured well, is permeability, is viscosity is pressure at an infinite distance from the well, and is wellbore pressure of the hydraulically fractured well. From (4.18), we can see that when a hydraulically fractured well is producing at constant wellbore pressure, its flow rate is directly proportional to well length, permeability, and pressure difference.
If is a constant, the productivity index is given by
Using (2.17) and the results in Table 1, we can calculate the pressure at a point as follows:
And the points, which satisfy form equipotential threads when the well is producing at constant wellbore pressure. Figure 2 shows equipotential threads with different fluid potential values. It can be found that the equipotential threads can be approximately taken as a family of ellipses whose focuses are the two endpoints of the hydraulic fractured well. This conclusion is in accordance with the regular assumption in [6].
5. Conclusions
This paper proposes a new model using a weakly singular integral equation for the productivity of infiniteconductivity hydraulically fractured wells producing at constant wellbore pressure. A Galerkin method and a midrectangle algorithm are constructed to solve this integral equation. And their numerical results are similar to each other. Then the numerical solutions of this equation are utilized in the productivity formula, which provide fast analytical tools to evaluate production performance of infiniteconductivity fractured wells. Furthermore, the following conclusions are obtained.(1)The numerical results for the weakly singular integral equation are convergent. Moreover, the point convergence intensity is symmetric. And the point convergence intensity is monotonically decreasing from the endpoints to the midpoint when a hydraulically fractured well is producing at constant wellbore pressure. (2)The numerical approximation for the productivity is convergent when .(3)The equipotential threads, which are generated from the numerical results, can be approximately taken as a family of ellipses whose focuses are the two endpoints of the hydraulically fractured well. This conclusion is in accordance with the regular assumption in [6].
References
 K. H. Guppy, S. Kumar, and V. D. Kagawan, “Pressuretransient analysis for fractured wells producing at constant pressure,” SPE Formation Evaluation, vol. 3, no. 1, pp. 169–178, 1988. View at: Google Scholar
 H. CincoLey and F. SamaniegoV, “Transient pressure analysis for fractured wells,” Society of Petroleum Engineers of AIME, vol. 33, no. 9, pp. 1749–1766, 1981. View at: Google Scholar
 E. Ozkan, U. Ohaeri, and R. Raghavan, “Unsteady flow to a well produced at a constant pressure in a fractured reservoir,” SPE Formation Evaluation, vol. 2, no. 2, pp. 186–200, 1987. View at: Google Scholar
 K. H. Guppy and H. CincoLey, “Transient flow behavior of a vertically fractured well producing at constant pressure,” Society of Petroleum Engineers of AIME, vol. 21, no. 12, pp. 9962–9963, 1981. View at: Google Scholar
 I. S. Nashawi, “Constantpressure analysis of infiniteconductivity fractured gas wells influenced by nonDarcy flow effects,” in Proceedings of the SPE 15th Middle East Oil and Gas Show and Conference (MEOS '07), pp. 260–271, Manama, Bahrain, March 2007. View at: Google Scholar
 F. Kuchuk and W. F. Brigham, “Transient flow in elliptical system,” Society of Petroleum Engineers of AIME, vol. 19, no. 6, pp. 401–410, 1979. View at: Google Scholar
 I. H. Sloan and A. Spence, “The Galerkin method for integral equations of the first kind with logarithmic kernel: theory,” IMA Journal of Numerical Analysis, vol. 8, no. 1, pp. 105–122, 1988. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 I. H. Sloan and A. Spence, “The Galerkin method for integral equations of the first kind with logarithmic kernel: applications,” IMA Journal of Numerical Analysis, vol. 8, no. 1, pp. 123–140, 1988. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 I. G. Graham, “Galerkin methods for second kind integral equations with singularities,” Mathematics of Computation, vol. 39, no. 160, pp. 519–533, 1982. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 K. OravPuurand, A. Pedas, and G. Vainikko, “Nyström type methods for Fredholm integral equations with weak singularities,” Journal of Computational and Applied Mathematics, vol. 234, no. 9, pp. 2848–2858, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. Huang, T. Lü, and Z. C. Li, “Mechanical quadrature methods and their splitting extrapolations for boundary integral equations of first kind on open arcs,” Applied Numerical Mathematics, vol. 59, no. 12, pp. 2908–2922, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Z. Rui, H. Jin, and L. Tao, “Mechanical quadrature methods and their splitting extrapolations for solving boundary integral equations of axisymmetric Laplace mixed boundary value problems,” Engineering Analysis with Boundary Elements, vol. 30, no. 5, pp. 391–398, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products, Academic Press, San Diego, Calif, USA, 7th edition, 2007.
 M. Fogiel, Handbook of Mathematical, Scientific, and Engineering, Research and Education Association, Piscataway, NJ, USA, 1994.
Copyright
Copyright © 2012 Chaolang Hu 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.