CC BY
90
ORIGINAL RESEARCH ARTICLE DOI: https://doi.org/10.18599/grs.2018.2.108-114
Hydraulic fracturing efficiency evaluation in the vicinity of a single well for a reservoir with two fractures
Yu.I. Yakhina
Termosim LLC, Kazan, Russian Federation E-mail: [email protected]
Abstract. The solution of fluid flow problem in an unbounded homogeneous reservoir penetrated by a single well with two intersecting fractures with taking into account their hydrodynamic resistance is constructed and investigated. A general integral representation of the perturbed pressure field is obtained using the instantaneous point-source. As a particular case, the quasi-stationary operating mode of the well is considered. The accuracy and limitations of the obtained asymptotic solutions are estimated. A comparative analysis obtained results is done. Well productivity and the pseudoskin factor for the fractures system are determined, explicit analytical expressions for these characteristics are constructed. In the course of computational experiments, the interaction of intersecting fractures at different opening angles are investigated. An estimate of the efficiency of repeated hydraulic fracturing of the productive formation is obtained. It is shown that the maximum flow is achieved for a perpendicular arrangement of the fractures, and the distribution of the outflow (inflow) along the flat vertical fracture essentially depends on its relative filtration resistance.
Keywords: hydraulic fracturing, crack azimuth of repeated fracturing, efficiency of repeated fracturing, quasi-stationary operation of well, pseudoskin factor of repeated fracturing, filter resistance in cracks, instantaneous point-source, fracture hydrodynamics
Recommended citation: Yakhina Yu.I. (2018). Hydraulic fracturing efficiency evaluation in the vicinity of a single well for a reservoir with two fractures. Georesursy = Georesources, 20(2), pp. 108-114. DOI: https:// doi.org/10.18599/grs.2018.2.108-114
Introduction
In order to involve undrained and weakly drained mineral resources in development, hydraulic fracturing is used. This technology contributes to the creation of a high-conductive fracture, which allows increasing the productivity of producing wells or injectivity of injection wells. According to expert assessments of hydrodynamic studies, the hydraulic fracturing of the reservoir is the most effective geological and technical measure that guarantees an increase in the efficiency of poorly permeable reservoirs development. One of the first studies of stationary single-phase inflow to a well with a single fracture of hydraulic fracturing is presented in the paper (Prats, 1961). Theoretical calculations show that the fracturing allows several times to increase the well production by activating the weakly draining zones, the volume of extraction or injection depending on the conductivity and fracture length (Cinco-Ley et al., 1978; Meehan et al., 1989).
A review of a number of other publications and an analysis of the efficiency of fracturing in a well with a single vertical fracture in a circular feed loop was carried
© 2018 Yu.I. Yakhina. Published by Georesursy LLC This is an open access article under the CC BY 4.0 license (https://creativecommons.org/licenses/by/4.0/)
out (Morozov, 2016) on the basis of various estimates of productivity and the generalized pseudoskin factor concept (Economides et al., 2002).
With a significant reduction in the production rates with respect to the initial regimes, a second fracturing is performed. At the same time, to ensure full coverage of the reservoir and introduction of new reserves, a method of reorientation of the azimuth is used, which makes it possible to realize the development of the fracture in a direction different from the first hydraulic fracturing.
Repeated hydraulic fracturing is one of the most common methods of intensifying the extraction of mineral reserves and increasing reservoir productivity, which makes the task of analyzing and predicting its effectiveness topical. In a particular case, the problem of a stationary inflow to a well with several vertical fractures with uniformly distributed inflow along the fractures in the region bounded by the circular feed circuit was considered in the work (Raghavan, Joshi, 1993). A similar problem of a steady uniformly distributed inflow to two fractures with different geometric characteristics at a given intersection angle was solved numerically in (Lihtarev, Pestrikov, 2010).
In the framework of this study, theoretical issues related to the analysis of the unsteady hydrodynamic regime of one and two intersecting fractures of different
SCIENTl FK AND TECHNICAL JOURNAL
M GEDRESDURCES www.geors.ru
azimuth are considered taking into account their hydraulic resistance in an unbounded homogeneous reservoir. The efficiency of re-fracturing at different opening angles was estimated. Explicit analytical expressions are obtained for the productivity of the well and the pseudoskin factor.
General integral representation of the disturbance in the reservoir pressure field
During the operation of the well, fluid movement in the near-wellbore zone occurs due to the potential energy of the reservoir elastic state, as a result of lowering pressure at the bottom of the well. At the same time, the volume of liquid increases, and the porosity, and possibly the permeability, decreases due to the expansion of the formation skeleton. In the case of fluid injection in the formation, the flow of injected fluid is maintained by increased pressure at the bottom, and the filtration process develops in the opposite direction.
The main characteristic of filtration processes is the field of reservoir pressure, which, in the case of unsteady elastic filtration in a porous medium with distributed sources, is described by the piezoconductivity equation (Charnyy, 1963):
dP _ q*
dt mQp '
1*1 \y\M <°o,t>0,
where k = £/(^m0) is the conductivity coefficient, which determines the propagation velocity of pressure perturbations P(x,y,z,t) in a porous medium with permeability k, porosity m0, elastic capacity coefficient p, liquid viscosity y; q*(x,y,z,t) is the density of distributed sources, which is the amount of fluid entering the volume unit of the medium per unit time t; A is the differential Laplace operator in a system of Cartesian coordinates (x,y,z).
Let there be no perturbation of pressure in the reservoir at the initial moment of time, then the solution of the piezoconductivity equation for zero initial conditions can be written in a general integral form using the method of instant point sources (Carslaw, Jaeger, 1964; Tikhonov, Samarskii, 1999): t +00
P(X,y,Z,t) = j jfj
0 -oo
m0p
— t)cLÇ dr] dÇ dr,
7], Ç, x,y,z,t 1
(1)
where G*(^,n,Z,x,y,z,t - t) is the influence function of the instantaneous point source, which is the pressure perturbation at the point (x,y,z) at time t, caused by a single instantaneous point source at the point at time t = t:
GXÇ.rj, ^x,y,z,t-t) =
2^/îix(t—T)
The total inflow or outflow of fluid Q in a well with fractures per unit time is, by definition, equal to:
Q(t) = Jff q*(f,rj,(,t) dÇ dr] d(.
We use the obtained representation (1) to study the filtration process when liquid is injected in a homogeneous isotropic reservoir with a thickness h in the vicinity of an injection well with two flat vertical fractures created as a result of hydraulic fracturing and intersecting at an angle a (Fig. 1). In order to simplify further reasoning, we assume that fractures of length 2x/ and width Sf are symmetrically disposed about the axis of the well, have the same geometry and identical hydrodynamic properties. The flow rate Q(t) pumped into the well is distributed evenly along the thickness of the reservoir (fracture height) with the surface density q(x,t). Here, q is the amount of liquid injected per unit of formation thickness from a unit of fracture length per unit time. The density of the distributed source along the fracture, taking into account the smallness of the opening of the discontinuity, d is determined by the relation: q*(x,y,z,t) ~ q(x,t)/ Sf (at q = 0 outside the fracture |y| > df /2).
In view of the linearity of the problem under consideration and the identity of the fractures , the perturbation of the pressure field created by them can be regarded as the sum of the perturbations caused by single fractures , one of which is located along the x axis as shown in Fig. 1, and the other, a repeated fracture , is directed along the axis with coordinates:
x'\ /cos a — sin a\ /x\ v'I Vsina cosa/\y/-
— sin I cos a / \yj
Assuming homogeneity of the reservoir and uniformity of vertical injection, the general integral representation of the pressure field (1) after integration over n is transformed to the form:
0 -Xf
(*-fl2+y2 . e 4x(t-r) +e
(*'-?) V2
4x(t-T)
(*-f)2+(y-q)2+(z-P2
g 4X(t-T) .
(2)
Fig. 1. Injection well with suppressing flat fractures
SCIENTIFIC AND TECHNICAL JOURNAL
GEORESURSY M
The total flow rate Q takes the form: xf
Q(t) = 2h J q(x,t)dx. ~xf
Single fracture hydrodynamics
Hydraulic fracturing of the productive formation allows creating a fracture with a high permeability kf >> k, which leads to a decrease in the hydraulic resistance of the bottomhole zone and, as a consequence, to an increase in the filtration surface. For the same reason, the transient processes of pressure redistribution in the Pf fracture occur very rapidly, and the fluid flow at each instant of time can be described by a stationary filtration equation (Meehan et al., 1989):
kfSf d2Pf u dx2
= q; 0 < x < Xf, t > 0.
(3)
By assumption, the homogeneous fracture consists of two differently directed «wings» symmetrical with respect to the well, which allows us to consider the liquid filtration region in only one direction of the fracturing (for example, in the interval 0 < x < xf), as shown in Fig. 2.
At the bottom of the well, there is a perturbation of the fluid pressure Pw, relative to the reservoir, and the inflow from the end is considered negligible. Thus, equation (3) is supplemented by boundary conditions:
_ dPf Pf(.X,t)\x=0 — Pw,
= 0.
X-Xf
The solution of the problem formulated above, which describes the process of filtration in a fracture, can be represented in the form:
X xf
Jfq(f,t)d{ + xJ q((,t)df
V-
kf5f
(4)
Relation (4) relates the pressure distribution along the
0 < x < xf.
fracture to the liquid flow density q(x,t), and should be considered together with equation (2) at y = 0: t xf
qtf.r)
0 -Xf
iïzSl
(x'-flV 2
e 4x(t-T)+e «ct-T)
df dz,
(5)
where 0 < x < xf, t > 0, x' = x cosa and y' = x sina.
To determine the distributed inflow along the fracture, combining equations (4) and (5), we obtain the integral equation:
4 %(f,T) -iizfli fr'-flV2'
-JLf [ 3ÏL ' 4nk) J (t-
0 -Xf
T)
JJ_ kfSf
e 4K(t-T)+e 4«(t-T)
xf
J fq((,t)df + x J q(f,t)df
d^dr+
(6)
where 0 < x < xf and t > 0.
Quasistationary well regime with intersecting fractures
As preliminary estimates show, under unchanged operating conditions of the well (constant bottomhole pressure Pw or constant flow rate Q), the quasistationary filtration regime in the reservoir, characterized by the time-invariant profile of the relative distribution of the inflow along the fracture
q(x) = q(x, t)/< q > ,
where <q>- Q(t)/Qhxfi is set at xt/xf > 100.
In this case, the integral equation (6) takes the form:
Xf
+ Ei
+ Ei\-
+
4xt
(x'-02 + y
4xt
4xt
I
kfSf
df +
xf
J fq(f)df+x J q(0 df
(7)
Fig. 2. Geometrical model of a wing of a fracture
uu
f e~u
Here Ei(—x) = — -du - is integral exponential
J u
X
function, for small values of the argument (Abramovitz, Stegun, 1979):
Ei(-x) « ln(1.781%).
Accordingly, in the quasi-stationary approximation, instead of (7) for large values of time t, we have:
mOCO
Pw ~ 4hx
+ In
+ In
Mo x 7
((x-Q2\ ((x' + 02+y'2\ \2.2458 xt) n\ 2.2458 xt )
qx' + 02 + y'2\ Kx'-Q2+y'2\ V 2.2458 xt ) n\ 2.2458 xt )
kfSf
xf
J fq(Odf + x J q(0 df
df +
(8)
o x
where x' = x cosa h y' = x sina.
The numerical solution of the distributed flow rate of liquid along a fracture in the presence of a filtration resistance
Here and below, the following typical values for the injection conditions of the liquid (water) are accepted in the calculations: the pressure at the bottom is Pw = 4.053 106 Pa (40 atm), the fracture half-length is xf = 100 m, the fracture width is Sf = 0.2 m, the piezoconductivity coefficient is x = 2.78 m2/s (104 m2/ hour), the permeability of the reservoir is k = 10-13 m2 (100 mD), the viscosity of the fluid is y = 0.001 Pas (1 cps), the operating time of the well is t = 360000 sec (100 hours).
The results of the computational experiments presented in Fig. 3a have shown that the distribution of the inflow along the fracture essentially depends on its relative filtration resistance, kxf / (kf Sf). At the end of the fracture , a sharp increase in the outflow is observed, which is associated with a change in the character of the liquid flowfrom the linear to the quasiradial one.
In turn, for a known injection density q(x,t), which takes into account the filtration resistance in the fractures, the quasi-stationary pressure distribution in the reservoir can be calculated on the basis of the asymptotic representation of the general equation (2):
0
V 4>it J \ 4xt J
V 4*t /
where x' = x cosa - y sina and y' = x sina + y cosa.
A comparative analysis (Fig. 4) shows that the pressure drop along the fractures increases substantially with increasing filtration resistance.
Let us consider the distribution of the inflow into the reservoir with a uniform pressure distribution along the fractures , Pf = Pw, in the absence of filtration resistance (kf ^ ®). An analysis of the solution of the inverse problem of the outflow distribution at a given pressure shows that when the opening angle between fractures increases, the pressure gradients in the near wellbore zone of the reservoir decrease and the local outflow of liquid decreases, as shown in Fig. 5. And at the remote end part of the fracture , the reverse effect is observed - the flow density increases due to the decrease in interference between the fractures.
When the opening angle tends to 0°, the pressure in the well near the wellbore is equalized and the injection density is correspondingly reduced - two fractures begin to work as one and the outflow from the central part increases.
5<H
40 -
| 3<H
E
20-
10-
a)
10 20 30 -to 50 60 70 SO 90 X, in
3 534E
« 32-
b)
10 20 3» 40 50 60
x, m
70 SO 90
Fig. 3. The distribution of flows to the fracture (a) and the pressure along the fracture (b) in the vicinity of the well at a mutually perpendicular arrangement of fractures and different permeabilities: 1 - kf = 5104 mD and 2 - kf = 104 mD.
SCIENTIFIC AND TECHNICAL JOURNAL
www.geors.ru GEDRESURSY
40—I
-200
q)
100 100 200
40-5Î-
E ,
S3
CC 25-1 2015-1
-100 ,) b) y, m
200
200
x, m
Fig. 4. Perturbation of reservoir pressure in the vicinity of the well with mutually perpendicular arrangement offractures with permeability: (a) - kf = 5104 mD (kx / (kd) ~ 1) and (b) - k = 104 mD (kx / (kd) ~ 5)
100 so 80 70
q)
10 20 30 40 30 60
x, m
b)
Fig. 5. Distribution of the inflow along the primary fracture (a) at a constant depression of 40 atm at different opening angles and perturbation of the reservoir pressure (b) in the vicinity of the well with mutually perpendicular arrangement of the fractures
Well productivity and pseudoskin factor
Summarizing the traditional concept of the skin factor (Economides et al., 2002; Morozov, 2016), we define the pseudoskin factor of the well with intersecting fractures of the finite permeability for a given depression. From the relation (8) for x = 0 it follows immediately:
-Xf
Pw = ~
fiQCt)
4 hxf
IMlln(T=L=W
J nk VV2.2458*t/
2nkhxf
fxf
f m J In ^
\o
Itl I
2.2458*t\ \
rà J dt /
where r is the radius of the well. As a result, the
w '
productivity of the well will be given by the formula: Q 2nhk 1
S + ln
(f
2458xt^
(9)
where pseudoskin factor is:
xf
u
In the particular case of a single fracture of infinite conductivity (jq = 1) the pseudoskin factor in (9) is equal to S = - ln(Xf/e). In addition, according to the estimates obtained in (Charnyy, 1963), on the basis of the method of successive change of stationary states, the expression V2.2458>rt can be considered as the radius of the mobile current circular feed loop re. As a result, the above relation (9) can be rewritten in the equivalent form:
vQ
Pw =
2nkh
In
(Jl.)
which coincides with the solution (Raghavan, Joshi, 1993).
Calculation results of pseudoskin factor (9), shown in Fig. 6, show that the minimum pseudoskin factor and the maximum injection (inflow) volume increase is achieved with the perpendicular arrangement of fractures and essentially depends on their resistance to filtration.
X , 111
10 20 30 40 50 60 70 80
-6
a)
90
b)
Xf, 111
Fig. 6. Pseudoskin factor at a given bottomhole pressure Pw with permeability (a) kf= 106 mD (kxf / (kf Sf) ~ 0,05) and (b) f = 104 mD (kxf / (kf S) ~ 5)
Repeated hydraulic fracturing allows to increase the pseudoskin factor at high fracture permeability up to 5-10%, as shown in Fig. 6a. Efficiency of re-fracturing according to Fig. 6b increases to 10-20% in the case of small fractures permeability, which is in agreement with the numerical solutions obtained in (Lihtarev, Pestrikov, 2010).
Conclusion
As is known, hydraulic fracturing of the reservoir is a complex, energy-intensive and expensive technological process for intensifying field development. Every year, this technology is increasingly used in the operation of both producing and injection wells. However, the consequences of applying repeated fracturing are not always positive, which makes the task of assessing its effectiveness urgent.
In the study of the hydrodynamic regime of a well with intersecting fractures penetrating a homogeneous reservoir, it is shown that under unchanged conditions of well operation (steady pressure at the bottom Pw or constant flow rate Q), a quasi-stationary operation mode of the well with a constant injection profile (inflow) in fractures is established. The distribution of the inflow along the fracture essentially depends on its relative filtration resistance, kxf / (kf Sf ). At the same time, as the opening angle is increased, the pressure gradients in the near wellbore zone of the reservoir decrease and the local outflow (inflow) of the liquid decreases. At the remote end part of the fractures , the reverse effect is observed - the flux density increases due to the decrease in interference between the fractures. The maximum total increase in inflow into the reservoir is achieved with a perpendicular arrangement of the fractures .
Repeated hydraulic fracturing of the productive formation allows to increase the pseudoskine factor
at a high permeability (k
106 mD, kx /(k S ) •
' 0,05)
fractures up to 5-10%. The effectiveness of re-fracturing increases to 10-20% in the case of small permeability (k~ 104 mD, kx /(k S ) ~5).
Further research continuation can be aimed at solving both more general problems of unsteady filtration and inverse problems: justifying the methods of hydrodynamic exploration of wells in order to determine the opening angle between two hydraulic fracturing fractures and other hydrodynamic characteristics of the fracture system.
Acknowledgements
The author expresses gratitude to the scientific adviser, Professor A.N. Salamatin for valuable advice in the planning and implementation of research and recommendations on the formatting of work. The author thanks Professor G.T. Bulgakova for choosing the topic and area of research and expresses gratitude to the reviewer for valuable comments and suggestions that contributed to the improvement of this article.
References
Abramowitz M., Stegun A. (1979). Handbook of Mathematical Functions. Moscow: Nauka, pp. 58-59. (In Russ.)
Carslaw H.S, Jaeger J.C. (1964). Conduction of Heat in Solids. Moscow: Nauka, pp. 251-259. (In Russ.)
Charnyy I. A. (1963). Underground hydro gas dynamics. Moscow: Gostoptekhizdat, pp. 250-313. (In Russ.)
Cinco-Ley H., Samaniego V.F., Domínguez A.N. (1978). Transient pressure behavior for a well with a finite-conductivity vertical fracture. SPE Journal, 18(4), pp. 253-264.
Economides M., Oligney R., Valkó P. (2002). Unified fracture design: bridging the gap between theory and practice. Alvin, Texas: Orsa Press, 25 p.
Lihtarev A. V., Pestrikov A. V. (2010). Opredelenie kojefficienta produktivnosti skvazhiny, peresechennou dvumua treshhinami gidrorazryva plasta raznogo azimuta, na osnove matematicheskogo modelirovaniua [Determination of the production efficiency of a well crossed by two induced hydraulic fractures of a different azimuth, based on mathematical modeling]. Nauchno-tehnicheskiy vestnikOAO "NK"Rosneft'", pp. 12-14. (In Russ.)
SCIENTIFIC AND TECHNICAL JOURNAL
www.geors.ru GEDRESURSY
Meehan D. N., Horne R.N., Ramey H.J. (1989). Interference testing of finite conductivity hydraulically fractured wells. SPEAnnual Technical Conference and Exhibition. Society of Petroleum Engineers, pp. 137-152.
Morozov P.E. (2016). Psevdoskin-faktor i optimal'naya provodimost' vertikal'noy treshhiny gidavlicheskogo razryva plasta [Pseudoskine factor and optimal conductivity of vertical induced hydraulic fracture]. Mezhd. nauchno-prakt. konf.: Innovacii v razvedke i razrabotke neftyanyh i gazovyh mestorozhdeniy [Proc. Sci. and Pract. Conf.: Innovations in explorationa and development of oil and gas fields], Kazan, pp. 53-56. (In Russ.)
Prats M. (1961). Effect of vertical fractures on reservoir behavior-incompressible fluid case. SPE Journal, 1(02), pp. 105-118.
Raghavan R., Joshi S.D. (1993). Productivity of multiple drainholes or fractured horizontal wells. SPE Form. Eval., 8(01), pp. 11-16.
Tikhonov A. N., Samarskii A. A. (1999). Equations of Mathematical Physics. Moscow: MSU, pp. 477-485. (In Russ.)
About the Author
Yulia I. Yakhina - Junior Software Engineer, Engineering-Analytical Group, Department of Scientific Research, Termosim LLC
59/1 Magistralnaya str., Kazan, 420108, Russian Federation
E-mail: [email protected]
Manuscript received 29 January 2018;
Accepted 28 April 2018; Published 30 June 2018
