MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 428596 10.1155/2012/428596 428596 Research Article Productivity Formulae of an Infinite-Conductivity Hydraulically Fractured Well Producing at Constant Wellbore Pressure Based on Numerical Solutions of a Weakly Singular Integral Equation of the First Kind Hu Chaolang 1 Lu Jing 2 He Xiaoming 3 Chen Kue-Hong 1 College of Mathematics Sichuan University Sichuan, Chengdu 610041 China scu.edu.cn 2 Department of Petroleum Engineering The Petroleum Institute, P.O. Box 2533 Abu Dhabi UAE pi.ac.ae 3 Department of Mathematics and Statistics Missouri University of Science and Technology Rolla, MO 65401 USA mst.edu 2012 04 06 2012 2012 13 04 2012 29 05 2012 2012 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.

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 infinite-conductivity 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 steady-state 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 infinite-conductivity 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 low-permeability reservoirs. Now, it has become an integral part of most well completions.

For hydraulically fractured wells, unsteady-state 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 . Cinco-Ley and Samaniego-V studied the early-time pressure data for a well intercepted by a finite-conductivity vertical fracture . Ozkan et al. described the characteristics of a well producing at a constant pressure in a natural fractured reservoir . Guppy and Cinco-Ley presented semianalytic solutions for unsteady-state 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 . Nashawi discussed a semianalytical equation which incorporates the effects of non-Darcy flow in the fracture . Kuchuk and Brigham presented analytical solutions to elliptical flow problems that are applicable to infinite-conductivity vertically fractured wells, elliptically shaped reservoirs, and anisotropic reservoirs producing at constant rate or pressure .

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 .

This paper proposes a new productivity model of an infinite-conductive hydraulically fractured well producing at constant pressure. The weakly singular integral equation  arising from this model has attracted much attention for its numerical solutions .

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 infinite-conductivity 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. Infinite-Conductivity Hydraulically Fractured Well Model Producing at Constant Wellbore Pressure

In this section, we will construct a new productivity model of an infinite-conductivity 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.

Schematic of ideal fracture.

The porous media domain is as follows: (2.1)Ω={(x,y,z)x2+y2<Re2,0<z<H}, where Re is cylinder radius and Ω is the cylindrical body. Since both the vertical well and the hydraulic fracture are fully penetrating, we may use two-dimensional model to study the pressure behavior. Assume the fracture length is L,L=2xf, the fracture is taken as a line sink in the two-dimensional space, and the coordinates of the two ends are (-xf,0) and (xf,0). Supposing that point (x,0) is on the fracture, in order to obtain the pressure at point (x,y) caused by the point (x,0), we have to obtain the basic solution of the following partial differential equation in Ω: (2.2)Kμ2Px2+Kμ2Py2=-q(x)δ(x-x)δ(y), where P(x,y) is the pressure at point (x,y). K is the formation permeability, μ is viscosity, q(x) is the flow rate at point (x,0), δ(y),δ(x-x) are Dirac functions and (2.3)x[-xf,xf], that is, the hydraulic fracture is between [-xf,xf].

Let (2.4)P|r=Pi>0, where P|r 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 P0 is uniform. That is (2.5)P(x,0)=P0=constant,x[-xf,xf].

Thus, if q(x) is known, the total productivity is given by (2.6)Q=0H-xfxfq(x)dxdz=H-xfxfq(x)dx.

In order to simplify the previously mentioned equations, we take the following dimensionless transforms: (2.7)xD=x2xf,yD=y2xf,PD=Pi-PPi,PD0=Pi-P0Pi.

Obviously, we obtain (2.8)xDx=12xf,Px=(PxD)(xDx)=12xfPxD,2Px2=x(Px)=14xf22PxD2.

Since (2.9)PD=Pi-PPi, then (2.10)2Px2=-Pi4xf22PDxD2.

In the same manner, we have (2.11)2Py2=-Pi4xf22PDyD2,δ(x-x)δ(y)=δ[2xf(xD-xD)]δ(2xfyD)=(2xf)-2δ(xD-xD)δ(yD). So (2.12)Kμ2Px2+Kμ2Py2=-KμPi4xf2(2PDxD2+2PDyD2)=-q(xfxD)(2xf)-2δ(xD-xD)δ(yD). That is, (2.13)2PDxD2+2PDyD2=qD(xD)δ(xD-xD)δ(yD), where (2.14)qD(xD)=q(2xfxD)μKPi,xD[-12,12].

Note that after taking the dimensionless transforms, producing section [-xf,xf] of the fractured well is changed to [-1/2,1/2], and the initial conditions are changed as follows: (2.15)PD|rD0+,qD(xD)=0,xD[-12,12]  .

It is well known that the fundamental solution of two-dimensional Laplace equation is as follows [13, 14]: (2.16)g(x,y;x,y)=-ln[(x-x)2+(y-y)2]1/22π.

Thus, if qD(xD) is known, the solution to (2.13) is given by (2.17)PD(xD,yD)=-1/21/2g(xD,yD;xD,0)qD(xD)dxD.

If the hydraulically fractured well is infinite conductivity, then we have (2.18)PD(xD,0)=PD0,xD[-12,12], where PD0 is a constant.

From (2.17), we can obtain the following equation: (2.19)-1/21/2g(xD,0;xD,0)qD(xD)dxD=PD0, and it can be simplified as follows: (2.20)(-12π)-1/21/2ln|xD-xD|qD(xD)dxD=PD0. Let (2.21)u(xD)=qD(xD)2πPD0,xD[-12,12].

Consequently, we obtain the following first kind weakly singular integral equation: (2.22)-1/21/2ln|x-y|u(x)dx=-1,y[-12,12].

Once u(x) is obtained from the previous equation, we can use it to compute (2.23)qD(x)=2πPD0u(x),x[-12,12].

Using (2.24)qD(xD)=q(2xfxD)μKPi,xD[-12,12], we can obtain (2.25)q(2xfxD)=KPiμqD(xD),xD[-12,12].

The total productivity of the hydraulically fractured well is as follows: (2.26)Q=H-xfxfq(x)dx=2xfH-1/21/2q(2xfxD)dxD=2xf(KPiHμ)-1/21/2qD(xD)dxD=4xfπPD0(KPiHμ)-1/21/2u(x)dx=4xfπ(Pi-P0Pi)(KPiHμ)-1/21/2u(x)dx,Q=4xfπ(Pi-P0)KHμ-1/21/2u(x)dx.

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 -1/21/2u(x)dx is a constant. So the total productivity formula is as follows: (2.27)Q=4xfπ(Pi-P0)KHμ-1/21/2u(x)dx4xfπ(Pi-P0)KHμhj=0n-1uj4xfπ(Pi-P0)KHμ0.7213=9.0641xf(Pi-P0)KHμ=4.53205(Pi-P0)KHLμ, where L is the length of the hydraulically fractured well.

Galerkin approximations of (2.22) with h=1/32.

i 1 2 3 4
c i 2.6610334017 0.9490861410 0.8841883020 0.7384058892
c i + 4 0.6653019950 0.6118426439 0.5732934677 0.5440894309
c i + 8 0.5215521538 0.5039599049 0.4902132354 0.4795754716
c i + 12 0.4715448802 0.4657790288 0.4620502406 0.4602190183
c i + 16 0.4602190183 0.4620502406 0.4657790288 0.4715448802
c i + 20 0.4795754716 0.4902132354 0.5039599049 0.5215521538
c i + 24 0.5440894309 0.5732934677 0.6118426439 0.6653019950
c i + 28 0.7384058892 0.8841883020 0.9490861410 2.6610334017

Galerkin approximations of (3.11) with different h.

i 1 2 3 4
Q ( 1 / 2 i ) 0.6932919850 0.7068397093 0.7139761064
Q ( 1 / 2 i + 4 ) 0.7176334503 0.7194835363 0.7204138104 0.7208802386
Q ( 1 / 2 i + 8 ) 0.7211137731 0.7212306202 0.7212890637 0.7213182904

Mid-rectangle approximations of (2.22) with h=1/32.

i 1 2 3 4
c i 2.4652280029 1.1168311123 0.8795040869 0.7498254718
c i + 4 0.6700249966 0.6154959679 0.5759489664 0.5461808975
c i + 8 0.5215521538 0.5039599049 0.4902132354 0.4807028032
c i + 12 0.4725784047 0.4667483804 0.462979470 0.4611289341
c i + 16 0.4611289341 0.462979470 0.4667483804 0.4725784047
c i + 20 0.4807028032 0.4902132354 0.5039599049 0.5215521538
c i + 24 0.5440894309 0.5732934677 0.6154959679 0.6700249966
c i + 28 0.7498254718 0.8795040869 1.1168311123 2.4652280029

Mid-rectangle approximations of (4.16) with different h.

i 1 2 3 4
Q ( 1 / 2 i ) 0.6859163596 0.7044679397 0.7134577603
Q ( 1 / 2 i + 4 ) 0.7177072748 0.7196865358 0.7205981981 0.7210138248
Q ( 1 / 2 i + 8 ) 0.7212012467 0.7212846926 0.7213212664 0.7213369747
3. Galerkin Method for Steady-State 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 u(x) of (2.22) is known, the approximation of (2.26) is also presented.

Let (3.1)h=1n,xj=-12+jh,(j=0,,n),Γ=[-12,12],Γj=[xj-1,xj],(j=1,,n).

Clearly, (3.2)Γ=j=1nΓj.

Let (3.3)ej(x)={1,xΓj,0,xΓj,Sh=span{ej(x),j=1,,n}.

Define the energy inner product (3.4)[u,v]=(Ku,v)=-Γln|x-y|u(x)v(y)dsxdsy,u,vL2(Γ) and the energy norm (3.5)u2=[u,u].

The weak formulation is to find uL2(Γ) such that (3.6)[u,v]=(g,v),vL2(Γ).

And the Galerkin formulation is to find (3.7)uh(x)=j=1ncjej(x)Sh such that (3.8)[uh,vh]=(g,vh),vhSh.

Then the matrix representation of the approximation is (3.9)Ac¯=g¯, where (3.10)A=(aij)n×n,aij=-ΓiΓjln|x-y|dsxdsy,(i,j=1,,n),c¯=[c1,c2,,cn]T,g¯=[g1,g2,,gn]T,gj=(g,ej)=Γg(x)ej(y)dsxdsy,(j=1,,n).

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: (3.11)Q(h)=hj=1ncj.

In fact(3.12)Q(h)=-1/21/2uh(x)dx=-1/21/2j=1ncjej(x)dx=j=1ncj-1/21/2ej(x)dx=j=1ncjh, which completes the proof of Formula (3.11).

In the following, we will present the error estimate of the Galerkin approximation. Since ShL2(Γ), then (3.6) and (3.8) lead to (3.13)[uh-u,v]=0,vSh.

Hence (3.14)uh=Phu, where Ph:L2(Γ)Sh is the orthogonal projection with respect to the energy inner product. Consequently, there holds (3.15)uh-u=infvShu-v. Using (3.15), the following theorem, which is recalled from [7, 8], is valid.

Theorem 3.1 (see [<xref ref-type="bibr" rid="B7">7</xref>, <xref ref-type="bibr" rid="B8">8</xref>]).

If Γ is sufficiently smooth, g is a smooth function on Γ, and Sh is a subspace which is construct of piecewise smooth constant functions, then there is a constant c>0 such that (3.16)uh-uch2/3,|(uh-u,g)|ch3.

Now we turn to the implementation issue of (3.9). Clearly, (3.17)limh0hlnh=0.

Using (3.17), we can obtain (3.18)aij=-ΓiΓjln|x-y|dsxdsy=12(j-i-1)2h2ln|(j-i-1)h|-(j-i)2h2ln|(j-i)h|-32h2+12(j-i+1)2h2ln|(j-i+1)h|,gj=(g,ej)=Γg(x)ej(y)dsxdsy=Γg(x)dsxΓej(y)dsy=Γg(x)dsxh.

In the numerical experiment, we choose g(x)=-1 and Γ=[-1/2,1/2] in the previous algorithm in order to obtain the numerical solutions of (2.22). Table 1 shows the Galerkin approximations with h=1/32, from which we can see that the solution is symmetric, that is, ci=cn+1-i,i=1,,n. And Table 2 shows the total productivity of fractured well with different h, which is convergent when h0. More observations and conclusions from the numerical results will be presented at the end of Section 4.

4. Midrectangle Algorithm for Steady-State 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 g(x)C[t-a,t+a]2, there holds (4.1)t-at+ag(x)ln|x-t|dx=2g(t)a(lna-1)+O(a3lna).

Proof.

Let (4.2)h=an,xi=t-a+ih,(i=0,1,,n-1,n,n+1,,2n). Obviously, (4.3)xn=t,xi+x2n-i=2t. By Taylor expansion, we can obtain (4.4)g(xi)+g(x2n-i)=2g(t)+g(ζi)2(xi-t)2. That is, (4.5)g(xi)+g(x2n-i)=2g(t)+g(ζi)2((n-i)h)2. Hence (4.6)t-at+ag(x)ln|x-t|dx=limh0[hi=0n-1g(xi)ln|xi-t|+hi=n+12ng(xi)ln|xi-t|]=limh0[hi=0n-1g(xi)ln(a-ih)+hi=n+12ng(xi)ln(ih-a)]=limh0[hi=0n-1g(xi)ln(a-ih)+hi=n+1-2n2n-2ng(xi+2n)ln((i+2n)h-a)]=limh0[hi=0n-1g(xi)ln(a-ih)+hi=1-n0g(xi+2n)ln((i+2n)h-a)]=limh0[hi=0n-1g(xi)ln(a-ih)+hi=-(1-n)-0g(x-i+2n)ln((-i+2n)h-a)]=limh0[hi=0n-1g(xi)ln(a-ih)+hi=n-10g(x-i+2n)ln((-i+2n)h-a)]=limh0[hi=0n-1g(xi)ln(a-ih)+hi=0n-1g(x-i+2n)ln((-i+2n)h-a)]=limh0[hi=0n-1g(xi)ln(a-ih)+hi=0n-1g(x2n-i)ln(-ih+2nh-a)]=limh0[hi=0n-1g(xi)ln(a-ih)+hi=0n-1g(x2n-i)ln(a-ih)]=limh0[hi=0n-1(g(xi)+g(x2n-i))ln(a-ih)]=limh0[hi=0n-1(2g(t)+g(ζi)2((n-i)h)2))ln(a-ih)]=limh0[hi=0n-1(2g(t)ln(a-ih))]+limh0[hi=0n-1g(ζi)2((n-i)h)2ln(a-ih)]=limh0[hi=0n-1(2g(t)ln((n-i)h)]+limh0[hi=0n-1g(ζi)2((n-i)h)2ln((n-i)h)]=I1+I2. Here (4.7)I1=limh0[hi=0n-1(2g(t)ln((n-i)h)]=2g(t)limh0[hi=0n-1(ln((n-i)h)]=2g(t)limh0[hi=1nln(ih)]=2g(t)0alnxdx=2g(t)[xlnx-x]|0a=2g(t)[alna-a], where we have used the result of limx0xlnx=0. Now we consider I2. Let mg′′(ζi)M,(i=0,1,,n-1). There exists such a constant k that (4.8)I2=limh0[hi=0n-1g(ζi)2((n-i)h)2ln((n-i)h)]=k2limh0[hi=0n-1((n-i)h)2ln((n-i)h)]=k2limh0[hi=1n(ih)2ln(ih)]=k20ax2lnxdx=k2(x33lnx-x39)|0a=k2(a33lna-a39).

Plugging I1 and I2 into (4.6), the proof is completed.

Let (4.9)h=1n,xj=-12+jh,(j=0,1,,n),yj=-12+(j+12)h,(j=0,1,,n-1).

Clearly, (4.10)yj=12(xj+xj+1),(j=0,1,,n-1).

Let uj be the approximation of u(yj), and let yi=(1/2)(xi+xi+1),(i=0,1,,n-1) be the singular points in (2.22). Then (4.11)-1/21/2u(x)ln|x-yi|dx=j=0n-1xjxj+1u(x)ln|x-yi|dx=j=0i-1xjxj+1u(x)ln|x-yi|dx+xixi+1u(x)ln|x-yi|dx+j=i+1n-1xjxj+1u(x)ln|x-yi|dx=j=0i-1h[u(yj)ln|yj-yi|+O(h2)]+2u(yi)h2(lnh2-1)+O((h2)3lnh2)  +j=i+1n-1h[u(yj)ln|yj-yi|+O(h2)]=j=0,jin-1hu(yj)ln|yj-yi|+hu(yi)lnh2e+O(h2). Hence we obtain the following theorem.

Theorem 4.2.

For    all    u ( x ) C [ - 1 / 2 , 1 / 2 ]       ,    h = 1 / n , x j = ( - 1 / 2 ) + j h , j = 0,1 , , n ,  yj=-1/2+(j+1/2)h,(j=0,1,,n-1), there holds (4.12)-1/21/2u(x)ln|x-yi|dx=j=0,jin-1hu(yj)ln|yj-yi|+u(yi)hlnh2e+O(h2),(i=0,1,,n-1).

From Theorem 4.2, we can construct the following midrectangle algorithm for (2.22): (4.13)-1=-1/21/2u(x)ln|x-yi|dxj=0,jin-1(hln|yj-yi|)uj+(hlnh2e)ui,(i=0,1,,n-1).

The error estimate of this method is given by Theorem 4.2. And the matrix representation of (4.13) is (4.14)Au¯=-I, where (4.15)u¯=[u0,u1,,un-1]T,I=[1,1,,1]T,A=(aij)i,j=0n-1,ai,j={hlnh2ei=jhln(|j-i|h)ij,(i,j=0,1,,n-1).

After u¯ in (4.14) is obtained, the total productivity of the hydraulically fractured well can be approximated as follows: (4.16)Q(h)=hj=0n-1uj.

Table 3 shows the midrectangle approximations with h=1/32. And Table 4 shows the total productivity of fractured well with different h. 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 (4.17)0.7213-1/21/2u(x)dxQ(h)Q(h).

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: (4.18)Q=0.7213×2π×L×KHμ×(Pi-P0).

Here L is the length of the hydraulically fractured well, K is permeability, μ is viscosity Pi is pressure at an infinite distance from the well, and P0 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 (Pi-P0) is a constant, the productivity index is given by (4.19)Jv=0.7213×2π×L×KHμ=4.53205×KHLμ.

Using (2.17) and the results in Table 1, we can calculate the pressure at a point (x,y) as follows: (4.20)P(x,y)(h2π)i=031uiln((x-xi+0.5)2+y2).

And the points, which satisfy (4.21)P(x,y)(h2π)i=031uiln((x-xi+0.5)2+y2)=constant 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 .

Profile of equipotential threads of a hydraulically fractured well.

5. Conclusions

This paper proposes a new model using a weakly singular integral equation for the productivity of infinite-conductivity 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 infinite-conductivity fractured wells. Furthermore, the following conclusions are obtained.

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.

The numerical approximation for the productivity is convergent when h0.

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 .

Guppy K. H. Kumar S. Kagawan V. D. Pressure-transient analysis for fractured wells producing at constant pressure SPE Formation Evaluation 1988 3 1 169 178 2-s2.0-0023977980 Cinco-Ley H. Samaniego-V F. Transient pressure analysis for fractured wells Society of Petroleum Engineers of AIME 1981 33 9 1749 1766 2-s2.0-0019611240 Ozkan E. Ohaeri U. Raghavan R. Unsteady flow to a well produced at a constant pressure in a fractured reservoir SPE Formation Evaluation 1987 2 2 186 200 2-s2.0-0023362893 Guppy K. H. Cinco-Ley H. Transient flow behavior of a vertically fractured well producing at constant pressure Society of Petroleum Engineers of AIME 1981 21 12 9962 9963 2-s2.0-0019340471 Nashawi I. S. Constant-pressure analysis of infinite-conductivity fractured gas wells influenced by non-Darcy flow effects Proceedings of the SPE 15th Middle East Oil and Gas Show and Conference (MEOS '07) March 2007 Manama, Bahrain 260 271 2-s2.0-34250891556 Kuchuk F. Brigham W. F. Transient flow in elliptical system Society of Petroleum Engineers of AIME 1979 19 6 401 410 2-s2.0-0018323344 Sloan I. H. Spence A. The Galerkin method for integral equations of the first kind with logarithmic kernel: theory IMA Journal of Numerical Analysis 1988 8 1 105 122 10.1093/imanum/8.1.105 967846 ZBL0636.65143 Sloan I. H. Spence A. The Galerkin method for integral equations of the first kind with logarithmic kernel: applications IMA Journal of Numerical Analysis 1988 8 1 123 140 10.1093/imanum/8.1.123 967847 ZBL0636.65144 Graham I. G. Galerkin methods for second kind integral equations with singularities Mathematics of Computation 1982 39 160 519 533 10.2307/2007329 669644 ZBL0496.65068 Orav-Puurand K. Pedas A. Vainikko G. Nyström type methods for Fredholm integral equations with weak singularities Journal of Computational and Applied Mathematics 2010 234 9 2848 2858 10.1016/j.cam.2010.01.033 2652130 ZBL1193.65226 Huang J. T. Li Z. C. Mechanical quadrature methods and their splitting extrapolations for boundary integral equations of first kind on open arcs Applied Numerical Mathematics 2009 59 12 2908 2922 10.1016/j.apnum.2009.06.006 2560824 ZBL1175.65142 Rui Z. Jin H. Tao L. Mechanical quadrature methods and their splitting extrapolations for solving boundary integral equations of axisymmetric Laplace mixed boundary value problems Engineering Analysis with Boundary Elements 2006 30 5 391 398 2-s2.0-33646135783 10.1016/j.enganabound.2006.01.005 ZBL1187.65135 Gradshteyn I. S. Ryzhik I. M. Table of Integrals, Series, and Products 2007 7th San Diego, Calif, USA Academic Press xlviii+1171 2360010 Fogiel M. Handbook of Mathematical, Scientific, and Engineering 1994 Piscataway, NJ, USA Research and Education Association